DOC《实验报告》

close all;

Xb=impseq(0,0,1);

Ha=stepseq(1,1,10);

Hb=impseq(0,0,3)+2.5*impseq(1,0,3)+2.5*impseq(2,0,3)+impseq(3,0,3);

i=0;

while(s);

if(s==1) l=1;

while(l) A=input('请输入振幅:\n[100,1000]');

a=input('请输入衰减:\n[100,600]');

w=input('请输入角频率(rad/s):\n[100,600]');

fs=input('请输入采样频率:\n[100,1200]');

Xa=FF(A,a,w,fs);

i=i+1;

string=['fs=',num2str(fs)];

figure(i) DFT(Xa,50,string);

l=input('是否继续实验:1(继续实验1),0(进行其他实验内容):\n');

end else if(s==2) kk=input('系统和系统响应分析,请选择时域信号的类型:1(内容2 a),2(内容2 b),3(2 b中xc(n)的长度改为5),0(退出):\n');

while(kk) if(kk==1) m=conv(Xb,Hb);

N=5;

i=i+1;

figure(i) string=('hb(n)');

Hs=DFT(Hb,4,string);

i=i+1;

figure(i) string=('xb(n)');

DFT(Xb,2,string);

string=('y(n)=xb(n)*hb(n)');

else if(kk==2) m=conv(Ha,Ha);

N=19;

string=('y(n)=ha(n)*ha(n)');

else if(kk==3) Xc=stepseq(1,1,5);

m=conv(Xc,Ha);

N=14;

string=('y(n)=xc(n)*ha(n)');

end end end i=i+1;

figure(i) DFT(m,N,string);

kk=input('请再选择信号类型:1,2,3,0(进行其他实验内容):\n');

end else if(s==3) A=1;

a=0.5;

w=2.0734;

fs=1;

Xa1=FF(A,a,w,fs);

i=i+1;

figure(i) string=('The xa1(n)(A=1,a=0.4,T=1)');

[Xs,w]=DFT(Xa1,50,string);

i=i+1;

figure(i) string=('hb(n)');

Hs=DFT(Hb,4,string);

Ys=Xs.*Hs;

y=conv(Xa1,Hb);

N=53;

i=i+1;

figure(i) string=('y(n)=xa(n)*hb(n)');

[yy,w]=DFT(y,N,string);

i=i+1;

figure(i) subplot(2,2,1) plot(w/pi,abs(yy));

axis([-2

2 0 2]);

xlabel('w/pi');

ylabel('|Y(jw)|');

title('FT[x(n)*h(n)]');

subplot(2,2,3);

plot(w/pi,abs(yy));

axis([-2

2 0 2]);

xlabel('w/pi');

ylabel('|Ys(jw)|');

title('FT[xs(n)].FT[h(n)]');

end end end s=input('请选择实验步骤:\n 1(1 时域采样序列分析),2(2 系统和系统响应分析),3(3 卷积定理的验证),0(退出):\n');

end 相应子程序function [c,l]=DFT(x,N,str) n=0:N-1;

k=-200:200;

w=(pi/100)*k;

l=w;

c=x*(exp(-j*pi/100)).^(n'*k);

magX=abs(c);

angX=angle(c);

subplot(2,2,1);

t=max(x);

n=0:N-1;

stem(n,x,'');

xlabel('n');

ylabel('xa(n)');

title('信号的原形');

text((0.3*N),(0.8*t),str);

hold on n=0:N-1;

m=zeros(N);

plot(n,m);

subplot(2,2,3);

plot(w/pi,magX);

xlabel('w/pi');

ylabel('|Y(jw)|');

title('上图的信号傅氏变换(|jw|)');

function c=FF(A,a,w,fs) n=0:50-1;

c=A*exp((-a)*n/fs).*sin(w*n/fs).*stepseq(0,0,49);

function[x,n]=impseq(n0,n1,n2) n=[n1:n2];

x=[(n-n0)==0];

function[x,n]=stepseq(n0,n1,n2) n=[n1:n2];

x=[(n-n0>=0)];

七 程序结果及

图表: : 系统和系统响应分析实验图形如下 八 实验总结:在分析理想采样序列特性的实验中,采样频率不同时,相应理想采样序列的傅立叶变换频谱的数字频率度量都相同,模拟频率也相同,因为只有采用定性的方法才能比较出来采样频率大于角频率的2倍时可以出现混叠.在卷积定理验证的实验中,如果选用不同的频域采样点数M值求得结果是有差异的 题目 (十一) 一:实验题目: 对给出的信号进行谱分析. x1(n)=R4(n);

n+1 , 0