雷克子波的matlab编程

呵呵，你的for循环用错了。把t=-50:50提到for循环前。

clear all

close all

dt=5*exp(-12)

freq=10^5

t=-50:1:50

lt=length(t)

V0=1

for k=1:1:lt

v(k)=(1-2*((((4/freq)/dt+1)/2-t(k)+1)*freq*dt)^2)*exp(-((((4/freq)/dt+1)/2-t(k)+1)*freq*dt)^2)

end

plot(t,v)

title('Ricker-雷克子波')

xlabel('时间（ns）')

ylabel('幅值')

还有你的dt和freq不太匹配，不过程序你可以自己改变的。有图片

%% 自己写了如下程序，望采纳！clcclearclose all

%Ricker子波

f=40dt=0.004

nw=6./f/dt

nw=2*floor(nw/2)+1

nc=floor(nw/2)

i=1:nw

alpha=(nc-i+1).*f*dt*pi

beta=alpha.^2

ricker=(1.-beta.*2).*exp(-beta)

subplot(311)

plot(ricker)

title('Ricker子波')

t=1:size(ricker,2)%绘制x(t)的振幅谱和相位谱

k=-25:25

W=(pi/12.5)*k

Fv=ricker*(exp(-j*pi/12.5)).^(t'*k)

magFv=abs(Fv)%绘制振幅谱

subplot(312)

plot(magFv)

title('Ricker子波的振幅谱')

angFv=angle(Fv)

subplot(313)

plot(angFv)

title('Ricker子波的相位谱')

helpdlg('欢迎采纳！有问题请继续提问！谢谢','提醒')

---