一、題目?

二、matlab數(shù)字仿真

clear; close all; clc %% 基本參數(shù)定義 gamma=45; %給定的相位裕度 data1=5; %補(bǔ)償修正量 data2=5; %補(bǔ)償修正量 s=tf('s'); K=100;%開環(huán)增益 num=K; den1=[1,0]; den2=[0.1,1]; den3=[0.025,1]; den=conv(den1,conv(den2,den3));%分母部分 G0=tf(num,den);%校正前傳函 %% 超前校正 [magG01,phG01,wG01]=bode(G0);%導(dǎo)出bode圖中的各個(gè)數(shù)據(jù) magG01=20*log10(squeeze(magG01)); phG01=squeeze(phG01);%刪除長(zhǎng)度為 1 的維度 [G0m1,P0m1,W0cg1,W0cp1]=margin(G0);%Gm為幅值裕度，Pm為相角裕度，Wcg為穿越頻率，Wcp為截止頻率 psim=gamma-P0m1+data1; %設(shè)置校正裝置最大相角 a=(1+sind(psim))/(1-sind(psim));%求解超前校正系數(shù)a mag_wm=-10*log10(a);%尋找穿越頻率 fun=@(w)(20*log10(abs(freqresp(G0,w)))-mag_wm)^2;%定義一個(gè)隱名函數(shù) lb=0.001;%定義bode圖的下屆范圍 ub=100;%定義bode圖的上屆范圍 wm=fminbnd(fun,lb,ub);%在[lu,ub]范圍中找到fun函數(shù)返回值w的最小值 w2=wm*sqrt(a);%求解校正裝置的第二轉(zhuǎn)折頻率 w1=w2/a;%求解校正裝置的第一轉(zhuǎn)折頻率 Gc1=tf([1/w1,1],[1/w2,1]);%校正傳函 G01=Gc1*G0;%串聯(lián)超前校正后的傳函%% 滯后校正 lb=0.001;%定義bode圖的下屆范圍 ub=100;%定義bode圖的上屆范圍 [magG0,phG0,wG0]=bode(G01,{lb,ub});%導(dǎo)出bode圖中的各個(gè)數(shù)據(jù) magG0=20*log10(squeeze(magG0)); phG0=squeeze(phG0);%刪除長(zhǎng)度為 1 的維度 [G0m,P0m,W0cg,W0cp]=margin(G01); gamma0_wc=gamma+data2;%設(shè)置校正裝置最大相角 fun=@(w)(180+angle(freqresp(G01,w))/(2*pi)*360-gamma0_wc)^2;%定義隱名函數(shù) wc=fminbnd(fun,lb,ub);%尋找隱名函數(shù)返回值最小的wc mag_wc=20*log10(abs(freqresp(G01,wc)));%求原系統(tǒng)wc處的幅值 b=10^(-mag_wc/20);%求解滯后校正系統(tǒng)系數(shù) w2=0.1*wc;%校正裝置第二轉(zhuǎn)折頻率，校正點(diǎn)往前十倍頻程 w1=b*w2;%矯正系統(tǒng)第一轉(zhuǎn)折頻率（1/b倍頻關(guān)系） Gc2=tf([1/w2,1],[1/w1,1]);%寫出滯后校正傳函Gc=Gc1*Gc2;%超前滯后校正傳遞函數(shù) G=Gc*G0;%校正后系統(tǒng)傳函 damp(G) Gc1 Gc2 %% 導(dǎo)出超調(diào)量，調(diào)節(jié)時(shí)間變量 Gf=feedback(G,1);%單位負(fù)反饋 max0=max(step(Gf))-1;%尋找超調(diào)量 max0=max0*100; S2 = stepinfo(Gf,'RiseTimeThreshold',[0.05 0.95]);%導(dǎo)出單位階躍響應(yīng)參數(shù) SSS=S2.SettlingTime;%導(dǎo)出調(diào)節(jié)時(shí)間%% 作圖 figure(1)%圖1 step(Gf)%調(diào)節(jié)后傳函單位負(fù)反饋后單位階躍響應(yīng) figure(2)%圖2 bode(G0,Gc,G,{lb,ub})%頻率范圍為0.001到100的bode圖 grid on%加網(wǎng)格線 [Gm,Pm,Wcg,Wcp]=margin(G);%Gm為幅值裕度，Pm為相角裕度，Wcg為穿越頻率，Wcp為截止頻率 title(['滯后超前校正: \gamma^*=',num2str(gamma),...'\gamma=',num2str(Pm),'\circ, h=',num2str(-20*log10(1/Gm)),...' dB at w_c=',num2str(Wcp),' rad/s.']);%加圖頭 legend('原系統(tǒng)G0','校正傳函Gc','已校正系統(tǒng)G')%圖例 disp('校正函數(shù)為：'),Gc disp('校正后函數(shù)為：'),G disp('超調(diào)量max0：'),max0 disp('調(diào)節(jié)時(shí)間：'),SSS %命令行輸出內(nèi)容

?三、工程實(shí)現(xiàn)與電路設(shè)計(jì)

寫在實(shí)驗(yàn)報(bào)告冊(cè)上了，需要的私聊

總結(jié)

以上是生活随笔為你收集整理的利用滞后——超前系统解决高阶随动系统（课程设计）的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。