线性系统的频域分析MATLAB实验.docVIP

武汉工程大学实验报告 专业 班号 组别 指导教师 陈艳菲 姓名 同组者 实验名称 线性系统的频域分析 实验日期 第 次实验 实验目的 实验内容 实验结果及分析 实验心得与体会 一、实验目的 1．熟练掌握用MATLAB语句绘制频域曲线。 2．掌握控制系统频域范围内的分析校正方法。 3．掌握用频率特性法进行串联校正设计的思路和步骤。 二，实验内容。 某单位负反馈控制系统的开环传递函数为，试设计一超前校正装置， 使 校正后系统的静态速度误差系数，相位裕量，增益裕量 绘制伯德图程序，以及计算穿越频率，相位裕量 ans = 相位 Inf 9.0406 频率 Inf 3.1425 e=5; r=50; r0=9; [gm1,pm1,wcg1,wcp1]=margin(num0,den0); phic=(r-r0+e)*pi/180; [gm1,pm1,wcg1,wcp1]=margin(num0,den0); alpha=(1+sin(phic))/(1-sin(phic))[gm1,pm1,wcg1,wcp1]=margin(num0,den0); alpha =6.1261 [gm1,pm1,wcg1,wcp1]=margin(num0,den0);[gm1,pm1,wcg1,wcp1] num0=20; den0=[2,1,0]; w=0.1:1000;margin(num0,den0) grid; 原系统的伯德图： num/den = 1.2347 s + 1 ------------- 0.20154 s + 1 校正之后的系统开环传递函数为: num/den = 6.1734 s + 5 ------------------------------------------- 0.20154 s^4 + 1.6046 s^3 + 3.4031 s^2 + 2 s alpha =6.1261; [il,ii]=min(abs(mag1-1/sqrt(alpha))); wc=w( ii); T=1/(wc*sqrt(alpha)); numc=[alpha*T,1]; denc=[T,1]; [num,den]=series(num0,den0,numc,denc); [gm,pm,wcg,wcp]=margin(num,den); printsys(numc,denc) disp(D￡?y??oóμ??μí3?a?·′?μYoˉêy?a:);printsys(num,den) [mag2,phase2]=bode(numc,denc,w); [mag,phase]=bode(num,den,w); subplot(2,1,1);semilogx(w,20*log10(mag),w,20*log10(mag1),--,w,20*log10(mag2),-.); grid; ylabel(·ù?μ(db)); title(--Go,-Gc,GoGc); subplot(2,1,2); semilogx(w,phase,w,phase1,--,w,phase2,-,w,(w-180-w),:); grid; ylabel(?à??(0)); xlabel(?μ?ê(rad/sec)); title([D￡?y?°￡o·ù?μ?￡á?=,num2str(20*log10(gm1)),db,?à???￡á?=,num2str(pm1),0; D￡?yoó￡o·ù?μ?￡á?=,num2str(20*log10(gm)),db,?à???￡á?=,num2str(pm),0]); 矫正后系统的伯德图 矫正之前系统单位阶跃响应 矫正之后系统的单位阶跃响应： 比较矫正前后系统的响应情况：可以看出超前矫正使系统的调节时间变短，响应更加迅速，但是超调量偏大，对改善系统的动态性能起到了巨大的作用。 某单位负反馈控制系统的开