昆明理工大学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 向量化运算的特点。
matlab实验指导书
MATLAB实验指导书指导老师许承东实验一MATLAB基本操作实验目的1、熟悉MATLAB的工作环境；2、掌握MATLAB常用的操作运算符和一些基本操作；3、学会编写M文件。
实验内容e sin3t,其中t的取值范围为[0,4π]。
1、绘制衰减图线y=5.2/t（1）启动MATLAB，如图1.1所示。
图1.1 MATLAB的工作环境（2）直接在命令窗口输入以下代码：（3）程序执行后显示的衰减振荡曲线如图1.2所示。
图1.2 衰减振荡曲线（4）生成M文件。
在历史命令窗口中选中上面所写代码，单击鼠标右键，在弹出菜单中选择Create M-File菜单项，即可创建为M文件，将文件命名为quxian.m保存。
2、向量化和循环结构的比较（1）从工具栏中单击New M-file图标，或从菜单中选择File/New/M-File创建新的M文件，如图1.3所示。
图1.3 创建新的M文件（2）在M文件编辑器中输入以下代码：（3）保存运行。
将文件名改为vectorize_contrast.m保存。
单击Run 命令或直接按F5执行。
（4）执行结果实验二MATLAB基本图形绘制实验目的1、掌握MATLAB二维图形的绘制；2、掌握MATLAB三维图形的绘制。
实验内容1、二维图形的绘制（1）从工具栏中单击New M-file图标，或从菜单中选择File/New/M-File创建新的M文件。
（2）在M文件编辑器中输入以下代码：（3）保存运行。
将文件名改为example_plot.m保存。
单击Run命令或直接按F5执行。
（4）二维图形绘制结果如图2.1所示。
图2.1 二维图形绘制结果2、三维曲面绘制（1）从工具栏中单击New M-file图标，或从菜单中选择File/New/M-File创建新的M文件。
（2）在M文件编辑器中输入以下代码：（3）保存运行。
将文件名改为matlab_script.m保存。
单击Run命令或直接按F5执行。
实验二。典型二阶系统时-频域分析
实 验 报 告
学号： 姓名： 成绩：
一、 实验名称：典型环节的时域分析和频域分析
二、 实验目的：
(1) 了解、掌握matlab 模拟典型环节的基本方法，包括：比例环节、积分环节、一阶微分环节、惯性环节和振荡环节等。
(2) 熟悉各种典型环节的阶跃响应曲线和频域响应曲线 (3) 了解参数变化对动态特性的影响
三、 实验要求：
(1) 一人一机，独立完成实验内容 。
(2) 根据实验结果完成实验报告，并用A4纸打印后上交。
四、 时间： 五、 地点：
实验报告：
一、比例环节的时域分析和频域分析 比例环节的传递函数：()G s k
二、积分环节的时域分析和频域分析
积分环节的传递函数：1()G s s =
(1) 当k=1:3:10时，绘制系统()k
G s =的阶跃响应曲线，分析曲线特点。
三、一阶微分环节的时域分析和频域分析
一阶微分环节的传递函数：G（s）=TS+1-
四、惯性环节的时域分析和频域分析
惯性环节的传递函数：G(s)=
1
s 1
T ( T>0 )
五、典型二阶系统的时域分析和频域分析
典型二阶系统的传递函数：G(s)=
2
2
2s ()
2n
n
n
Y R s s w s w w
ξ=++（）
关键参数：阻尼比ζ ，和自然频率ωn。
《MATLAB程序设计》实验指导书
三、
实验仪器和设备
1、 计算机一台。 2、 MATLAB7.0 以上集成环境。
四、
预习要求
2
《MATLAB 程序设计》实验指导书
1、 复习 MATLAB 的启动与Байду номын сангаас出，熟悉 MATLAB 运行环境。 2、 复习 MATLAB 中矩阵的生成以及矩阵运算的基本原理。
五、
实验内容及步骤
实验内容：
1、 求下列表达式的值 1） z1
5
《MATLAB 程序设计》实验指导书
2） 函数调用 函数文件编制好后，就可调用函数进行计算了。函数调用的一般格式为 [输出实参表]=函数名（输入实参表） 注意：函数调用时各实参出现的顺序、个数，应与函数定义时形参的顺序、个 数一致，否则会出错。函数调用时，实参先传递给形参，然后再执行函数功能。 5、 选择结构 1） if 语句 a） 单分支 if 语句 if 条件 语句块 end b） 双分支 if 语句 if 条件 语句块 1 else 语句块 2 end c） 多分支 if 语句 if 条件 1 语句块 1 elseif 条件 2 语句块 2 …… elseif 条件 n 语句块 n else 语句块 n+1 end 2） switch 语句 switch 表达式 case 结果表 1 语句块 1 case 结果表 2 语句块 2 ……
二、
实验原理
1、 M 文件 用 MATLAB 语言编写的程序，称为 M 文件，它们的扩展名均为.m。M 文件根据 调用方式的不同分为两类，命令文件（Script file）和函数文件（Function file） 。 2、 建立新的 M 文件 启动 MATLAB 文本编辑器有 3 种方法： 1） 单击工具栏上的“New M-File”命令按钮。 2） 从 MATLAB 主窗口的“File”菜单中选择“New”菜单项，再选择“M-file”命令。 3） 在 MATLAB 命令窗口输入命令“edit”。 3、 打开已有的 M 文件 1） 在当前目录窗口选中要打开的 M 文件，双击鼠标左键。 2） 单击 MATLAB 主窗口工具栏上的“Open File”命令按钮，再从弹出的对话框中 选择所需打开的 M 文件。 3） 从 MATLAB 主窗口的“File”菜单中选择“Open”命令，在“Open”对话框中选中 所需打开的文件。 4） 在 MATLAB 命令窗口输入命令“edit 文件名” 。 4、 函数文件 1） 函数文件的基本结构 函数文件由 function 语句引导，其基本结构为： function 输出形参表=函数名（输入形参表） 注释说明部分 函数体语句 注意：函数名的命名规则与变量名相同。当输出形参多于一个时，应用方括号 括起来。
MATLAB实验指导书
MATLAB实验指导书电子与计算机科学技术学院实验一MATLAB编程环境一、实验目的1．熟悉MATLAB编程环境二、实验环境1．计算机2．MATLAB7.0集成环境三、实验说明1．首先应熟悉MATLAB7.0运行环境，正确操作2．实验学时：2学时四、实验内容和步骤1．实验内容（1）命令窗口的使用。
（2）工作空间窗口的使用。
（3）工作目录、搜索路径的设置。
（4）命令历史记录窗口的使用。
（5）帮助系统的使用。
（6）了解各菜单的功能。
2．实验步骤（1）分析实验内容，写出程序大致框架或完整的程序代码。
（2）进入MATLAB7.0集成环境。
（3）编辑程序并进行保存。
（4）运行程序，若有错误，修改错误后再次运行，如此反复进行到不显示出错为止。
（5）检查程序输出结果。
五、实验报告要求1．整理实验结果，填入相应表格中2．小结实验心得体会实验二矩阵基本运算（一）一、实验目的1．通过实验，进一步熟悉MATLAB编程环境2．通过实验掌握建立矩阵的几种方法3．通过实验理解常用的矩阵运算二、实验环境1．计算机2．MATLAB7.0集成环境三、实验说明1．首先应熟悉MATLAB7.0运行环境，正确操作2．自主编写程序，必要时参考相关资料3．实验学时：2学时四、实验内容和步骤1．实验内容（1）设A和B是两个同维同大小的矩阵，问：1）A*B和A.*B的值是否相等？2）A./B和B.\A的值是否相等？3）A/B和B\A的值是否相等？4）A/B和B\A所代表的数学含义是什么？（2）写出完成下列操作的命令。
1）将矩阵A第2—5行中第1，3，5列元素赋给矩阵B。
2）删除矩阵A的第7号元素。
3）将矩阵A的每个元素值加30。
4）求矩阵A的大小和维数。
5）将向量t的0元素用机器零来代替。
6）将含有12个元素的向量x转换成3*4矩阵。
7）求一个字符串的ASCII码。
8）求一个ASCII码所对应的字符。
（3）下列命令执行后，L1、L2、L3、L4的值分别是多少/A=1:9;B=10-A;L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);2．实验步骤（1）分析实验内容，写出程序大致框架或完整的程序代码。
MATLAB实验二
MATLAB实验二李彤自动化04班学号：201041803042一、实验目的：1. Learn to design branch and loop statements program2. Be familiar with relational and logical operators3. Practice 2D plotting二、实验内容：1. Assume that a,b,c, and d are defined, and evaluate the following expression.a=20; b=-2; c=0; d=1;(1)a>b; 1(2)b>d; 0(3)a>b&c>d; 0(4)a==b; 0(5)a&b>c; 0(6)~~b; 1a=2; b=[1 -2;-0 10]; c=[0 1;2 0]; d=[-2 1 2;0 1 0];(1)~(a>b)0 00 1(2)a>c&b>c;1 00 1⑶ c<=d;??? Error using ==> <= Matrix dimensions must agree. a=2; b=3; c=10; d=0;(9)a*b^2>a*c ; 0(10)d|b>a; 1(11)(d|b)>a; 0(12)isinf(a/b) ; 0(13)isinf(a/c) ; 1(14)a>b&ischar(d) ; 1(15)isempty(c); 02. Write a Matlab program to solve the function1()ln1y xx=-, where x is a number <1. Use an if structure toverify that the value passed to the program is legal. If the value of x is legal, caculate y(x). If not ,write a suitable error message and quit.Solution:x=input('Enter the coefficient x='); if x<1y=log(1/(1-x));fprintf('y=%f\n',y)elsedisp('The value of x is illegal');endEnter the coefficient x=0.1 y=0.105361Enter the coefficient x=1 The value of x is illegal!3. Write out m. file and plot the figures with grids1Assume that the complex function f(t) is defined by the equationf(t)=(0.5-0.25i)t-1.0Plot the amplitude and phase of function for 0 4.t ≤≤ Solution: t=0:0.1:4;x=sqrt((0.5.*t-1).^2+(0.25.*t).^2); y=atan((0.25.*t)./(1-0.5.*t)); plot(t,x); hold on;plot(t,y);4. Write the Matlab statements required to calculate y(t) from the equation22350()350t t y t t t -+≥⎧=⎨+<⎩for value of t between –9 and 9 in steps of 0.5. Use loops and branches to perform this calculation. Solution:for t=-9:0.5:9; if t>=0y=-3*t^2+5; elsey=3*t^2+5; endfprintf('y=%f\n',y); endy=248.000000 y=221.750000 y=197.000000 y=173.750000 y=152.000000 y=131.750000 y=113.000000 y=95.750000 y=80.000000 y=65.750000y=53.000000 y=41.750000 y=32.000000 y=23.750000 y=17.000000 y=11.750000 y=8.000000 y=5.750000 y=5.000000 y=4.250000y=2.000000 y=-1.750000 y=-7.000000 y=-13.750000 y=-22.000000 y=-31.750000 y=-43.000000 y=-55.750000 y=-70.000000 y=-85.750000y=-103.000000 y=-121.750000 y=-142.000000 y=-163.750000 y=-187.000000 y=-211.750000 y=-238.0000005. Write an m.file to evalue the equation 2()32y x x x =-+for all values of x between 0.1 and 3, insteps of 0.1. Do this twice, once with a for loop and once with vectors. Plot the resulting function using a 4.0 thick dashed blue line.Solution:for x=0.1:0.1:3; y=x.^2-3*x+2; plot(x,y,’bo ’); hold on; endx=0.1:0.1:3; y=x.^2-3*x+2;plot(x,y,'b--','LineWidth',4.0);。
MATLAB实验指导书
实验一：MATLAB的基本操作一、实验目的：1.熟悉MATLAB操作环境2．熟悉MATLAB的帮助3．熟悉MATLAB矩阵的操作4．熟悉MATLAB基本运算二、实验内容：1.自行设计试题，说明三角函数是按弧度还是按角度计算。
2.找出取整的函数，并自行设计试题，说明它们的区别。
3.自行设计试题，说明如何给矩阵按列符值。
4.自行设计试题，建一个5X5矩阵，然后将第一行乘1，第二行乘2，第三行乘3，第四行乘4，第五行乘5。
5.设计编写程序（命令），解方程3x5-7x4+5x2+2x-18=0。
6.自行设计试题，说明eval()函数的功能。
三、实验原理：ceil(x)= -4 -2 0 2 5 7fix(x) = -4 -2 0 1 4 6floor(x) = -5 -3 -1 1 4 6round(x) = -5 -2 0 1 5 72、系统的在线帮助help 命令：（1）.当不知系统有何帮助内容时，可直接输入help以寻求帮助:>> help（回车）（2）.当想了解某一主题的内容时，如输入：>> help syntax（了解Matlab的语法规定）（3）.当想了解某一具体的函数或命令的帮助信息时，如输入：>> help sqrt （了解函数sqrt的相关信息）lookfor命令现需要完成某一具体操作，不知有何命令或函数可以完成，如输入：>> lookfor line （查找与直线、线性问题有关的函数）3、常量与变量系统的变量命名规则：变量名区分字母大小写；变量名必须以字母打头，其后可以是任意字母，数字，或下划线的组合。
此外，系统内部预先定义了几个有特殊意义和用途的变量，见下表：数值型向量（矩阵）的输入(1)．任何矩阵（向量），可以直接按行方式．．．输入每个元素：同一行中的元素用逗号（，）或者用空格符来分隔；行与行之间用分号（；）分隔。
所有元素处于一方括号（[ ]）内；例1-1：>> Time = [11 12 1 2 3 4 5 6 7 8 9 10]>> X_Data = [2.32 3.43；4.37 5.98](2)．系统中提供了多个命令用于输入特殊的矩阵：上面函数的具体用法，可以用帮助命令help得到。
MATLAB实验二(修改)
实验二 信号的表示及其基本运算一、实验目的1、掌握连续信号及其MATLAB 实现方法；2、掌握离散信号及其MA TLAB 实现方法3、掌握离散信号的基本运算方法，以及MA TLAB 实现4 熟悉应用MATLAB 实现求解系统响应的方法4、了解离散傅里叶变换的MA TLAB 实现5、了解IIR 数字滤波器设计6、了解FIR 数字滤波器设计1二、实验设备计算机，Matlab 软件三、实验内容（一）、 连续信号及其MATLAB 实现1、 单位冲激信号()0,0()1,0t t t dt εεδδε-⎧=≠⎪⎨=∀>⎪⎩⎰ 例1.1：单位冲击信号的MATLAB 实现程序如下：t1=-4;t2=4;t0=0;dt=0.01;t=t1:dt:t2;n=length(t);x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt;stairs(t,x);axis([t1,t2,0,1.2/dt]);2、 任意函数()()()f t f t d τδττ+∞-∞=-⎰例1.2：用MA TLAB 画出如下表达式的脉冲序列()0.4(2)0.8(1) 1.2() 1.5(1) 1.0(2)0.7(3)f n n n n n n n δδδδδδ=-+-+++++++ 3 单位阶跃函数1,0()0,0t u t t ⎧≥⎪=⎨<⎪⎩例1.3：用MA TLAB 实现单位阶跃函数clear all;t=-0.5:0.001:1;t0=0;u=stepfun(t,t0);plot(t,u)axis([-0.5 1 -0.2 1.2])4 斜坡函数0()()g t B t t =-例1.4：用MA TLAB 实现g(t)=3(t-1)clear all;t=0:0.01:3;B=3;t0=1;u=stepfun(t,t0);n=length(t);for i=1:nu(i)=B*u(i)*(t(i)-t0);endplot(t,u)axis([-0.2 3.1 -0.2 6.2])5 抽样信号 抽样信号Sa(t)=sin(t)/t 在MATLAB 中用 sinc 函数表示。
昆明理工大学计算机仿真实验
《计算机仿真》上机实验报告姓名：学号：专业：测控技术与仪器班级： 121班实验一常微分方程的求解及系统数学模型的转换一．实验目的通过实验熟悉计算机仿真中常用到的Matlab指令的使用方法，掌握常微分方程求解指令和模型表示及转换指令，为进一步从事有关仿真设计和研究工作打下基础。
二. 实验设备个人计算机，Matlab软件。
三. 实验准备预习本实验有关内容（如教材第2、3、5章中的相应指令说明和例题），编写本次仿真练习题的相应程序。
四. 实验内容1. Matlab中常微分方程求解指令的使用题目一：请用MATLAB的ODE45算法分别求解下列二个方程。
要求:1.编写出Matlab 仿真程序；2.画出方程解的图形并对图形进行简要分析；3.分析下列二个方程的关系。
1．2．1、仿真程序function f1=f1(t,x)f1=-x^2;[t,x]=ode45('f1',[0,20],[1]); plot(t,x);gridfunction f2=f2(t,x)f2=x^2;[t,x]=ode45('f2',[0,20],[-1]); figure (2);plot(t,x);grid2.方程解的图形图形进行简要分析3.二个方程的关系两函数关于X 轴对称题目二：下面方程组用在人口动力学中，可以表达为单一化的捕食者-被捕食者模式（例如，狐狸和兔子）。
其中1x 表示被捕食者， 2x 表示捕食者。
如果被捕食者有无限的食物，并且不会出现捕食者。
于是有1'1x x ，则这个式子是以指数形式增长的。
大量的被捕食者将会使捕食者的数量增长；同样，越来越少的捕食者会使被捕食者的数量增长。
而且，人口数量也会增长。
请分别调用ODE45、ODE23算法求解下面方程组。
要求编写出Matlab 仿真程序、画出方程组解的图形并对图形进行分析和比较。
MATLAB 仿真程序：function f3=f3(t,x)f3=[x(1)-0.1*x(1)*x(2)+0.01*t ;-x(2)+0.02*x(1)*x(2)+0.04*t][t,x]=ode45('f3',[0,20],[30;20]);plot(t,x);grid方程组图形：2. Matlab 中模型表示及模型转换指令的使用题目三：若给定系统的的传递函数为1132106126)(23423+++++++=s s s s s s s s G 请用MATLAB 编程求解其系统的极零点模型。
MATLAB实验实验二指导书
实验二MATLAB的程序设计实验目的:（1）掌握函数M文件与命令M文件的建立与编辑；（2）掌握MA TLAB编程语言的基本控制结构与控制转移语句；（3）掌握MA TLAB程序的调试。
实验每组人数及学时：组人数1人，学时数：9学时实验要求：要求每位实验人员于课堂时间完成相关练习实验环境：装有Windows操作系统和MA TLAB软件的计算机。
实验类型：设计型实验理论基础：第三章1、程序流程控制语句for循环结构语法：for i＝初值：增量：终值语句1……语句nend说明：1．i＝初值：终值，则增量为1。
2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。
while 循环结构语法：while 逻辑表达式循环体语句end说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。
若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。
若表达式的值为假，则程序执行end之后的语句。
2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、以便程序能正常执行。
（执行循环体的次数不确定；每一次执行循环体后，一定会改变while后面所跟关系式的值。
）3、while循环也可以嵌套、其结构如下：while逻辑表达式1循环体语句1while逻辑表达式2循环体语句2end循环体语句3endif-else-end分支结构if 表达式1语句1else if 表达式2(可选)语句2else(可选)语句3endend说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。
2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。
3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。
switch-case结构语法：switch表达式case常量表达式1语句组1case常量表达式2语句组2……otherwise语句组nend说明：1．switch 后面的表达式可以是任何类型，如数字、字符串等。
昆明理工大学matlab通信仿真实验二
实验二题目：构建SIMULINK子系统和模块的封装实验目的：学习构建SIMULINK子系统，模块的封装。
传递函数模块的应用。
实验要求：学会构建SIMULINK的子系统：学习自建模块的封装，帮助文档编写，模型空间与Matlab 命令空间的交互。
实验内容：（1）学习信号与系统中的传递函数的应用：拉氏变换和Z变换。
a. 通过仿真求出H(s)=5/(2s+1)的阶跃响应和冲激响应。
构建模型得出波形。
对比采用默认仿真参数以及采用仿真步长为1/1000秒两种情况的波形。
你能够解释吗？程序zgq2_1.mdl仿真步长为default时的仿真运行如下：仿真步长为1/1000秒时的图如下：比较以上两幅仿真波形图，默认步长情况下的图的精度没有1/1000秒时的精细（或者是光滑），这就说明在仿真步长不一样的时候所要求展现的图形的精确度就不同。
b. 将200Hz的锯齿波通过1000次每秒的采样保持器之后得到时间离散信号，观察采样保持器前后的波形变化，解释为什么。
然后将时间离散信号通过传递函数H(z)=4/(z+0.3)的系统，求出响应波形。
仿真步长为1/1000秒。
需要设置所有模块的参数，特别是采样率一定要设置为1/1000秒。
参数设置如下：信号发生器的频率设置为200Hz，Zero-Order Hold的sample time 改为0.001s ，Discrete Transfer Fcn的sample time 改为0.001s,然后将Scope的Time range 改为0.1s.程序zgq2_2.（2）学习构建SIMULINK子系统：构建一个子系统，使得它具有将输入信号m(t)（如一个100Hz的正弦波）和一个常数C相加后再和一个1000Hz的幅度为A的正弦波相乘的功能。
即输出信号：事实上，这就是是个幅度调制子系统。
y(t)=A*[m(t)+C]*sin(2*pi*f*t)其中f=1000 Hz。
程序2_3.信号发生器的频率设置为100Hz ，Scope 的sampling time 改为0.00001s.（3）对子系统进行封装：请对（2）所建立所子系统进行封装（Mask ），编写参数输入对话框和帮助文档。
昆明理工大学MATLAB实验指导书(第二次实验)
************************ 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 。
《matlab原理与应用》实验指导书
《matlab原理与应用》实验指导书实验一 MATLAB集成环境使用实验性质：验证性实验级别：必做开课单位：信息与通信工程学院通信工程系学时：2学时一、实验目的：1、了解与MATLAB的运行环境相关的各类窗口及其相应的功能特点。
二、实验器材：安装matlab软件的计算机三、实验理论3.1 启动和退出M AT L A B不同的计算机系统，M AT L A B的启动也不一样。
在Windows 和M a c i n t o s h系统中，程序通常通过点击一个图标而启动。
在U N I X系统中，程序是通过在命令行系统提示符后键入如下字符启动：matlab如果上述工作有问题，可请教系统管理员。
当启动MATLAB时，如果matlab rc . m和startu p . m文件存在，则执行这些文件。
在这些文件中，为满足个人需要，用户可以给定命令以调整M AT L A B，例如，constants用于设置图形等。
在一个多用户系统上，系统管理员存储matlabrc . m文件，但你也能为自己的使用创建文件startup.m。
要退出MATLAB，键入quit或exit。
1、退出和中断e x i t，q u i t 结束M AT L A B会话。
程序完成，如果没有明确保存，则变量中的数据丢失。
Ctrl-c 中断一个MATLAB任务。
例如，当M AT L A B正在计算或打印时，中断一个任务，但会话并没有结束。
除此之外，对一些系统有指定的菜单选择。
例如，在Wi n d o w s和M a c i n t o s h系统中，在文件菜单下可以找到选项q u i t。
当编辑或执行M AT L A B时，下列的快捷键十分有用。
通常因为不同的平台使用不同的键，因此，给定了一些替换键。
在你的系统上试一下这些键，注意哪些键组合使用。
2、特殊的功能键↑或C t r l_p 恢复前面的命令。
↓或C t r l_n 恢复当前命令之后键入的命令。
- 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 是子图的编号。