- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
例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(the
文档评论(0)