資源描述：
《第12章例題源程序》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫(kù)。
1、例12.1三相不平衡交流電路分析。計(jì)算圖12.1所示三相不平衡交流電路各支路電流()并繪制相量圖。其中ra=rb=rc=5Ω,rab=6Ω,rbc=10Ω,rca=15Ω,E=220V。用MATLAB先求，進(jìn)而求各支路電流:a=[16,-5,-5;-5,20,-5;-5,-5,25];%系數(shù)矩陣b=[220;-110-110*sqrt(3)*i;-110+110*sqrt(3)*i];I=inv(a)*b;%解方程Ia=I(1)-I(3)%求各支路電流Ib=I(2)-I(1)Ic=I(3)-I(2)h=compass([Ia,Ib,Ic]);%繪制相量圖set(h,'LineWid
2、th',2);例12.2調(diào)諧振蕩電路分析。分析圖12.3所示的調(diào)諧振蕩電路（iL=f(v)=，β>0，γ>0），要求繪制振蕩波形和相軌跡?；谝陨蠣顟B(tài)方程建立函數(shù)文件vdpol.m：functionydot=vdpol(t,y)ydot(1)=0.1*(1-y(2)^2)*y(1)-y(2);%μ的值可以任意變化，此處取0.1ydot(2)=y(1);ydot=ydot';求解微分方程，并繪制振蕩波形(t,y)和相軌跡(y,dy/dt)：t0=0;tf=60;%確定積分區(qū)間y0=[0;0.25];%確定初始條件[t,y]=ode45('vdpol',[t0,tf],y0);%求解微
3、分方程subplot(1,2,1);plot(t,y(:,2));%繪制振蕩波形subplot(1,2,2);plot(y(:,2),y(:,1));%繪制相軌跡例12.3一典型線性反饋控制系統(tǒng)結(jié)構(gòu)如圖12.6所示。解法2：利用MATLAB控制系統(tǒng)工具箱中已經(jīng)定義的一些LTI仿真函數(shù)，編寫如下程序：G=tf(4,[1,2,3,4]);Gc=tf([1,-3],[1,3]);H=tf(1,[0.01,1]);G_o=Gc*G;%構(gòu)造開環(huán)系統(tǒng)的傳遞函數(shù)。G_c=feedback(G_o,H);%構(gòu)造閉環(huán)系統(tǒng)的傳遞函數(shù)。step(G_o);%求開環(huán)系統(tǒng)的階躍響應(yīng)并繪制相應(yīng)的曲線。figu
4、re;step(G_c);%求閉環(huán)系統(tǒng)的階躍響應(yīng)并繪制相應(yīng)的曲線。例12.45自行車輪飾物的運(yùn)動(dòng)軌跡。為了使平淡的自行車增添一份美感，同時(shí)，也為了增加自行車的安全系數(shù)，一些騎車的人及自行車廠家在自行車的輻條上安裝一款亮麗奪目的飾物，當(dāng)有這種飾物的自行車在馬路上駛過時(shí)，這飾物就象游龍一樣，對(duì)街邊的行人閃過一道波浪形的軌跡。這一波一閃的游龍，其軌跡是什么曲線？試畫出它的圖形。當(dāng)自行車在一個(gè)拋物線型的拱橋上通過時(shí)，或是在一拱一拱的正弦曲線上通過時(shí)，這軌跡是什么曲線?試畫出其圖形。程序如下：clear;%第(1)類情況的實(shí)現(xiàn)x0=0:0.01:2;R=0.1;r=0.075;x1=x0-r
5、*sin(x0/R);%計(jì)算f(x)=0時(shí)p點(diǎn)運(yùn)動(dòng)軌跡y1=R-r*cos(x0/R);subplot(3,1,1);plot(x1,y1,x0,0);%繪制運(yùn)動(dòng)軌跡曲線和f(x)曲線xlabel('x1');ylabel('y1');gridon%第(2)類情況的實(shí)現(xiàn)x0=-1:0.01:1;R=0.1;r=0.1;y0=0.2-0.2*x0.^2;%計(jì)算路面曲線fai=atan(-0.4*x0);%求φint=inline('sqrt(1+(-0.4*x).^2)');%定義θ的積分函數(shù)fork=1:length(x0)theta1(k)=quad(int,0,x0(k))/R
6、;%調(diào)用quad函數(shù)求θendx2=x0+R*0.4*x0./sqrt(1+(-0.4*x0).^2)-r*sin(theta1-fai);%運(yùn)動(dòng)軌跡方程y2=y0+R./sqrt(1+(-0.4*x0).^2)-r*cos(theta1-fai);subplot(3,1,2);plot(x2,y2,x0,y0)%繪制運(yùn)動(dòng)軌跡曲線和f(x)曲線xlabel('x2');ylabel('y2');gridon%第(3)類情況的實(shí)現(xiàn)x0=0:0.01:10;R=0.1;r=0.075;y0=0.3*sin(x0);%計(jì)算路面曲線fai=atan(0.3*cos(x0));%求φint=
7、inline('sqrt(1+(0.3*cos(x)).^2)');%定義θ的積分函數(shù)fork=1:length(x0)theta2(k)=quad(int,0,x0(k))/R;%調(diào)用quad函數(shù)求θend%p點(diǎn)運(yùn)動(dòng)軌跡方程x3=x0-0.3*R*cos(x0)./sqrt(1+(0.3*cos(x0)).^2)-r*sin(theta2-fai);y3=0.3*sin(x0)+R./sqrt(1+(0.3*cos(x0)).^2)-r*cos(theta2-fai
