信号与系统陈后金matlab,陈后金信号与系统matlab实验4.pdf

陈后金信号与系统matlab实验4

第五章非周期信号的频域分析

作者卢未来 日期 年 月 日

: :2012 5 10

5.6非周期信号的频域分析的MATLAB实现

5-19(174 )

例 页 试用数值方法近似计算三角波信号

f(t)= (1−|t|)p (t)的频谱。

2

编写函数

function y=sf1(t,w);

y=(t>=-1&t<=1).*(1-abs(t)).*exp(-j*w*t);

执行代码

%program5_1

w=linspace(-6*pi,6*pi,512);

N=length(w);F=zeros(1,N);

for k=1:N

F(k)=quad('sf1',-1,1,[],[],w(k));

end

figure(1);

plot(w,real(F));

xlabel('\omega');

ylabel('F(j\omega)');

figure(2);

plot(w,real(F)-sinc(w/2/pi).^2);

xlabel('\omega');

title('计算误差');

5-20(175 ) 1

例 页 试计算宽度和幅度均为 的方波信号p(t)在0~ f (Hz)频

1 m

谱范围内所包含的信号能量。

编写函数:

function y=sf2(t);

y=2*sinc(t).*sinc(t);

执行代码:

>>%program5_2

>>f=linspace(0,5,256);

>>N=length(f);w=zeros(1,N);

>>fork=1:N

w(k)=quad('sf2',0,f(k));

end

>>plot(f,w);

>>xlabel('Hz');

>>ylabel('E');

例5-21(176页)试画出α = ±0.9时

1

jΩ

F(e )= − jΩ

1−αe

的幅度频谱。

%program5_3

b=[1];

a1=[1-0.9];a2=[10.9];

w=linspace(0,2*pi,512);

h1=freqz(b,a1,w);

h2=freqz(b,a2,w);

plot(w/pi,abs(h1),w/pi,abs(h2),':');

legend('\alpha=0.9','\alpha=-0.9');