實驗四MATLAB用于頻域分析實驗?zāi)康睦肕ATLAB可以畫出頻率特性的精確圖形求取系統(tǒng)的穩(wěn)定裕量對系統(tǒng)進(jìn)行分析和設(shè)計實驗原理對于系統(tǒng)的傳遞函數(shù)，通過函數(shù)調(diào)用可畫出2種頻率特性圖。繪制伯德圖可用函數(shù)Bode(num,den),該函數(shù)表示在同一幅圖中，分上、下兩局部生成對數(shù)幅頻特性和相頻特性，橫軸均以rad/sec為單位。繪制奈氏圖可用函數(shù)nyquist(num,den)。通過這兩種形式的調(diào)用，可以計算G(jw)的實部和虛部，但是不能直接在屏幕上產(chǎn)生奈氏圖，需要通過調(diào)用plot(re,im)函數(shù)才可得到奈氏圖。對系統(tǒng)進(jìn)行頻率特性分析時，相位裕量和幅值裕量是衡量系統(tǒng)相對穩(wěn)定性的重要指標(biāo)，應(yīng)用MATLAB函數(shù)可以方便的求出系統(tǒng)的相位裕量和幅值裕量。[gm,pm,wcg,wcp]=margin(mag,phase,w);實驗代碼1、bode([1],[10])nyquist([1],[10])2、%伯德圖num=[1];zeta1=0.1;den1=[12*zeta11];zeta3=0.3;den3=[12*zeta31];zeta5=0.5;den5=[12*zeta51];zeta7=0.7;den7=[12*zeta71];zeta9=0.9;den9=[12*zeta91];[mag1,phase1,w1]=bode(num,den1);[mag3,phase3,w3]=bode(num,den3);[mag5,phase5,w5]=bode(num,den5);[mag7,phase7,w7]=bode(num,den7);[mag9,phase9,w9]=bode(num,den9);subplot(211);semilogx(w1,20*log10(mag1),w3,20*log10(mag3),w5,20*log10(mag5),w7,20*log10(mag7),w9,20*log10(mag9));subplot(212);semilogx(w1,phase1,w3,phase3,w5,phase5,w7,phase7,w9,phase9);%奈氏圖num=[1];zeta4=0.4;den4=[12*zeta41];zeta6=0.6;den6=[12*zeta61];zeta8=0.8;den8=[12*zeta81];[re1,im1]=nyquist(num,den4);[re2,im2]=nyquist(num,den6);[re3,im3]=nyquist(num,den8);plot(re1,im1,re2,im2,re3,im3);3、z=[-10];p=[0,-1/2];k=1/2;h=zpk(z,p,k);bode(h);gridon;4、h1=tf([2.33],[0.1621]);h2=tf([1],[0.03681]);h3=tf([1],[0.001671]);h=h1*h2*h3;[num,den]=tfdata(h);[mag,phase,w]=bode(num,den);subplot(211);semilogx(w,20*log10(mag));grid;subplot(212);semilogx(w,phase);grid;[gm,pm,wcg,wcp]=margin(mag,phase,w);實驗結(jié)果1、2、3、4、實驗總結(jié)通過本次