untitled.m

用SVD奇异值分解法和阻尼最小二乘算法解决地球物理反问题

x=0:0.02:1;
y=0:1:10;
n=length(x);
m=length(y);
G=zeros(m,n);
pi=3.14;
for i=1:m
    for j=1:n
        G(i,j)=exp(-i*x(j));
    end
end
M=zeros(n,1);
for i=1:n
    M(i,1)=1-0.5*cos(2*pi*x(i));
end
[U,S,V]=svd(G);
 D=V*pinv(S)*U'*M; 
 y=o*randn(M,N)+u