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('請(qǐng)選擇你要進(jìn)行的實(shí)驗(yàn)步驟:\n 1(時(shí)域采樣序列分析),2(系統(tǒng)和響應(yīng)分析),3(卷積定理的驗(yàn)證),0(退出)\n');
    while step<0&step>3
       step=input('請(qǐng)選擇你要進(jìn)行的實(shí)驗(yàn)步驟:\n 1(時(shí)域采樣序列分析),2(系統(tǒng)和響應(yīng)分析),3(卷積定理的驗(yàn)證),0(退出)\n');
    end
    switch step
        case 0,
            flag=0;
        case 1,         
%實(shí)驗(yàn)內(nèi)容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('請(qǐng)按《Enter》鍵以結(jié)束步驟1');
close all;
          case 2,     
%實(shí)驗(yàn)內(nèi)容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的長(zhǎng)度為10');
yc2=conv(xc2,ha);
figure;
DFT(yc2,length(yc2),'yc2(n)','|Yc2(jw)|','yc2=xc2*ha,xc的長(zhǎng)度為5');
input('請(qǐng)按《Enter》鍵以結(jié)束步驟2');
close all;
        case 3,
%實(shí)驗(yàn)內(nèi)容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');
%畫(huà)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('請(qǐng)按《Enter》鍵以結(jié)束步驟3');
close all;
    end
end