计算机仿真大作业姚雷阳10班2【参考】.docVIP

《计算机仿真》 MATLAB课程设计报告学院：自动化学院专业：自动化专业班级：20姓名：学号：2013年12月2日 Matlab程序文件的写入和读出、、Simulink建立子系统并对子系统进行封装，以及对PID控制器的参数进行整定这三个方面。 PID的参数整定要用到Simulink建立系统，并且要设计M文件，采用编程法求取系统阶跃响应的性能指标。 关键词： PID整定法 Simulink 子系统封装 M文件目录 一、基本操作 1 1、 1 2、 2 3、 2 4、 4 二、子系统封装 4 1、 5 2、观察输入输出波形 6 3、数据导出 7 三、PID 控制器参数整定 9 1、采用工程整定法分别设置控制器参数 91.1先建立原控制系统的Simulink仿真模型 91.2对原控制系统进行仿真 91.3在该控制系统中引入PID控制器 101.4求取系统临界稳定增益并设置参数 11 2、PID控制仿真并绘图 132.1 P控制仿真 132.2 PI控制仿真 152.3PID控制仿真 16 3、制作M文件求取系统阶跃响应性能指标 17 四、总结 19 五、参考文献 20一、基本操作 1、1.1 自己定义一组数据，并将其保存到文件data.dat。要求第一列为时间t(t为等差数列，0≤t≤200)；第二列为与t对应的201个幅值数据，作为信号f 1 (t )的幅值；第三列为按s的降幂排列的传递函数分子系数；第四列为按s的降幂排列的分母系数。第三列、第四列的数据个数不能超过5个。 MATLAB文本： a=0:1:200; b= random(Poisson,0:200,1,201); c=[9 8 6 2 1]; d=[7 5 4 3 1]; fid=fopen(d:\data.dat,wt); fprintf(fid,%4d,a); fprintf(fid,\n); fprintf(fid,%4d,b); fprintf(fid,\n); fprintf(fid,%4d,c); fprintf(fid,\n); fprintf(fid,%4d,d);fclose(fid);图1 data.dat 的读入 2、2.1读入data.dat数据，画出f 1 (t )的时域波形。 MATLAB文本： fid=fopen(d:\data.dat,rt); t=fscanf(fid,%d,201); f1=fscanf(fid,%d,201); fclose(fid); plot(t,f1)title(f1-t曲线图);xlabel(t);ylabel(f1); 图2 f1曲线图 3、 3.1 求取f 2 (t )，将结果保存到result.mat文件，画出其时域波形。 MATLAB文本： for n=1:10f2(n)=2*f1(n); end for n=11:100 if abs(f1(n))15f2(n)=f1(n)*f1(n); else f2(n)=f1(n-3); end end for n=101:201f2(n)=f1(n-3); end plot(t,f2)title(f2-t曲线图); xlabel(t);ylabel(f2); fid=fopen(d:\result.mat,wt);fprintf(fid,%4d,f2); fclose(fid); 图3 f2曲线图 4、4.1 按 data.dat中的第三列、第四列，求取其对应的传递函数，绘制其bode图。 fid=fopen(d:\data.dat,rt);f=fscanf(fid,%d,402);num=fscanf(fid,%d,5);den=fscanf(fid,%d,5); fclose(fid); num1=ctranspose(num);den1=ctranspose(den); Gs=tf(num1,den1) Gs = 9 s^4 + 8 s^3 + 6 s^2 + 2 s + 1 ------------------------------- 7 s^4 + 5 s^3 + 4 s^2 + 3 s + 1Continuous-time transfer function. bode(Gs) 图4 系统波特图 二、子系统封装 1、.1用 SIMULINK 建立系统：y = ax 3+ bx + c，x为输入，y为输出,a、b、c 为常数。对该系统进行封装，要求通过对话框能修改a、b、c的值。图5　 未封装系统图6 封装后系统图7 参数更改界面 2、观察输入输出波形 2.1若输入x为幅度为5、频率为0.25Hz的锯齿波，采用示波器显示输出y及输入x的波形。