3、in=f(t)%輸入信號in=(t>0)*1;%階躍信號%filenameex123.mL=2;%電感值C=0.32;%電容值forR=[1.459]%仿真電阻值分別為1.5,3,5歐姆的情況[t,x]=ode45('funcforexl23',[0,15],[0;0.5],[],R,L,C);%也可采用ode23,ode15s等求解figure(1);plot(t,x(:,1));holdon;xlabel('timesec');text(2,0.07,'leftarrowi_L(t)');grid;figure(2);plot(t,x(:,2))
4、;holdon;xlabel('timesec');text(2.1,0.75,'leftarrowu_C(t)');grid;End輸入輸出的傳遞函數(shù):① R=1.4時:② R=5時:③ R=9時:連續(xù)系統(tǒng)的傳遞函數(shù)如下:借助多項(xiàng)式乘法函數(shù)conv來處理:兩個向量分別用num和den表示。num=[b1,b2,…,bm,bm+1]den=[a1,a2,…,an,an+1]num=4*conv([1,2],conv([1,6,6],[1,6,6]));den=conv([1,0],conv([1,1],conv([1,1],conv([1,1],[1
5、,3,2,5]))));給出零極點(diǎn)增益模型:num=[1,11,30,0];den=[1,9,45,87,50];[z,p,k]=tf2zp(num,den)結(jié)果表達(dá)式:用Simulink系統(tǒng)建模仿真畫出階躍信號通過傳遞函數(shù)后的關(guān)系曲線。還可以進(jìn)一步修改信號源參數(shù),使用三角波、正弦波等作為激勵信號,觀察輸出信號的情況。<1>階躍信號<2>三角波<3>正弦波實(shí)驗(yàn)作業(yè)2第一題WaterfallScope(瀑布顯示圖)1、利用WaterfallScope(瀑布顯示圖)的仿真模型,觀察ChirpSignal掃頻信號源的頻譜時間函數(shù),顯示需要表達(dá)的三維函數(shù)瀑布圖
6、。更換掃頻信號源,選擇UniformRandomNumber或Band-LimitedWhiteNoise信號源,的頻譜時間函數(shù),顯示需要表達(dá)的三維函數(shù)瀑布圖。ChirpSignal掃頻信號源應(yīng)用WaterfallScope觀察掃頻儀輸出下圖所示是仿真系統(tǒng)顯示的瀑布圖。系統(tǒng)運(yùn)行時瀑布圖中的頻譜線群隨著掃頻儀的頻率變化而移動,十分生動直觀。表1表3所示是原理圖中的各模塊的主要參數(shù)。UniformRandomNumber信號源下結(jié)果圖Band-LimitedWhiteNoise信號源第二題:設(shè)計(jì)一個編程開關(guān)仿真系統(tǒng)框圖設(shè)計(jì)一個編程開關(guān)仿真系統(tǒng)框圖。圖中使用
7、三個信號產(chǎn)生器,分別產(chǎn)生幅度為1,頻率為10Hz、20Hz、30Hz的正弦波、鋸齒波和方波。VariableSelector(可變選擇器)根據(jù)Idx(指針端口)當(dāng)時的數(shù)據(jù),決定三個信號饋送到哪個端口。程序:x1=[123231312];>>x=[x1]'*[ones(1,30)];%生成一個9行30列的數(shù)據(jù)塊>>y=x(:);%變換成270行1列的數(shù)據(jù)塊>>y1=reshape(y,3,90);%變換成3行90列的數(shù)據(jù)塊>>y2=y1';%變換成90行3列的數(shù)據(jù)塊>>t1=1:90;t=t1';>>D=[ty2];%產(chǎn)生D供FromWorkspace讀
8、入>>sim('switchil')%啟動仿真編程開關(guān)仿真系統(tǒng)框圖:運(yùn)行結(jié)果:仿真實(shí)驗(yàn)作業(yè)3第