发布时间 : 星期一 文章Matlab考题题整理带答案 - 图文更新完毕开始阅读
WORD格式
根据专业方向特色和相关科研工作需求，经过与导师商量后，结合一个课题具 体任务，编写一份Matlab应用工作报告。报告由：a课题任务要求，b技术路 线，c程序，d运行结果，e总结、等部分构成，完成的报告经导师给出简单评 语并签字后缴来。
a,课题任务：研究了一种生物质，油菜秸秆对水溶液中金属离子铜的吸附行
为，分别从pH，用量，温度几个方面考察秸秆的吸附性，并对分析的最佳条件进行了探讨。同时从吸附热力学和吸附动力学角度探讨吸附机理。结果表明，
100ml溶液pH=5.30，秸秆用量0.75g时，秸秆对铜的吸附量可达到6mg/g左右。 b,技术路线：通过实验，获得一系列的数据，然后通过Matlab来做各种关 系图。从图中找到g各种关系式。
c,程序：x=[2.202.723.444.135.38] y=[2.393.836.076.396.84] plot(x,y); xlabel('pH'); ylabel('吸附量')
图1
x=[0.50.751.01.251.5]y=[6.056.195.334.694.02] plot(x,y);
xlabel('秸秆用量g'); ylabel('吸附量')
专业资料整理
WORD格式
图2
通过数据图，得到比较理想的实验条件pH和秸秆用量，接下来做动力
学和等温线。
>x=[0.1670.5123458]
y=[0.0332010.0860590.1697790.3220610.4807690.6441220.809061
10.]
plot(x,y); xlabel('时间t'); ylabel('时间/吸附量')
图3
x=[0.23630.154960.136190.129060.133730.13315]y=[0.252180.047070.020140.012670.008810.00706] plot(x,y);
xlabel('1/吸附量'); ylabel('1/平衡浓度')
图4
x=[0.626540.809770.865850.88920.873770.87564] y=[0.598291.32731.695891.897372.055032.15149] plot(x,y);
xlabel('Lg吸附量'); ylabel('Lg平衡浓度')
专业资料整理
WORD格式
图5
d,总结：从图1和图2，分析看可以得到比较理想的对于本次实验的pH和 秸秆用量。后面实验是在前面的基础上得到的。图3是吸附动力学反应速率图， 从图中可以看到线性拟合程度很好，符合二级反应速率方程。图4和图5是吸附 等温线作图，看以看出图4的线性拟合较图5的好，说明符合Langmuir吸附等 温模型。
[例2.1]已知SISO系统的状态空间表达式为(2-3)式，求系统的传递函数。 A=[010;001;-4-3-2];B=[1;3;-6];C=[100];D=0; [num,den]=ss2tf(a,b,c,d,u) [num,den]=ss2tf(A,B,C,D,1)
[例2.2]从系统的传递函数(2-4)式求状态空间表达式。 num=[153]; den=[1234];
[A,B,C,D]=tf2ss(num,den) [例2.3]对上述结果进行验证编程。 %将[例2.2]上述结果赋值给A、B、C、D阵； A=[-2-3-4;100;010]；B=[1;0;0]；C=[153]；D=0； [num,den]=ss2tf(A，B，C，D,1)
32
s2ss3 [例2.4]给定系统
G(s)，求系统的零极点增益模型和状态空间模型，
32
s0.5s2s1
并求其单位脉冲响应及单位阶跃响应。
解：
num=[1213];den=[10.521];
sys=tf(num,den)%系统的传递函数模型 Transferfunction: s^3+2s^2+s+3
----------------------------- s^3+0.5s^2+2s+1
sys1=tf2zp(num,den)%系统的零极点增益模型sys1=
sys2=tf2ss(sys)%系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式
专业资料整理
WORD格式
impulse(sys2)%系统的单位脉冲响应 step(sys2)%系统的单位阶跃响应
[例3.1]对下面系统进行可控性、可观性分析。 解：
a=[-1-22;0-11;10-1];b=[201]';c=[120] Qc=ctrb(a,b)%生成能控性判别矩阵 rank(Qc)%求矩阵Qc的秩 ans=3%满秩，故系统能控
Qo=obsv(a,c)%生成能观测性判别矩阵 rank(Qo)%求矩阵Qo的秩 ans=3%满秩，故系统能观测
[例3.2]已知系统状态空间方程描述如下：
试判定其稳定性，并绘制出时间响应曲线来验证上述判断。 解：
A=[-10-35-50-24;1000;0100;0010]; B=[1;0;0;0];C=[172424];D=[0]; [z,p,k]=ss2zp(A,B,C,D,1); Flagz=0; n=length(A); fori=1:n ifreal(p(i))>0 Flagz=1; end end
disp('系统的零极点模型为');z,p,k 系统的零极点模型为 ifFlagz==1
disp('系统不稳定'); elsedisp('系统是稳定的'); end 运行结果为: 系统是稳定的
step(A,B,C,D)％系统的阶跃响应。
专业资料整理