变参数的积分程序.txt

谐波小波程序变参数的积分程序多点人工生成地震波程序低通滤波程序经验模式分解降噪程序

function y =wjl 
a = -100:100;
coef = meshgrid(a,a);
[ii,jj]= size(coef);
sol = zeros(ii,jj);
for (i =1:ii)
for (j = 1:jj)
    sol(i,j)=quad(@fun,0,10,1e-7,[],coef(i,j))  ;
end
end
y = sol;
surf(sol)
shading interp

function [y] =fun(x,n)
y = besselj(n,x);