例12.1三相不平衡交流电路分析。计算图12.1所示三相不平衡.docVIP

PAGE 6 PAGE 6 例12.1 三相不平衡交流电路分析。计算图12.1所示三相不平衡交流电路各支路电流()并绘制相量图。其中ra= rb= rc=5Ω, rab=6Ω, rbc=10Ω, rca=15Ω,E=220V 。 用MATLAB先求，进而求各支路电流: a=[16,-5,-5;-5,20,-5;-5,-5,25]; %系数矩阵 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,LineWidth,2); 例12.2 调谐振荡电路分析。分析图12.3所示的调谐振荡电路（iL=f(v)=，β0，γ0），要求绘制振荡波形和相轨迹。 基于以上状态方程建立函数文件vdpol.m： function ydot=vdpol(t,y) ydot(1)=0.1*(1-y(2)^2)*y(1)-y(2); %μ的值可以任意变化，此处取0.1 ydot(2)=y(1); ydot=ydot; 求解微分方程，并绘制振荡波形(t,y)和相轨迹(y,dy/dt)： t0=0;tf=60; %确定积分区间 y0=[0;0.25]; %确定初始条件 [t,y]=ode45(vdpol,[t0,tf],y0); %求解微分方程 subplot(1,2,1);plot(t,y(:,2)); %绘制振荡波形 subplot(1,2,2);plot(y(:,2),y(:,1)); %绘制相轨迹 例12.3 一典型线性反馈控制系统结构如图12.6所示。 解法2：利用MATLAB控制系统工具箱中已经定义的一些LTI仿真函数，编写如下程序： G=tf(4,[1,2,3,4]); Gc=tf([1,-3],[1,3]); H=tf(1,[0.01,1]); G_o = Gc*G ; %构造开环系统的传递函数。 G_c = feedback(G_o, H); %构造闭环系统的传递函数。 step(G_o); %求开环系统的阶跃响应并绘制相应的曲线。 figure; step(G_c); %求闭环系统的阶跃响应并绘制相应的曲线。 例12.4 自行车轮饰物的运动轨迹。为了使平淡的自行车增添一份美感，同时，也为了增加自行车的安全系数，一些骑车的人及自行车厂家在自行车的辐条上安装一款亮丽夺目的饰物，当有这种饰物的自行车在马路上驶过时，这饰物就象游龙一样，对街边的行人闪过一道波浪形的轨迹。这一波一闪的游龙，其轨迹是什么曲线？试画出它的图形。当自行车在一个抛物线型的拱桥上通过时，或是在一拱一拱的正弦曲线上通过时，这轨迹是什么曲线?试画出其图形。 程序如下： clear; %第(1)类情况的实现 x0=0:0.01:2;R=0.1;r=0.075; x1=x0-r*sin(x0/R); %计算f(x)=0时p点运动轨迹 y1=R-r*cos(x0/R); subplot(3,1,1); plot(x1,y1,x0,0); %绘制运动轨迹曲线和f(x)曲线 xlabel(x1);ylabel(y1);grid on %第(2)类情况的实现 x0=-1:0.01:1;R=0.1;r=0.1; y0=0.2-0.2*x0.^2; %计算路面曲线 fai=atan(-0.4*x0); %求φ int=inline(sqrt(1+(-0.4*x).^2));%定义θ的积分函数 for k=1:length(x0) theta1(k)=quad(int,0,x0(k))/R; %调用quad函数求θ end x2=x0+R*0.4*x0./sqrt(1+(-0.4*x0).^2)-r*sin(theta1-fai);%运动轨迹方程 y2=y0+R./sqrt(1+(-0.4*x0).^2)-r*cos(theta1-fai); subplot(3,1,2); plot(x2,y2,x0,y0)