怎样在matlab中输出雷克子波

雷克子波及其振幅谱的matlab输出：

fm=50%主频

dt=0.001%时间域采样间隔

number=100%采样点数

t=-number/2+1:number/2

a=(1-2*(pi*fm*t*0.001).^2).*exp(-(pi*fm*t*0.001).^2)

subplot(2,1,1)

plot(t,a)

title('Ricker-雷克子波')

xlabel('时间t（ms）')

ylabel('幅值A')

for i=1:100

f(i)=10*(i-1)%计算频率

end

Y=abs(fft(a))%fourier变换，取振幅谱

subplot(2,1,2)

plot(f,Y)

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

xlabel('频率f（hz）')

ylabel('振幅森亩胡谱')

%注耐旦：时间域采样间隔为0.001s，采样点数为100点，总的时间长度为0.1s，则频率域采样间隔此拦为10hz。

W=fft(w,N)

thisfilter(1:17)=0.01

thisfilter(18:499)=1

thisfilter(500:N)=0.01%%%W对应的50Hz的峰值最大，找到对应的两个最简桐大峰并咐搭值的位绝拿置

thispy3=thisfilter.*W%滤波

你把AM,FM的信号时域表达式直接用语句写出来不就成了

比如:

t=0:0.1:100

y=10cos(2.*pi.*t+cos(2.*pi.*t))

不就是裤春扰一个调频波

至于有胡旦没有直接的函森消数.我倒是没见到过

---