昆明理工大学MATLAB实验指导书(第二次实验)
Matlab第二次实验报告
Matlab第二次实验报告实验目的：1、了解plot函数和subplot函数的基本用法和matlab绘图的基本原理。
2、了解图形的属性设置。
比如画图的颜色，画图采用的线性标识符等。
二：实验基本知识1：1. 单窗口单曲线绘图；2. 单窗口多曲线绘图；3. 单窗口多曲线分图绘图；4. 多窗口绘图；5.可任意设置颜色与线型；6.图形加注功能；7.fplot——绘制函数图函数；8.ezplot——符号函数的简易绘图函数9：subplot—将画图区域分块函数。
实验内容：将高数课本后的18个图画出上机练习程序第1——6图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1);x=-50:1:50;y=a.*x.^3;subplot(3,2,1);plot(x,y,'r');title('y=a*x^3');xlabel('X');ylabel('Y');x=0:1:50;y=sqrt(a.*x.^3);subplot(3,2,2);plot(x,y,x,-y,'r');title('y^2=a*x^3');xlabel('X');ylabel('Y');x=-3:0.1:3;y=exp(-x.^2);subplot(3,2,3);plot(x,y,'r');title('y=e^x^2');xlabel('X');ylabel('Y');x=-3:0.1:3;y=8*a^3./(x.^2+4*a^2);subplot(3,2,4);plot(x,y,'r');title('y=8*a^3./(x.^2+4*a^2)'); xlabel('X');ylabel('Y');x=0:0.01:5;y=sqrt(x.^3./(a*2.-x)); subplot(3,2,5);plot(x,y,x,-y,'r');title('y^2*(2a-x)=x^3'); xlabel('X');ylabel('Y');%t=-1:0.01:5;%x=3*a.*t./(1+t.^3);%y=3*a*t.^2./(1+t.^3);subplot(3,2,6);ezplot('x.^3+y.^3-3*3*x.*y'); %plot(x,y);title('x^3+y^3-3axy=0'); xlabel('X');ylabel('Y');输入a=3：显示结果第7——10图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0:0.1:2*pi;x=a*cos(t).^3;y=a*sin(t).^3;subplot(2,2,1);plot(x,y,'b')title('x^1.5+y^1.5=a^1.5'); xlabel('X');ylabel('Y');t=-2*pi:0.1:2*pi;x=a.*(t-sin(t));y=a.*(1-cos(t));subplot(2,2,2);plot(x,y,'b')title('°ÚÏß');xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*(1-cos(t));subplot(2,2,3);plot(p.*cos(t),p.*sin(t),'b'); title('p=a(1-cos£¨t£©)'); xlabel('X');ylabel('Y');t=0:0.01:2*pi;p=a.*t;subplot(2,2,4);plot(p.*cos(t),p.*sin(t),'b'); title('p=at');xlabel('X');ylabel('Y');输入a=1，第11题图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º') figure(1);t=-2*pi:pi/100:2*pi;p=exp(a.*t);plot(p.*cos(t),p.*sin(t)); title('p=e^at');xlabel('X');ylabel('Y');输入a=0.1，显示结果：clear,clca=input('ÇëÊäÈëaµÄÖµ£º')figure(1)t=0.2*pi:pi/100:100*pi;p=a./t;plot(p.*cos(t),p.*sin(t));title('p=at');xlabel('X');ylabel('Y');输入a=1，显示结果：第13——14题图：clear,clca=input('ÇëÊäÈëaµÄÖµ£º');figure(1);t=0:pi/100:2*pi;p=(a^2.*sin(2*t)).^0.5;subplot(1,2,1);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*sin(2t)');xlabel('X');p=(a^2.*cos(2*t)).^0.5;subplot(1,2,2);plot(p.*cos(t),p.*sin(t),-p.*cos(t),-p.*sin(t),'g'); title('p^2=a^2*cos(2t)');xlabel('X');ylabel('Y');输入a=1，显示结果：第15-16题图：clear;clc;a=input('ÇëÊäÈëa:');figure(1)t=0:pi/200:pi;p=a.*cos(3*t);subplot(1,2,1);plot(p.*cos(t),p.*sin(t),'r--');title('p=a*cos(3t)');xlabel('X');p=a.*sin(3*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'r--'); title('p=a*sin(3t)');xlabel('X');ylabel('Y');输入a=1，显示结果：第17——18题图：clear;clc;a=input('ÇëÊäÈëa:');figure(1);t=-pi:pi/200:pi;subplot(1,2,1);p=a.*sin(2*t);plot(p.*cos(t),p.*sin(t),'k-.'); title('p=a*sin(2t)');xlabel('X');ylabel('Y');p=a.*cos(2*t);subplot(1,2,2);plot(p.*cos(t),p.*sin(t),'k-.');title('p=a*cos(2t)');xlabel('X');ylabel('Y');输入a=1，显示结果：实验心得：这次实验主要是学习matlab软件的二维绘图功能：matlab软件还有丰富的图形修饰功能，如改变线条的形式和颜色；除此之外其还可以以多种形式如在一个图形中显示多个函数图形，一个figure中显示多个坐标系，还有一个程序中就可以用多个figure一起来表示不同的函数……在画图过程中你可能因为角度的区间取值不一样，画出来的图形可能跟书本上给出的图形有很大的差异，但没关系，只要你思路是对的，区间可以自己慢慢改动知道跟书本上给出的图形想接近。
实验二MATLAB程序设计含实验报告
实验二MATLAB 程序设计一、 实验目的1.掌握利用if 语句实现选择结构的方法。
2.掌握利用switch 语句实现多分支选择结构的方法。
3.掌握利用for 语句实现循环结构的方法。
4.掌握利用while 语句实现循环结构的方法。
5.掌握MATLAB 函数的编写及调试方法。
二、 实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。
点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、 实验内容1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。
并输入几组典型值加以检验。
（提示：提示输入使用input 函数）2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。
其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。
要求：（1）用switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
（提示：注意单元矩阵的用法）3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。
重复此过程，最终得到的结果为1。
如：2?13?10?5?16?8?4?2?16?3?10?5?16?8?4?2?1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。
请为关键的Matlab 语句填写上相关注释，说明其含义或功能。
Matlab实验指导书含答案
实验一：M a t l a b操作环境熟悉一、实验目的1．初步了解Matlab操作环境。
2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
➢求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)^3*(2-x);➢求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3) y=1/sqrt(1-x^2);(4) y=(x-1)/(x+1)/(x+2);➢求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x^2));➢代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin(x)^2+cos(x)^2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
➢从y=x^2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2(4) y4=x^2+2(5) y5=x^4(6) y6=x^2/23．两个函数之间的操作➢求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5➢乘积(1) exp(-x)*sin(x)(2) sin(x)*x➢商(1) sin(x)/cos(x);(2) x/(1+x^2);(3) 1/(x-1)/(x-2);➢求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab应该注意函数表达式的文本式描述。
现代控制工程MATLAB实验指导书V
昆明理工大学机电学院机械工程专业研究生实验《现代控制工程》MATLA实验指导书昆明理工大学机电工程学院现代控制工程》MATLAB 实验指导书一、MATLAB 简介1 MATLAB 概述MATLAB 是MA Trix LABoratory 的缩写，早期主要用于现代控制中复杂的矩阵、向量的各种运算。
由于MA TLAB 提供了强大的矩阵处理和绘图功能，很多专家因此在自己擅长的领域用它编写了许多专门的MATLAB 工具包( toolbox )，如控制系统工具包( control systems toolbox )；系统辨识工具包( system identification toolbox )；信号处理工具包( signal processing toolbox )；鲁棒控制工具包( robust control toolbox )；最优化工具包(optimization toolbox )等等。
由于MA TLAB 功能的不断扩展，所以现在的MATLAB 已不仅仅局限与现代控制系统分析和综合应用，它已是一种包罗众多学科的功能强大的“技术计算语言( TheLanguage of Technical Computing )”。
MathWorks 公司于1992 年推出了具有划时代意义的MATLAB 4.0 版本，并推出了交互式模型输入与仿真系统SIMULINK ，它使得控制系统的仿真与CAD 应用更加方便、快捷，用户可以方便地在计算机上建模和仿真实验。
1997 年MathWorks 推出的MA TLAB 5.0 版允许了更多的数据结构，1999 年初推出的MA TLAB 5.3 版在很多方面又进一步改进了MATLAB 语言的功能。
2000 年底推出的MATLAB 6.0 。
最新版本是MA TLAB7.0 。
MATLAB 以矩阵作为基本编程单元，它提供了各种矩阵的运算与操作，并有较强的绘图功能。
MATLAB 集科学计算、图像处理、声音处理于一身，是一个高度的集成系统，有良好的用户界面，并有良好的帮助功能。
matlab第二次试验上机操作内容
实验名称：数组、矩阵和基本数据类型
一、实验目的：
1、熟悉数组和矩阵的基本运算；
2、熟悉数值数组、字符数组的创建；
3、熟悉数值数组、字符数组的赋值和寻访。
二、实验内容：
1、对数组和矩阵做各种基本运算
2、创建数值和字符数组；
3、对字符数组进行赋值和寻访，同时实现数值和字符之间的转换。
2)生成2*5的随机矩阵，利用逻辑索引找出其中〉0.4的element
3)对R=[-0.5,-0.3+0.4*i,-0.3-0.4*i];分别求实部，虚部，幅度和相位并记录结果
4)对x=[6.3,-6.3,6.8,-6.8]分别进行取整（ceil，floor，fix，round）运算，并记录结果
5)对 ， 分别进行数组乘，除和矩阵乘除，并记录结果
三、实验平台：MATLAB集成系统
四、设计流程：
启动软件，输入指令，观察结果
五、程序清单：
此处写程序内容
六、调试和测试结果：
此处写程序的执行结果和实验过程中的调试经过、出现的错误和对应的解决方法
七、教师评语与成绩评定
此处由老师填写
实验的详细内容如下：
1)继续第一次没有完成的内容
（输入command前首先在commadwindow输入diary exp_2以保存实验结果供写实验报告用）
6)利用A=magic(4)生成4阶magic数组，对其做提取对角元素（diag），变形（reshape），左右反转，上下反转，旋转90°和重复（repmat）等操作并记录结果
2 2 3;3 3 3;1 1 1];分别作各种关系运算（包括大于，小于，等于，大于等于，小于等于和非）
8)对A=[1 0;2 3];B=[1 1;20];进行各种逻辑运算
MATLAB实验二
3. 三维图形
• 三、视点
➢view(az,el)-------------设置观察点方向
az为方位角，el为仰角。 • 方位角为视点位置在XY平面上的投影与X轴形成的角度，
正值表示逆时针，负值表示顺时针。 • 仰角为XY平面的上仰或下俯角，正值表示视点在XY平
面上方，负值表示视点在XY平面下方。
figure grid gtext hold subplot text title xlabel ylabel
创建图形窗口 放置坐标网格线 用鼠标放置文本 保持当前图形窗口内容 创建子图 放置文本 放置图形标题 放置X轴坐标标记 放置Y轴坐标标记
3. 三维图形 • 一、 plot3函数
• 将二维函数plot的功能扩展到三维空间，绘制三维图形。 • 函数格式：plot3(x1,y1,z1,c1,x2,y2,z2,c2,…)
theta=[0:0.01:2*pi];
rho=sin(2*theta).*cos(2*theta);
polar(theta,rho); 绘制极坐标图命令
title('polar plot');
例9、程序：
theta=linspace(0, 2*pi);
r=cos(4*theta);
polar(theta, r);
所组成的画面。
4. 动画设计
• 【例14】 播放一个不断变化的眼球程序段。
• m=moviein(20); 建立一个20个列向量组成的矩阵
• for j=1:20
•
plot(fft(eye(j+10))) %绘制出每一幅眼球图并保存到m矩阵中
• m(:,j)=getframe;
MatLAB第二次实验
自动控制理论实 验 报 告姓 名 班 级 成 绩_______________学 号 日 期 合作者________________________________________________________________________________________________一、实验题目：利用MATLAB 求系统的的瞬态性能指标 二、实验目的：1、熟悉并掌握MatlAB 软件在工程计算中的应用。
2、求系统的的瞬态性能指标。
三、实验仪器：电脑、MatlAB 软件； 四、实验内容用impulse 函数求单位脉冲响应，step 求求单位阶跃响应。
程序：%a_1.mt=[0:0.01:0.8] %num_G=[50];tao=0;den_G=[0.05 1+50*tao 50];G1=tf(num_G ,den_G);tao=0.0125;den_G=[0.05 50*tao tao 50];G2=tf(num_G ,den_G); tao=0.025;den_G=[0.05 1+50*tao 50];G3=tf(num_G ,den_G); %[y1,T]=impulse(G1,t);[y1a,T]=step(G1,t); [y2,T]=impulse(G2,t);[y2a,T]=step(G2,t); [y3,T]=impulse(G3,t);[y3a,T]=step(G3,t); %subplot(1,2,1)plot(T,y1,'r-',T,y2,'g-',T,y3,'b-')legend('tao=0','tao=0.0125','tao=0.025') xlabel('t(sec)') ylabel('x(t)'); grid on;subplot(1,2,2)plot(T,y1a,'r-',T,y2a,'g-',T,y3a,'b-')legend('tao=0','tao=0.0125','tao=0.025') grid on;xlabel('t(sec)') ylabel('x(t)');二、用lsim函数求任意输入响应程序：%3_2.mclc;t=[0:0.01:1];u=sin(2*pi*t);%tao=0.025;num_G=[50];den_G=[0.05 1+50*tao 50];G=tf(num_G,den_G);%y=lsim(G,u,t);%plot(t,u,'r-',t,y,'g-',t,u','b',t,y','r',t,u'-y,'b-','linewidth',1)legend('u(t)','xo(t)','u(t)','xo(t)','e(t)=u-y')grid;xlabel('t(sec)')ylabel('x(t)');三、利用MATLAB求系统的的瞬态性能指标程序：close all;%t=0:0.001:1;%yss=1;dta=0.02;%num_G=[50];tao=0;den_G=[0.05 1+50*tao 50];G1=tf(num_G,den_G);tao=0.0125;den_G=[0.05 1+50*tao 50];G2=tf(num_G,den_G);tao=0.025;den_G=[0.05 1+50*tao 50];G3=tf(num_G,den_G);y1=step(G1,t);y2=step(G2,t);y3=step(G3,t);%r=1;while y1(r)<yss;r=r+1;endtr1=(r-1)*0.001;%[ymax,tp]=max(y1);tp1=(tp-1)*0.001;%mp1=(ymax-yss)/yss;%s=1001;while y1(s)>1-dta & y1(s)<1+dta;s=s-1;endts1=(s-1)*0.001;%τ=0.0125r=1;while y2(r)<yss;r=r+1;endtr2=(r-1)*0.001;[ymax,tp]=max(y2);tp2=(tp-1)*0.001;mp2=(ymax-yss)/yss;s=1001;while y2(s)>1-dta &y2(s)<1+dta;s=s-1;endts2=(s-1)*0.001;%τ=0.025r=1;while y3(r)<yss;r=r+1;endtr3=(r-1)*0.001;[ymax,tp]=max(y3);tp3=(tp-1)*0.001;mp3=(ymax-yss)/yss;s=1001;while y3(s)>1-dta & y3(s)<1+dta;s=s-1;endts3=(s-1)*0.001;%[tr1 tp1 mp1 ts1;tr2 tp2 mp2 ts2;tr3 tp3 mp3 ts3]结果：ans =0.0640 0.1050 0.3509 0.3530 -------系统10.0780 0.1160 0.1523 0.2500 -------系统20.1070 0.1410 0.0415 0.1880 -------系统3任课教师：。
MATLAB实验指导书
实验一、Matlab 软件的安装与基本使用一、实验目的和要求：学会MATLAB 软件的安装方法，熟悉MA TLAB 的工作环境及其基本功能的设置，以及基本指令的使用，了解MATLAB 数值数组及其矩阵操作的使用、复数的输入方法，了解基本的绘图函数程序。
二、实验内容和步骤：1、完成MATLAB 软件的安装。
2、安装完成后启动MA TLAB 软件，熟悉MATLAB 的工作环境。
打开MATLAB 7，默认打开的窗口包括： (1)命令窗口(Command Window)； (2)命令历史窗口(Command History)； (3)工作间管理窗口(Workspace)； (4)当前路径窗口(Current Directory)此外，还有捷径键、编译窗口、图形窗口和帮助窗口等其他种类的窗口，结合课本实例，熟悉每个窗口的功能和基本使用。
3、（1）求23)]47(212[÷-⨯+的运算结果，注意命令后有无分号的作用。
命令及结果如下：（2）求20012x x v t a t =++的值，其中t=5s 、010x m =、015/v m s =、29.81/sec a m =-。
命令及结果如下：(3) 设3x =, 4v =，用MATLAB 命令求下列格式的值： 1）22log()x v +2)(2)x v -3) |sin 2|vx e4) eps 5)34x v --命令及结果如下：4、比较以下三种不同的输入矩阵1234246836912a ⎛⎫⎪= ⎪ ⎪⎝⎭的方法，运行并观察实验结果，掌握矩阵的基本输入方法。
>> a=[1 2 3 42 4 6 83 6 9 12]>> a = [1:4; 2:2:8; 3:3:12]>> a1 =[1:4]; a = [a1; 2*a1; 3*a1]5、输入复数矩阵⎥⎦⎤⎢⎣⎡-+++=i ii i24356223B ，并求其实部、虚部、模和幅角（分别用弧度和角度表示），了解MATLAB 向量化运算的特点。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
************************ MATLAB上机指导书************************昆明理工大学机电学院彭用新2015年3月实验三符号计算一、操作部分：在命令窗口执行命令完成以下运算，记录运算结果。
1.findsym：帮助我们获取系统定义的自变量f= sym('sin(a*x+b*y)'); findsym(f)2.numden（获取分子分母）, sym2poly,（获取多项式时系数）poly2sym（根据多项式系数获得符号表达式）[n,d]=numden(sym('x*x+y')+sym('y^2'))p=sym('2*x^3+3*x^2+4'); sym2poly(p)x=[2,3,0,4]; poly2sym(x)3. collect ：合并同类项；expand：展开多项式；horner: 分解成嵌套形式；factor：因式分解；simplify: 对表达式化简syms x y; collect(x^2*y+y*x-x^2-2*x)collect((x+y)*(x^2+y^2+1), y)syms x y; expand((x-2)*(x-4))syms x;horner(x^3-6*x^2+11*x-6)syms x;factor(x^3-6*x^2+11*x-6)syms x;simplify((x^2+5*x+6)/(x+2))4. finverse ：求得符号函数的反函数。
syms x y; finverse(1/tan(x))f= x^2+y; finverse(f,y)finverse(f)pose 求符号函数的复合函数syms x y; f = 1/(1 + x^2); g = sin(y); compose(f,g)6. subs :表达式替换。
syms a b;subs(a+b,a,4)subs(a+b,4)subs(cos(a)+sin(b),{a,b},{sym('alpha'),2}) subs('x^2+2*y',{'x','y'},{3,4})7.极限：limitsyms x t;limit(sin(x)/x)limit((x-2)/(x^2-4),2)limit((1+2*t/x)^(3*x),x,inf)limit(1/x,x,0,'right')8. 微分：diffsyms a b c x;f=sym('a*x^2+b*x+c')diff(f)diff(f,2)diff(f,a)diff(f,a,2)9. 积分：intsyms a b c xf=sym('a*x^2+b*x+c')int(f)int(f,x,0,2)int(f,a)int(int(f,a),x)10级数：symsum,taylorsyms ksymsum(1/k,k,1,inf)symsum(1/(k*(k+1)),k,1,inf)syms xtaylor(sin(x),x,10)subs(ans,x,pi/2)11.solve 解符号方程式f=sym('a*x^2+b*x+c');solve(f)f1=sym('x+y+z-10=0');f2=sym('x-y+z=0');f3=sym('2*x-y-z+4=0');solve(f1,f2,f3);[x,y,z]=solve(f1,f2,f3)12解微分符号方程式：dsolve,Dy代表dy/dt, D2y代表d2y/dt2 。
dsolve('Dy=5')dsolve('Dy=x','x')dsolve('D2y=1+Dy')dsolve('D2y=1+Dy','y(0)=1','Dy(0)=0')[x,y]=dsolve('Dx=x+y, Dy=2*x','x(0)=0','y(0)=1')二、使用部分：利用matlab求解，并记录运行结果1、计算符号表达式()cos()sin()f x x x x=+-在1xπ=-处的值2、设x为符号变量，42()21f x x x=++，32()635g x x x x=+++，试进行如下运算。
（1）()() f x g x+，（2）()() f x g x⨯，（3）对()f x进行因式分解3、指出下面各条语句的输出结果>> f=sym('2*u');>> subs(f,'u',2)>> f2=subs(f,'u','u+2')>> a=3;>> subs(f2,'u',a+2)>> subs(f2,'u','a+2')>> syms x y>> f3=subs(f,'u',x+y)>> subs(f3,[x,y],[1,2])>> subs(f3,[x,y],[x+y,x+y])>> f=sym('2*u');4、化简5、 因式分解：6555234-++-x x x x6、求极限 ：7、已知f(x)=ax2 +bx+c,求f(x)的微分8、已知f(x)=ax2 +bx+c,求f(x)的积分。
9、实验四 矩阵运算一、实验目的：掌握基本的矩阵运算及常用的函数。
二、实验内容：⎥⎦⎤⎢⎣⎡=654321a ⎥⎦⎤⎢⎣⎡-=531142b ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=201c ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=063258741d1、下列运算是否合法，为什么？如合法，结果是多少？ (1) result1 = a' (2) result2 = a * b (3) result3 = a + b (4) result4 = b * d(5) result5 = [b ; c' ] * d (6) result6 = a . * b (7) result7 = a . / b (8) result8 = a . * c (9) result9 = a . \ b (10) result10 = a . ^2 (11) result11 = a ^2 (12) result11 = 2 . ^ a2、用MATLAB 求下面的的方程组。
(1)⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⨯⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----017413231511222315921274321x x x x()f x =lim()xx x a x a→∞+-11111111............23k 22x33x4k(k+1)++++++++++和的和。
(2)⎪⎪⎩⎪⎪⎨⎧=-++=--=-++=++56533332821w z y x w y x w z y x z y x 3、已知⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----=1323151122231592127A (1)求矩阵A 的秩(rank)(2)求矩阵A 的行列式(determinant) (3)求矩阵A 的逆(inverse)(4)求矩阵A 的特征值及特征向量(eigenvalue and eigenvector) 4、关系运算与逻辑运算 已知a=20,b=-2,c=0,d=1(1) r1 = a > b(2) r2 = a > b & c > d (3) r3 = a == b* (-10) (4) r4 = ~b | c三、思考题109910101022222++++==---=∑ n n y ,求y=？（用format long 查看y 的值）预习 MATLAB 图形绘制基础一、实习目的1、进一步熟悉M 文件调试过程2、熟练掌握MATLAB 二维曲线的绘制3、掌握图形的修饰二、实习内容1、基本命令：plot （数据作图、函数作图）, ezplot （简易作图），fplot(简易作图) （1）plot(x ,y,‟s‟) s 用来表示色彩、线型、点型 例3.1 用plot 绘制y=sinx/x 的图形 x=-2*pi:pi/10:2*pi; y=sin(x)./x; %准备数据 plot(x,y) %用plot 作图 注意x=0处的图形 例3.2 x=-2*pi:pi/10:2*pi;x=x+(x==0)*eps;y=sin(x)./x; %（准备数据,处理分母零点))plot(x,y) %（用plot作图）：例3.3 画出()π,0包罗线。
t=(0:pi/100:pi);y1=sin(t);y=-sin(t);y2=sin(t).*sin(9*t);t3=pi*(0:9)/9;y3=sin(t3).*sin(9*t3);plot(t,y1,'r:',t,y,'r:',t,y2,'b',t3,y3,'bo'),axis([0,pi,-1,1]) gtext(…该点数值x‟)运行结果见图3-1。
（2）ezplot用作图：ezplot(…sin(x)/x‟,[-4*pi,4*pi])运行结果如图3-2（3）fplot用作图：fplot(…sin(x)/x‟,[-2*pi,2*pi])图3-1图3-22、基本绘图控制参数：设置线型、线色、数据点形（1）、曲线点型、线色允许设置值（2（3（4）、图形标注：坐标轴名(label)、图形标题（title ）、图例(legend)、图形说明（text,gtext ）,格式如下：title(…s‟) 书写图形标题 xlabel(…s‟)横坐标轴名 ylabel(…s‟) 纵坐标轴名legend(…s1‟,‟s2‟,…) 绘制曲线所用线型、色彩、或数据点型图例 text(x1,y1,‟s‟) 在图面（x1,y1）坐标处书写字符注释 gtext(…s‟) 用鼠标选择书写字符注释的位置 （5）、常用坐标控制命令（hold on 使当前轴及图形保持而不被刷新，准备接受此后将绘制的新曲线。
hold off 使当前轴及图形不再具备不被刷的性质。
hold 是否具备不被刷的性质的双向开关。
subplot(m,n,k) 图形窗中有m×n 幅子图，k 是子图的编号。