ccint.ma

国外免费地震资料处理软件包

ccin[x_ /; Abs[x] <= 1]:= 3/2 Abs[x]^3 - 5/2 Abs[x]^2 + 1;
ccin[x_ /; Abs[x] > 1 && Abs[x] <= 2]:= 
	-1/2 Abs[x]^3 + 5/2 Abs[x]^2 - 4 Abs[x] + 2;
ccin[x_ /; Abs[x] > 2] := 0;
ccinf[w_]:= 2 Integrate[Cos[w x] (3/2 x^3 - 5/2 x^2 + 1),{x,0,1}] +
	    2 Integrate[Cos[w x] (-1/2 x^3 + 5/2 x^2 - 4 x + 2),{x,1,2}];
Plot[ccin[x], {x, -2, 2}, PlotStyle -> {Thickness[0.01]}, 
    	Frame -> True, 
    	FrameLabel -> {None, None, 
        "Cubic Convolution Interpolation", None}];
Plot[ccinf[w],{w,-2Pi,2Pi},PlotStyle->{Thickness[0.01]},
	Frame->True,
 	FrameLabel->{None,None,"Spectrum",None}];
Show[GraphicsArray[{%%,%}]];
Display["junk_ma.eps",%,"EPS",ImageSize->648];