資源描述：
《信號(hào)與線性系統(tǒng)實(shí)驗(yàn)報(bào)告材料.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、課程實(shí)驗(yàn)報(bào)告課程名稱：信號(hào)與線性系統(tǒng)專業(yè)班級(jí)：計(jì)算機(jī)科學(xué)與技術(shù)11級(jí)10班學(xué)號(hào)：U201114445姓名：王濤指導(dǎo)教師：偉平報(bào)告日期：2013/7/1計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院實(shí)驗(yàn)一連續(xù)時(shí)間系統(tǒng)的時(shí)域分析%1_1_1.mT=0.01;tx=0:T:2;x=rectpuls((tx-1),2).*exp(-tx);th=0:T:2;h=rectpuls((th-1),2).*2;t=(0+0):T:(2+2);y=conv(x,h);figuresubplot(3,1,1);plot(tx,x)ylabel('輸入激勵(lì)');subplot(3,1,2);plot(t
2、h,h)ylabel('單位沖激響應(yīng)');subplot(3,1,3);plot(t,y)ylabel('輸出響應(yīng)');%1_1_2.mT=0.01;tx=-4:T:4;x=rectpuls((tx-0),8).*(1-abs(tx)/4);th=0:T:4;h=(th>=0);t=(-4):T:8;y=conv(x,h);figuresubplot(3,1,1);plot(tx,x)ylabel('輸入激勵(lì)');subplot(3,1,2);plot(th,h)ylabel('單位沖激響應(yīng)');subplot(3,1,3);plot(t,y)ylabel(
3、'輸出響應(yīng)');%1_2_1.mR=10;%電阻值C=4;%電容值U0=2;%電容電壓初始值S=[num2str(C),'*',num2str(R),'*Dy+','y=heaviside(t)'];%得到微分方程的字符串表達(dá)式init=['y(0)=',num2str(U0)];%得到初始條件的字符串表達(dá)式y(tǒng)=dsolve(S,init,'t');%求解微分方程，得到符號(hào)解t=2:0.01:2*pi;x=ones(1,length(t));figuresubplot(1,2,1);plot(t,x)title('激勵(lì)電壓');subplot(1,2,2);
4、%2號(hào)子圖顯示輸出響應(yīng)ysezplot(y,[0,2*pi])title('輸出響應(yīng)')%1_2_2.mR=0.5;%電阻值L=1;%電容值I0=2;%電感電流初始值S=[num2str(L),'*Dy+',num2str(R),'*y=heaviside(t)'];%得到微分方程的字符串表達(dá)式init=['y(0)=',num2str(I0)];%得到初始條件的字符串表達(dá)式y(tǒng)=dsolve(S,init,'t');%求解微分方程，得到符號(hào)解t=0:0.01:2*pi;x=ones(1,length(t));%得到時(shí)間圍t的階躍信號(hào)u(t)的離散抽樣序列xf
5、iguresubplot(1,2,1);%多子圖顯示，將圖形框分為1x2個(gè)子圖，1號(hào)子圖顯示階躍信號(hào)xplot(t,x)title('激勵(lì)電壓');subplot(1,2,2);%2號(hào)子圖顯示輸出響應(yīng)ysezplot(y,[0,2*pi])title('輸出響應(yīng)')實(shí)驗(yàn)分析：本次實(shí)驗(yàn)應(yīng)用連續(xù)系統(tǒng)時(shí)域分析的兩種方法，卷積計(jì)算和微分方程的求解。實(shí)驗(yàn)二信號(hào)的傅里葉分析實(shí)驗(yàn)%2_1_1.mN=10;%希望看到的諧波次數(shù)symstTk;%定義時(shí)間t、周期Ta、下標(biāo)k三個(gè)符號(hào)y=subs(sym('Heaviside(t+T/20)-Heaviside(t-T/20)'
6、),'T',TT);%鋸齒脈沖基礎(chǔ)波形，tao/T=1/20%修改上面的sym中的參數(shù)，將T/20改為T/40皆可%得到占空比為1/20的。A0=int(y,t,-TT/2,TT/2)/TT;%直流分量A0%利用符號(hào)法得到各頻率分量的復(fù)數(shù)Fourier系數(shù)Ak的符號(hào)表達(dá)式Ak=int(y*exp(-2*i*pi*k*t/TT),t,-TT/2,TT/2)/TT;fk=sym(Ak)*sym(exp(2*i*k*pi*t/TT));%計(jì)算k在區(qū)間[-N,N]的Fourier系數(shù)，存入向量a（長(zhǎng)度為2*N+1）form=-N:-1%a(m+N+1)=numeri
7、c(subs(Ak,k,m));a(m+N+1)=double(subs(Ak,k,m));enda(N+1)=double(sym(A0));form=1:Na(m+N+1)=double(subs(Ak,k,m));end%for%利用symsum命令進(jìn)行符號(hào)求和，得到綜合信號(hào)ff=symsum(fk,k,-N,-1)+A0+symsum(fk,k,1,N);%圖形顯示結(jié)果figuren=-N:N;as=abs(a)*2;%由Fourier系數(shù)得到頻譜幅度subplot(3,1,1)ezplot(y,[-TT,TT])ylabel('原函數(shù)');subp
8、lot(3,1,2)ezplot(f,[-TT,TT
