51qiang168
2011-12-19, 20:11
现在假设有一地震信号是是s(n)=exp(-0.0016*n.^2)*sin(n*pi*0.2)+exp(-0.0016*n.^2)*sin(n*pi*0.32)，其n只能取整数，总长度为200第一步要先产生地震信号s(t)，让后让他与一滤波因子h(t)=pi*t/2cos(w0)sin(w0*t)(假设)，利用褶积求和实现对地震信号的滤波。我在产生地震信号时编写的程序如下：
x=[];
for n=0:199
temp=exp(-0.0016*n.^2)*sin(n*pi*0.2)+exp(-0.0016*n.^2)*sin(n*pi*0.32);
x=[x temp];{在这里系统提示说有点错，但是不影响程序的运行｝
end
n=0:199;
plot(n,x);
hold on
图是能画出来，但是不知道对不？我总觉得有错，我就想请教怎样利用循环产生这个地震信号？
x=[];
for n=0:199
temp=exp(-0.0016*n.^2)*sin(n*pi*0.2)+exp(-0.0016*n.^2)*sin(n*pi*0.32);
x=[x temp];{在这里系统提示说有点错，但是不影响程序的运行｝
end
n=0:199;
plot(n,x);
hold on
图是能画出来，但是不知道对不？我总觉得有错，我就想请教怎样利用循环产生这个地震信号？