MATLAB-考试试题-(1)

MATLAB考试试题(1)产生一个1x10的随机矩阵，大小位于（-55），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）答案：a=10*rand(1,10)-5;b=sort(a,'descend')1.请产生一个100*5的矩阵，矩阵的每一行都是[12345]2.已知变量：A=’ilovematlab’；B=’matlab’,请找出：（A）B在A中的位置。（B）把B放在A后面，形成C=‘ilovematlabmatlab’3.请修改下面的程序，让他们没有for循环语句！A=[123;456;789];[rc]=size(A);fori=1:1:rforj=1:1:cif(A(i,j)8|A(i,j)2)A(i,j)=0;endendend4.请把变量A=[123;456;789]写到文件里(output.xls)，写完后文件看起来是这样的1234567895.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。7.将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。9.a=[12345],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB中最简单的方法计算b，注意最简单哦。1、求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+2w=-3x+8z-5w=9-6x+5y-2z+10w=-8求系数矩阵的秩；求出方程组的解。解：（1）a=[34-7-12];5-742;108-5;-65-210];c=[4;-3;9;-8];b=rank(a)b=4（2）d=a\cd=-1.4841,-0.6816,0.5337,-1.2429即：x=-1.4841;y=-0.6816;z=0.5337;w=-1.24292、设y=cos[0.5+((3sinx)/(1+x^2))]把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线；解：x=linspace(0,2*pi,101);y=cos(0.5+3.*sin(x)./(1+x.*x));plot(x,y)3、设f(x)=x^5-4x^4+3x^2-2x+6（1）取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。（提示：用polyval函数）解：p=[1-43-26];x=linspace(-2,8,100);y=polyval(p,x);plot(x,y);axis([-2,8,-200,2300]);为了便于观察，在y=0处画直线，图如下所示：与y=0直线交点有两个，有两个实根。（2）用roots函数求此多项式的根a=roots(p)a=3.0000,1.6956,-0.3478+1.0289i,-0.3478-1.0289i4、在[-10，10；-10，10]范围内画出函数的三维图形。解：[X,Y]=meshgrid(-10:0.5:10);a=sqrt(X.^2+Y.^2)+eps;Z=sin(a)./a;mesh(X,Y,Z);matlab试卷，求答案一、选择或填空（每空2分，共20分）1、标点符号（）可以使命令行不显示运算结果，（）用来表示该行为注释行。2、下列变量名中（）是合法的。(A)char_1;(B)x*y;(C)x\y;(D)end3、为～，步长为的向量，使用命令（）创建。4、输入矩阵，使用全下标方式用（）取出元素“”，使用单下标方式用（）取出元素“”。5、符号表达式中独立的符号变量为（）。6、M脚本文件和M函数文件的主要区别是（）和（）。7、在循环结构中跳出循环，但继续下次循环的命令为（）。(A)return;(B)break;(C)continue;(D)keyboad二、（本题12分）利用MATLAB数值运算，求解线性方程组(将程序保存为test02.m文件)三、（本题20分）利用MATALAB符号运算完成（将程序保存为test03.m文件）：（1）创建符号函数（2）求该符号函数对的微分；（3）对趋向于求该符号函数的极限；（4）求该符号函数在区间上对的定积分；（5）求符号方程的解。四、（本题20分）编写MATALAB程序，完成下列任务（将程序保存为test04.m文件）：（1）在区间上均匀地取20个点构成向量；（2）分别计算函数与在向量处的函数值；（3）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1andy2”。五、（本题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该函数计算时的和（将总程序保存为test05.m文件）。六、（本题13分）已知求解线性规划模型：的MATLAB命令为x=linprog（c,A,b,Aeq,beq,VLB,VUB）试编写MATLAB程序，求解如下线性规划问题（将程序保存为test06.m文件）：问题补充：卷子的地址看不见符号,能做就做了一些.1、标点符号（;）可以使命令行不显示运算结果，（%）用来表示该行为注释行。2、下列变量名中（A）是合法的。(A)char_1;(B)x*y;(C)x\y;(D)end3、为～，步长为的向量，使用命令（本题题意不清）创建。4、输入矩阵，使用全下标方式用（本题题意不清）取出元素“”，使用单下标方式用（本题题意不清）取出元素“”。5、符号表达式中独立的符号变量为（）。6、M脚本文件和M函数文件的主要区别是（变量生存期和可见性）和（函数返回值）。7、在循环结构中跳出循环，但继续下次循环的命令为（C）。(A)return;(B)break;(C)continue;(D)keyboad二、（本题12分）利用MATLAB数值运算，求解线性方程组(将程序保存为test02.m文件)三、（本题20分）利用MATALAB符号运算完成（将程序保存为test03.m文件）：（1）创建符号函数symsx（2）求该符号函数对的微分；（3）对趋向于求该符号函数的极限；（4）求该符号函数在区间上对的定积分；（5）求符号方程的解。四、（本题20分）编写MATALAB程序，完成下列任务（将程序保存为test04.m文件）：（1）在区间上均匀地取20个点构成向量；（2）分别计算函数与在向量处的函数值；（3）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1andy2”。五、（本题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该函数计算时的和（将总程序保存为test05.m文件）。六、（本题13分）已知求解线性规划模型：的MATLAB命令为x=linprog（c,A,b,Aeq,beq,VLB,VUB）试编写MATLAB程序，求解如下线性规划问题（将程序保存为test06.m文件）：[例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)[例2.4]给定系统125.032)(2323sssssssG，求系统的零极点增益模型和状态空间模型，并求其单位脉冲响应及单位阶跃响应。解：num=[1213];den=[10.521];sys=tf(num,den)%系统的传递函数模型Transferfunction:s^3+2s^2+s+3-----------------------------s^3+0.5s^2+2s+1sys1=tf2zp(num,den)%系统的零极点增益模型sys1=sys2=tf2ss(sys)%系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式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:nifreal(p(i))0Flagz=1;endenddisp('系统的零极点模型为');z,p,k系统的零极点模型为ifFlagz==1disp('系统不稳定');elsedisp('系统是稳定的');end运行结果为:系统是稳定的step(A,B,C,D)％系统的阶跃响应资源与环境工程学院2008级硕士研究生《MatLab及其应用》试题注意，每题的格式均须包含3个部分a.程序（含程序名及完整程序）：b.运行过程：c.运行结果：(1)求解线性规划问题：0744357421321321321321x,xxxxxxxxxx.t.sxxxZmin问各xi分别取何值时，Z有何极小值。（10分）答：fprintf('线性规划问题求解\n');f=[-4;1;7];A=[3,-1,1;1,1,-4;];b=[4,-7]';Aeq=[1,1,-1];beq=[5]';lb=[0,0,];ub=[];x=linprog(f,A,b,Aeq,beq,lb,ub);xz=f'*x;fprintf('MINz=%f\n',z);运行结果：线性规划问题求解Optimizationterminatedsuccessfully.x=2.25006.75004.0000MINz=25.750000(2)编写一个函数，使其能够产生如下的分段函数：xxxxxxf65.0620.251.525.0)(，，，，并调用此函数，绘制曲线范围的，在2)()(2][0xfxfx。（10分）答：functiony=f(x)ifx=2y=0.5*x;elseifx6y=0.5;elsey=1.5-0.25*x;endend运行结果x=2f(x)=1x=0:0.05:2;y=diag(A2(x)'*A2(x+2));plot(x,y);xlabel('\bfx');ylabel('\bfy');(3)将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。（10分）①]22[)cos(，，xx（曲线图）；②4)y2,-4x(-242),(2222；yxyxf（曲面图）。答：subplot(2,2,2)