shiyan_1.m

一些matlab底下的小函数

function shiyan_1()
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);
flag=1;
while(flag==1)
    step=input('请选择你要进行的实验步骤:\n 1(时域采样序列分析),2(系统和响应分析),3(卷积定理的验证),0(退出)\n');
    while step<0&step>3
       step=input('请选择你要进行的实验步骤:\n 1(时域采样序列分析),2(系统和响应分析),3(卷积定理的验证),0(退出)\n');
    end
    switch step
        case 0,
            flag=0;
        case 1,         
%实验内容1
A=444.128;
a=50*sqrt(2)*pi;
w=50*sqrt(2)*pi;
xa=FF(A,a,w,1000);
DFT(xa,50,'xa(n)','|Xa(jw)|','fs=1000');
xa=FF(A,a,w,500);
figure;
DFT(xa,50,'xa(n)','|Xa(jw)|','fs=500');
xa=FF(A,a,w,200);
figure;
DFT(xa,50,'xa(n)','|Xa(jw)|','fs=200');
input('请按《Enter》键以结束步骤1');
close all;
          case 2,     
%实验内容2
xc1=stepseq(1,1,10);
xc2=stepseq(1,1,5);
DFT(xb,2,'xb(n)','|Xb(jw)|','');
figure;
[Hb,wb]=DFT(hb,4,'hb(n)','|Hb(jw)|','');
yb=conv(xb,hb);
figure;
DFT(yb,length(yb),'yb(n)','|Yb(jw)|','yb=xb*hb');
yc1=conv(xc1,ha);
figure;
DFT(yc1,length(yc1),'yc1(n)','|Yc1(jw)|','yc1=xc1*ha,xc的长度为10');
yc2=conv(xc2,ha);
figure;
DFT(yc2,length(yc2),'yc2(n)','|Yc2(jw)|','yc2=xc2*ha,xc的长度为5');
input('请按《Enter》键以结束步骤2');
close all;
        case 3,
%实验内容3
A2=1;
a2=0.4;
w2=2.0734;
xa2=FF(A2,a2,w2,1);
[Xa2,wa2]=DFT(xa2,length(xa2),'xa2(n)','|Xa2(jw)|','A=1,a=0.4,w=2.0734,T=1');
figure;
[Hb,wb]=DFT(hb,4,'hb(n)','|Hb(jw)|','');
ya=conv(xa2,hb);
figure;
DFT(ya,length(ya),'ya(n)','|Ya(jw)|','ya=xa2*hb');
%画Xa2·Hb＝FT[xa2]·FT[hb]的图形
Ya2=abs(Xa2.*Hb);
figure;
subplot(1,2,2);
plot(wa2/pi,Ya2);
xlabel('w/pi');
ylabel('Ya2');
title('Ya2＝FT[xa2]·FT[hb]');
input('请按《Enter》键以结束步骤3');
close all;
    end
end