MATLAB实验报告
《MATLAB/Simulink与控制系统仿真》实验报告
专业：
班级：
学号：
姓名：
指导教师：
实验1、MATLAB/Simulink 仿真基础及控制系统模型的建立
一、实验目的
1、掌握MATLAB/Simulink 仿真的基本知识；
2、熟练应用MATLAB 软件建立控制系统模型。 二、实验设备
电脑一台；MATLAB 仿真软件一个 三、实验内容
1、熟悉MATLAB/Smulink 仿真软件。
2、一个单位负反馈二阶系统，其开环传递函数为2
10()3G s s s
=
+。用Simulink 建立该控制系统模型，用示波器观察模型的阶跃响应曲线，并将阶跃响应曲线导入到MATLAB 的工作空间中，在命令窗口绘制该模型的阶跃响应曲线。
图 1系统结构图
图 2示波器输出结果图
3、某控制系统的传递函数为
()()()1()Y s G s X s G s =+，其中2
50()23s G s s s
+=+。用Simulink 建
立该控制系统模型，用示波器观察模型的阶跃响应曲线，并将阶跃响应曲线导入到MA TLAB 的工作空间中，在命令窗口绘制该模型的阶跃响应曲线。
图 3系统结构图图 4 示波器输出结果图
图 5 工作空间中仿真结果图形化输出
4、一闭环系统结构如图所示，其中系统前向通道的传递函数为
3
20.520()0.11220s G s s s s s
+=+++ ，而且前向通道有一个[-0.2，0.5]的限幅环节，图中用N 表示，反馈通道的增益为1.5，系统为负反馈，阶跃输入经1.5倍的增益作用到系统。用Simulink 建立该控制系统模型，用示波器观察模型的阶跃响应曲线，并将阶跃响应曲线导入到MATLAB 的工作空间中，在命令窗口绘制该模型的阶跃响应曲线。
图 6 系统结构图
图 7 示波器输出结果
实验2 MATLAB/Simulink 在控制系统建模中的应用
一、实验目的
1、掌握MATLAB/Simulink 在控制系统建模中的应用； 二、实验设备
电脑一台；MA TLAB 仿真软件一个 三、实验内容
1、给定RLC 网络如图所示。其中，为输入变量，为输出变量。求解这个系统的传递函数模型，零极点增益模型以及状态空间模型（假设，，，
）。
传递函数模型2
21
)(2++=
s s s G
程序代码如下：
clear all ;%清除工作空间的变量 num=[0,1]; %定义分子多项式 den=[1 2 2];%定义分母多项式
sy_tf=tf(num,den);%建立传递函数模型
[z,p,k]=tf2zp(num,den) %从传递函数模型获取系统的零极点增益 sy_zpk=zpk(z,p,k);%建立系统的零极点增益模型
[A,B,C,D]=zp2ss(z,p,k); %从零极点增益模型获取系统的状态空间模型 sys_ss=ss(A,B,C,D) %建立系统的状态空间模型 step(sy_tf) %求解系统的阶跃响应 grid on ； %添加栅格
程序运行结果
z =Empty matrix: 0-by-1 p =-1.0000 + 1.0000i -1.0000 - 1.0000i k =1
a = x1 x2
x1 -2 -1.414 x2 1.414 0 b = u1 x1 1
()i u t 0()u t 11R =Ω21R =Ω1C F =1L H
=
x2 0
c = x1 x2
y1 0 0.7071
d = u1
y1 0
Continuous-time model.
图 8 系统的阶跃响应曲线
2、已知某双环调速的电流环系统的结构图如图所示。试采用Simulink动态结构图求其线性模型。
图 9simulink中的系统动态模型
将图2模型存为“Samples_4_14.mdl”文件
在MA TLAB命令窗口运行以下命令，得到一个线性状态空间模型（A,B,C,D）。
[A,B,C,D]=linmod('Samples_4_14'); %提取simulink模型的状态空间模型
输出结果如下
A =1.0e+003 *
-0.0781 0 0 0 1.7964
0 -0.5000 0 0 0
0.0141 0 -0.5000 0 0
0 0.5000 -0.5000 0 0
0 0.1600 -0.1600 0.0250 -0.0599
B =0
1
C = 195.3125 0 0 0 0
D = 0
在MA TLAB命令窗口运行以下命令
[num,den]=ss2tf(A,B,C,D);%将状态空间模型转换为传递函数模型
pritfsys(num,den,'s'); %以传递函数模型形式显示出来
输出结果：
num/den = 4.5475e-013 s^4 + 5.8208e-011 s^3 + 56137724.5509 s^2 + 32454622005.9881 s + 2192879865269.464
-----------------------------------------------------------------------------
s^5 + 1138.0052 s^4 + 392683.3832 s^3 + 43221369.7605 s^2 + 3506268712.5749 s
+ 157887350299.4013
实验3 MATLAB/Simulink在时域分析法中的应用
一、实验目的
1、掌握时域分析中MA TLAB/Simulink函数的应用；
2、掌握MATLAB/Simulink在稳定性分析中的应用。
二、实验设备
电脑一台；MA TLAB仿真软件一个
三、实验内容
1、某随动系统的结构如图所示。利用MA TLAB完成如下工作：（1）对给定的随动系统建立数学模型；（2）分析系统的稳定性，并且绘制阶跃响应曲线；（3）计算系统的稳态误差；（4）大致分析系统的总体性能，并给出理论上的解释。
图 1 系统的结构图
解：利用ＭＡＴＬＡＢ求解的基本步骤如下
1．求取系统传递函数
clc;
clear all;
num1=[20];den1=[1 2 0];
sys1=tf(num1,den1); %二阶系统的传递函数
num2=[0.1 0];den2=[0 1];
sys2=tf(num2,den2); %微分环节传递函数
sys_inner=feedback(sys1,sys2); %内环反馈的传递函数
sys_outer=feedback(sys_inner,1)%外环反馈的传递函数
20
--------------
s^2 + 4 s + 20
得到系统的传递函数
2042^20
G(s)++=
s s
2．进行稳定性分析
den=[1 4 20];
roots(den) %求闭环系统特征多项式的根 pzmap(sys_outer); %利用pzmap 命令绘制系统的零极点图 grid on ;
程序运行结果 ans =
-2.0000 + 4.0000i -2.0000 - 4.0000i
由结果可知，系统特征根都具有负实部，因此闭环系统是稳定的。 系统零极点分布图如图2所示
图 2 系统零极点分布图
3．求阶跃响应
num=[20];den=[1 4 20];
[y,t,x]=step(num,den) %计算闭环系统的阶跃响应 plot(x,y); %绘制阶跃响应曲线 grid on ;
如下图3，横坐标表示响应时间，纵坐标表示系统输出
图3系统阶跃响应曲线
图4系统阶跃响应曲线
4．分析系统的响应特性
%计算系统的超调量
y_stable=1; %阶跃响应的稳态值
max_response=max(y); %闭环系统阶跃响应的最大值sigma=(max_response-y_stable) %阶跃响应的超调量
程序运行结果
sigma =0.2076
系统稳态误差为0，波形图如下
图 5 系统误差曲线图
%计算系统的上升时间
for i=1:length(y) %遍历响应曲线
if y(i)>y_stable %如果某个时刻系统的输出值大于稳态值 break ; %循环中断 end end
tr=x(i) %计算此时对应的时间，就是阶跃响应的上升时间 %计算系统的峰值时间
[max_response,index]=max(y); %查找系统阶跃响应的最大值
tp=x(index) %计算此时对应的时间，就是阶跃响应的峰值时间 %计算系统的调整时间---------》取误差带为2% for i=1:length(y) %遍历响应曲线
if max(y(i:length(y)))<=1.02*y_stable %如果当前响应值在误差带内 if min(y(i:length(y)))>=0.98*y_stable break ; %循环退出 end end end
ts=x(i) %计算此时对应的时间，就是阶跃响应的调整时间
程序运行结果 tr = 0.5245 tp = 0.7730 ts = 1.8773
即上升时间为0.52秒，峰值时间0.77秒，并且系统在经过1.88秒后进入稳态。 2、已知某二阶系统的传递函数为2
22)(n
n n
s s s G ωζωω++=
，（1）将自然频率固定为1=n ω，5,3,2,1,...,1.0,0=ζ，分析ζ变化时系统的单位阶跃响应；（2）将阻尼比ζ固定为55.0=ζ，
分析自然频率n ω变化时系统的阶跃响应（n ω变化范围为0.1~1）。
ω=1，ζ取不同（1）解：利用ＭＡＴＬＡＢ建立控制系统的数学模型，并且同时显示
n
值时的阶跃响应曲线
clc;
clear;
t=linspace(0,20,200);%设置仿真时间
omega=1;%设置二阶系统的自然频率
omega2=omega^2;%计算自然频率的平方
zuni=[0,0.1,0.2,0.5,1,2,3,5];
num=omega2;
for k=1:8
den=[1 2*zuni(k)*omega omega2];
sys=tf(num,den);
y(:,k)=step(sys,t);
end
figure(1);
plot(t,y(:,1:8));
grid;
gtext('zuni= 0');
gtext('zuni= 0.1');
gtext('zuni= 0.2');
gtext('zuni= 0.5');
gtext('zuni= 1');
gtext('zuni= 2');
gtext('zuni= 3');
gtext('zuni= 5');
图10 固定自然频率，阻尼比变化时系统的阶跃响应曲线
结论：当固定频率后，改变阻尼比，在ζ<1不会改变阶跃响应的震荡频率；而当ζ>1时，阶跃响应曲线不再震荡，系统阻尼。
ω从0.1变化到1是系统的阶跃响应曲线
（2）绘制ζ=0.55，
n
clc;
clear;
t=linspace(0,20,200);%设置仿真时间
zuni=0.55;%设定阻尼系数
omega=[0.1 0.2 0.4 0.7 1];%设定自然频率向量
omega2=omega.*2; %计算自然频率的平方
for k=1:5%循环五次，分别计算在五种不同的自然频率下系统的阶跃响应
num=omega2(k);
den=[1 2*zuni*omega(k) omega2(k)];
sys=tf(num,den);%系统传递函数
y(:,k)=step(sys,t); %计算当前自然频率下，二阶系统的阶跃响应值
end
figure(2);
plot(t,y(:,1:5)); %在一幅图像上依次绘出上述5条阶跃响应曲线
grid;
gtext('omega=0.1');
gtext('omega=0.2');
gtext('omega=0.4');
gtext('omega=0.7');
gtext('omega=1.0');
图 11 固定阻尼系数，自然频率变化时系统的阶跃响应曲线
ω从0.1变化到1时，系统震荡频率加快，上升时间减少，结论：当自然频率
n
过渡过程时间减少；系统响应更加迅速，动态性能更好。
ω决定了系统阶跃响应的震荡频率。
自然频率
n
实验4 MATLAB/Simulink 在根轨迹分析法中应用
一、实验目的
1、掌握MATLAB/Simulink 绘制根轨迹函数；
2、掌握MATLAB/Simulink 绘制根轨迹的方法。 二、实验设备
电脑一台；MA TLAB 仿真软件一个 三、实验内容
1、已知单位负反馈控制系统的开环传递函数)
4)(1()
1()(+-+=
s s s s k s G 。（1）画出这个
系统的根轨迹；（2）确定使闭环系统稳定的增益值k ；（3）分析系统的阶跃响应性能；（4）利用rltool 对系统的性能进行分析。
解：利用ＭＡＴＬＡＢ求解的基本步骤如下 1． 建立系统的数学模型
代码 clc; clear; num=[1 1];
den=conv([1 0],conv([1 -1],[1 4]));
sys=tf(num,den) %控制系统的开环传递函数模型 程序运行结果
Transfer function: s + 1 ----------------- s^3 + 3 s^2 - 4 s
结果输出是用来绘制部分根轨迹的那部分传递函数 ２．绘制根轨迹图
rlocus(sys); %绘制系统的根轨迹曲线
grid on ; title('根轨迹图 '); 得到系统根轨迹如图1
图6系统根轨迹图
利用“rlocfind“命令计算选定点处的增益和其他闭环极点。
[k,poles]=rlocfind(sys) %计算选定点处的增益和其他闭环极点
运行结果
k = 6.2809
poles = -2.9488
-0.0256 + 1.4592i
-0.0256 - 1.4592ii
当增益Ｋ>6时，闭环系统的极点都位于虚轴的左部，处于稳定。
3．使用rltool进行分析
执行命令“rltool(sys) “得到根轨迹分析界面图形，图2所示
图7控制系统根轨迹分析与设计器
利用rltool进行工具分析系统的阶跃响应。设定系统增益为20，可得到如图3的结果
分析，系统稳定，并且系统误差为0。系统的穿越频率为1.41，相角稳定裕度为17，剪切频率为1.38。
图8 K=20时系统的阶跃响应
实验5 MATLAB/Simulink在频域分析法中的应用
一、实验目的
1、掌握MATLAB绘制伯德图和乃奎斯特曲线；
2、熟练应用MATLAB分析稳定裕度。
二、实验设备
电脑一台；MA TLAB仿真软件一个
三、实验内容
1、已知晶闸管-直流电机开环系统结构图如图所示。试用Simulink动态结构图进行频域分析并求频域性能指标。
解：利用simulink求解步骤如下
步骤1 在simulink中建立系统动态模型，如下图
图表9 系统动态模型图
步骤2 求取系统的线性空间状态模型，并求取频域性能指标
在MA TLAB运行以下命令
[A,B,C,D]=linmod('Samples_7_9'); %提取simulin系统的线性空间状态模型
sys=ss(A,B,C,D);
margin(sys); %求取频域性能指标
图表10 系统的开环bode图和频域性能指标从图中可以看出：
幅值裕度GM=26.4dB，穿越频率为152rad/sec
相位裕度PM=54deg，穿越频率为25.5rad/sec
实验6 MATLAB_Simulink 在控制系统校正中的应用
一、实验目的
1、掌握建立控制系统的数学模型及设计系统的串联校正装置；
2、了解校正前后系统性能的比较。 二、实验设备
电脑一台；MA TLAB 仿真软件一个 三、实验内容
1、某单位负反馈控制系统的开环传递函数)
2)(1()(++=
s s s k
s G ，设计一个串联的校
正装置，使校正后的系统静态速度误差系数1
10-≥s ，相角裕度
45≥，增益裕量dB 10≥。
解：求解步骤如下
步骤1计算得系统开环传递函数)
2)(1(20
)(++=s s s s G
步骤2 建立控制系统的数学模型 代码如下
clc; clear;
num_open=[0 20];
den_open=conv(conv([1 0],[1 1]),[1 2]); sys_open=tf(num_open,den_open)
运行结果
Transfer function: 20 ----------------- s^3 + 3 s^2 + 2 s 步骤3 分析系统的动态特性
[Gm,Pm,Wcg,Wcp]=margin(sys_open) %计算相角裕量和增益裕量 margin(sys_open); 运行结果 Gm = 0.3000
Pm = -28.0814 Wcg = 1.4142 Wcp = 2.4253
频率响应特性曲线如下
图表 11 闭环系统的频率响应曲线
计算结果显示，未校正系统增益裕量只有—10.5，相角裕量为—28.0814，相角穿越频率为1.4142，幅值穿越频率为2.4253。系统尚不稳定需要串联校正环节。
步骤4 设计系统的串联校正装置
先设计止滞后环节。假定校正后的系统增益穿越频率为1，并且取零极点之比为10，则滞后环节的传递函数为
01
.0s 1
.0s ++
%设计串联校正器的滞后环节
num_zhihou=[1 0.1]; den_zhihou=[1 0.01];
sys_zhihou=tf(num_zhihou,den_zhihou); %滞后环节的传递函数模型
sys_new=sys_open*sys_zhihou %加入滞后环节后系统的开环传递函数 margin(sys_new); %绘制加入滞后环节后系统的Bode 曲线
图表 12 加入滞后环节的系统频率响应曲线
MATLAB实验报告
MATLAB程序设计语言 实 验 报 告 专业及班级：电子信息工程 姓名：王伟 学号：1107050322 日期 2013年6月20日
实验一 MATLAB 的基本使用 【一】 实验目的 1.了解MATALB 程序设计语言的基本特点，熟悉MATLAB 软件的运行环境； 2.掌握变量、函数等有关概念，掌握M 文件的创建、保存、打开的方法，初步具备将一般数学问题转化为对应计算机模型处理的能力； 3.掌握二维图形绘制的方法，并能用这些方法实现计算结果的可视化。 【二】 MATLAB 的基础知识 通过本课程的学习，应基本掌握以下的基础知识： 一. MATLAB 简介 二. MATLAB 的启动和退出 三. MATLAB 使用界面简介 四. 帮助信息的获取 五. MATLAB 的数值计算功能 六. 程序流程控制 七. M 文件 八. 函数文件 九. MATLAB 的可视化 【三】上机练习 1. 仔细预习第二部分内容，关于MATLAB 的基础知识。 2. 熟悉MATLAB 环境，将第二部分所有的例子在计算机上练习一遍 3. 已知矩阵???? ??????=??????????=123456789,987654321B A 。求A*B ，A .* B ，比较二者结果是否相同。并利用MATLAB 的内部函数求矩阵A 的大小、元素和、长度以 及最大值。 程序代码： >> A=[1 2 3;4 5 6;7 8 9]; >> B=[9 8 7;6 5 4;3 2 1]; >> A*B ans =
30 24 18 84 69 54 138 114 90 >> A.*B ans = 9 16 21 24 25 24 21 16 9 两者结果不同 >> [m,n]=size(A) m = 3 n = 3 >> b=sum(A) b = 12 15 18 >> a=length(A) a = 3 >>max(A) ans =
MATLAB实验报告50059
实验一MATLAB操作基础 实验目的和要求： 1、熟悉MATLAB的操作环境及基本操作方法。 2、掌握MATLAB的搜索路径及设置方法。 3、熟悉MATLAB帮助信息的查阅方法 实验内容： 1、建立自己的工作目录，再设置自己的工作目录设置到MA TLAB搜索路径下，再试 验用help命令能否查询到自己的工作目录。 2、在MA TLAB的操作环境下验证课本；例1-1至例1-4，总结MATLAB的特点。 例1-1
例1-2 例1-3 例1-4
3、利用帮助功能查询inv、plot、max、round等函数的功能。 4、完成下列操作： （1）在matlab命令窗口输入以下命令： x=0:pi/10:2*pi; y=sin(x); （2）在工作空间窗口选择变量y，再在工作空间窗口选择回绘图菜单命令或在工具栏中单击绘图命令按钮，绘制变量y的图形，并分析图形的含义。
5、访问mathworks公司的主页，查询有关MATLAB的产品信息。 主要教学环节的组织： 教师讲授实验目的、开发环境界面、演示实验过程，然后同学上机练习。 思考题： 1、如何启动与退出MA TLAB集成环境？ 启动： （1）在windows桌面，单击任务栏上的开始按钮，选择‘所有程序’菜单项，然后选择MA TLAB程序组中的MA TLABR2008b程序选项，即可启动 MATLAB系统。 （2）在MA TLAB的安装路径中找到MA TLAB系统启动程序matlab.exe，然后运行它。 （3）在桌面上建立快捷方式后。双击快捷方式图标，启动MA TLAB。 退出： （1）在MA TLAB主窗口file菜单中选择exitMATLAB命令。 （2）在MA TLAB命令窗口中输入exit或quit命令。 （3）单击MATLAB主窗口的关闭按钮。 2、简述MATLAB的主要功能。 MATLAB是一种应用于科学计算领域的数学软件，它主要包括数值计算和符 号计算功能、绘图功能、编程语言功能以及应用工具箱的扩展功能。 3、如果一个MATLAB命令包含的字符很多，需要分成多行输入，该如何处理？
matlab实验报告
MATLAB 数学实验报告 指导老师： 班级： 小组成员： 时间：201_/_/_
Matlab 第二次实验报告 小组成员： 1 题目：实验四，MATLAB 选择结构与应用实验 目的：掌握if 选择结构与程序流程控制，重点掌握break，return , pause语句的应用。 问题：问题1：验证“哥德巴赫猜想” ，即：任何一个正偶数（n>=6）均可表示为两个质数的和。要求编制一个函数程序，输入一个正偶数，返回两个质数的和。 问题分析：由用户输入一个大于6 的偶数，由input 语句实现。由if 判断语句判断是否输入的数据符合条件。再引用质数判断函数来找出两个质数，再向屏幕输出两个质数即可。 编程：function [z1,z2]=gede(n); n=input('please input n')
if n<6 disp('data error'); return end if mod(n,2)==0 for i=2:n/2 k=0; for j=2:sqrt(i) if mod(i,j)==0 k=k+1; end end for j=2:sqrt(n-i) if mod(n-i,j)==0 k=k+1; end end if k==0 fprintf('two numbers are') fprintf('%.0f,%.0f',i,n- i) break end
end end 结果分析 如上图，用户输入了大于6的偶数返回两个质数5和31,通过 不断试验，即可验证哥德巴赫猜想。 纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰，更快的解决问题。 2题目：实验四，MATLAB选择结构与应用实验 目的：用matlab联系生活实际，解决一些生活中常见的实际问 题。
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告 一、实验目的 1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。 2.学会运用MATLAB表示常用连续时间信号的方法 3.观察并熟悉一些信号的波形和特性。 4.学会运用MATLAB进行连续信号时移、反折和尺度变换。 5.学会运用MATLAB进行连续时间微分、积分运算。 6.学会运用MATLAB进行连续信号相加、相乘运算。 7.学会运用MATLAB进行连续信号的奇偶分解。 二、实验任务 将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。 三、实验内容 1.MATLAB软件基本运算入门。 1). MATLAB软件的数值计算： 算数运算 向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。 矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开； 矩阵的不同行之间必须用分号”;”或者ENTER分开。2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。 2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名” 2.MATLAB软件简单二维图形绘制 1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y) 2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p 表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p) 3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin]) 4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’) 5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’) 6).输出：grid on 举例1：
MATLAB全实验报告
《数学实验》报告 实验名称 Matlab 基础知识 学院 专业班级 姓名 学号 2014年 6月
一、【实验目的】 1.认识熟悉Matlab这一软件，并在此基础上学会基本操作。 2.掌握Matlab基本操作和常用命令。 3.了解Matlab常用函数，运算符和表达式。 4.掌握Matlab工作方式和M文件的相关知识。 5.学会Matlab中矩阵和数组的运算。 二、【实验任务】 P16 第4题 编写函数文件，计算 1! n k k = ∑，并求出当k=20时表达式的值。P27第2题 矩阵A= 123 456 789 ?? ?? ?? ?? ?? ，B= 468 556 322 ?? ?? ?? ?? ?? ，计算A*B,A.*B,并比较两者的区别。 P27第3题 已知矩阵A= 52 91 ?? ?? ?? ，B= 12 92 ?? ?? ?? ，做简单的关系运算A>B,A==B,A
P27第3题 >> A=[5 2;9 1];B=[1 2;9 2]; >>A>B >>A==B >>A> (A==B)&(A
Matlab数学实验报告一
数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦
实验目的：用二分法与Newton迭代法求解非线性方程的根； 用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 编程实现二分法及Newton迭代法； 学会使用Matlab函数solve、fzero、fsolve求解非线性方程（组）的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容： 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方） (1)在区间（0，1)内用二分法； (2)用迭代法x=(2-exp(x))/10，取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0，则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正，则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负，则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法； format long
syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果： （1）二分法： symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 （2）迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次，而用迭代法求得同样精度的解仅用5次，但由于迭代法一般只具有局部收敛性，因此通常不用二分法来求得非线性方程的精确解，而只用它求得根的一个近似解，再用收敛速度较快的迭代法求得其精确解。
MATLAB实验报告实验二
实验二 MATLAB矩阵及其运算 学号：3121003104 姓名：刘艳琳专业：电子信息工程1班日期：2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. （1）新建一个.m文件，验证书本第15页例2-1； （2）用命令方式查看和保存代码中的所有变量；
（3）用命令方式删除所有变量； （4）用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式
长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 （1）w=sqrt(2)*(1+0.34245*10^(-6)) （2）x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; （3）y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; （4）z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];
4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20]，对其进行如下操作：（1）输出A在[ 7, 10]范围内的全部元素； （2）取出A的第2，4行和第1，3，5列； （3）对矩阵A变换成向量B，B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]； （4）删除A的第2，3，4行元素； (1) (2)
数学实验“几种常见的求积分近似解的方法”实验报告(内含matlab程序)
西京学院数学软件实验任务书
实验二十一实验报告 一、实验名称：Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 二、实验目的：进一步熟悉Romberg 积分法，Gauss 型积分法，高斯-勒让德积分法，高斯-切比雪夫积分法，高斯-拉盖尔积分法，高斯-埃尔米特积分法。 三、实验要求：运用Matlab/C/C++/Java/Maple/Mathematica 等其中一种语言完成程序设计。 四、实验原理： 1．Romberg 积分法： 龙贝格积分法是用里查森外推算法来加快复合梯形求积公式的收敛速度，它的算法如下，其中()i m T 是通过一系列逼近原定积分的龙贝格分值. 计算(0)1[()()]2 b a T f a f b -= + 对1,2,3,k n = ，计算下列各步： 21()(1)1 111 1(21)()[()]222k k k k k j b a j b a T T f a ---=---=++∑
对1,2,,m k = 和,1,2,,1i k k k =-- ，计算111 441 m i i i m m m m T T T --+-=- 随着计算的步骤的增加，()i m T 越来越逼近积分()b a f x dx ?。 2．Gauss 型积分法： 高斯积分公式的思想是用n 个不等距的节点123,,,n x x x x 对被积函数进行插值，然后对插值后的函数进行积分，其积分公式为： 1 1 1 ()()n k k k f x dx A f x -=≈∑? 如果积分区间不是[1,1]-，则需转换到此区间： 11()()222 b a b a b a b a f x dx f t dt ---+= +? ? 其中系数k A 、节点k x 与n 的关系如下表所示： 3．高斯-切比雪夫积分法： 第一类切比雪夫积分形式为： 1 1 ()()n k k k f x dx A f x -=≈∑? 其中k A n π= ，21cos 2k k x n π-= 4．高斯-拉盖尔积分法： 高斯-拉盖尔公式有两种形式： 1 ()()n x k k k e f x dx A f x +∞ -=≈∑?
matlab数学实验报告5
数学实验报告 制作成员班级学号 2011年6月12日
培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度
差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x )
matlab实验报告
2015秋2013级《MATLAB程序设计》实验报告 实验一班级：软件131姓名：陈万全学号：132852 一、实验目的 1、了解MATLAB程序设计的开发环境，熟悉命令窗口、工作区窗口、历史命令等窗口的使用。 2、掌握MATLAB常用命令的使用。 3、掌握MATLAB帮助系统的使用。 4、熟悉利用MATLAB进行简单数学计算以及绘图的操作方法。 二、实验内容 1、启动MATLAB软件，熟悉MATLAB的基本工作桌面，了解各个窗口的功能与使用。 图1 MATLAB工作桌面 2、MATLAB的常用命令与系统帮助: (1)系统帮助 help：用来查询已知命令的用法。例如已知inv是用来计算逆矩阵，键入help inv即可得知有关inv命令的用法。 lookfor：用来寻找未知的命令。例如要寻找计算反矩阵的命令，可键入lookfor inverse，MATLAB即会列出所有和关键字inverse相关的指令。找到所需的
命令後，即可用help进一步找出其用法。 (2)数据显示格式: 常用命令：说明 format short 显示小数点后4位（缺省值） format long 显示15位 format bank 显示小数点后2位 format + 显示+，-，0 format short e 5位科学记数法 format long e 15位科学记数法 format rat 最接近的有理数显示 (3)命令行编辑：键盘上的各种箭头和控制键提供了命令的重调、编辑功能。 具体用法如下： ↑----重调前一行（可重复使用调用更早的） ↓----重调后一行 →----前移一字符 ←----后移一字符 home----前移到行首 end----移动到行末 esc----清除一行 del----清除当前字符 backspace----清除前一字符 (4)MATLAB工作区常用命令： who--------显示当前工作区中所有用户变量名 whos--------显示当前工作区中所有用户变量名及大小、字节数和类型 disp(x) -----显示变量X的内容 clear -----清除工作区中用户定义的所有变量 save文件名-----保存工作区中用户定义的所有变量到指定文件中 load文件名-----载入指定文件中的数据
MATLAB实验报告
数字信号处理及MATLAB 实验报告 班级： 学号： 姓名：
4.7.2 例4,2 设x(n)是由两个正弦信号及白噪声的叠加，试用FFT文件对其作频谱分析。程序清单 %产生两个正弦加白噪声 N=256; f1=.1;f2=.2;fs=1; a1=5;a2=3; w=2*pi/fs; x=a1*sin(w*f1*(0:N-1))+a2*sin(w*f2*(0:N-1))+randn(1,N); %应用FFT求频谱 subplot(2,2,1); plot(x(1:N/4)); title('原始信号'); f=-0.5:1/N:0.5-1/N; x=fft(x); y=ifft(x); subplot(2,2,2); plot(f,fftshift(abs(x))); title('频域信号'); subplot(2,2,3); plot(real(x(1:N/4))); title('时域信号');
例4.3 设x（n）为长度N=6的矩形序列，用MATLAB程序分析FFT取不同长度时x（n）频谱的变化。N=8,32,64,时x（n）的FFT MATLAB实现程序如下。 x=[1,1,1,1,1,1]; N=8; y1=fft(x,N); n=0:N-1; subplot(3,1,1);stem(n,abs(y1),'.k');axis([0,9,0,6]); N=32; y2=fft(x,N); n=0:N-1; subplot(3,1,2);stem(n,abs(y2),'.k');axis([0,40,0,6]); N=64; y3=fft(x,N); subplot(3,1,3);stem(n,abs(y3),'.k');axis([0,80,0,6]);
MATLAB实验报告
实验一 MATLAB 环境的熟悉与基本运算 一、实验目的及要求 1．熟悉MATLAB 的开发环境； 2．掌握MATLAB 的一些常用命令； 3．掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容 1.熟悉MATLAB 的开发环境： ① MATLAB 的各种窗口： 命令窗口、命令历史窗口、工作空间窗口、当前路径窗口。 ②路径的设置： 建立自己的文件夹，加入到MATLAB 路径中，并保存。 设置当前路径，以方便文件管理。 2.学习使用clc 、clear ，了解其功能和作用。 3.矩阵运算： 已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A.*B ，并比较结果。 4.使用冒号选出指定元素: 已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素；A 中所有列第2，3行的元素； 5.在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ 6.关系及逻辑运算 1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15]，求: y=a==b ，并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0]，求: x&y+x>y ，并分析结果 7.文件操作 1)将0到1000的所有整数，写入到D 盘下的文件 2)读入D 盘下的文件，并赋给变量num
8.符号运算 1)对表达式f=x 3 -1 进行因式分解 2)对表达式f=（2x 2*(x+3)-10)*t ，分别将自变量x 和t 的同类项合并 3）求 3(1)x dz z +? 三、实验报告要求 完成实验内容的3、4、5、6、7、8，写出相应的程序、结果
matlab实验心得体会.doc
matlab实验心得体会 篇一：matlab实验心得总结 通过《matlab仿真》实验使我学习掌握了许多知识。首先是对matlab 有了一个全新的认识，其次是对matlab的更多操作和命令的使用有了更高的掌握，最重要的事对matlab的处理能力有了一个更高的飞跃尤其是对相关函数的使用及相关问题的处理。 就对matlab相关的命令操作而言，通过这次实验的亲身操作和实践，学习掌握了许多原本不知道的或者不太熟悉的命令。比如说相关m文件的建立，画图用到的标注，配色，坐标控制，同一张图里画几幅不同的图像，相关参数的设置以及相关函数的调用格式等等。就拿建立一个数学方程而言，通过设置不同的参数达到所需要的要求和结果，而且还可以在不同的窗口建立不同的函数而达到相同的效果，比如说可以再命令窗口和m文件中通过不同的命令设置的到相同的所需的效果图。而自己对于矩阵及闭环传递函数的建立原本所掌握的知识几乎为零，而通过这次实验使我彻底的掌握了相关的命令操作和处理的方法，在这里我们不仅可以通过建立函数和参数来达到目标效果，而且还可以通过可视化的编程达到更快更方便，更简洁的效果。就拿可视化编程而言原本根本就只是听说而已罢了，从来就没有亲身去尝试过，然而现在自己却可以和容易的通过搭建不同功能木块来实现相关的函数及功能。这些在原本根本就不敢相信，然而通过《matlab仿真》的学习和实验亲身操作这些原本看似不可能的操作在此就变的轻而易举的事了。
再此我不得不题到的事指导老师教我们怎么去搭建构造相关闭环传递函数的实验，这个实验几乎在我们的这次实验中占据了非常大的比重，在后面的几个大一点的实验中几乎都是涉及这个方面的内容，我现在想说的事怎么去搭建相关的函数和功能模块对我们来说几乎已经不是什么难事了，就拿怎么去对模块功能的实现以及分析确实是个重点和难点。通过对同一个模块分析其对应的不同的参数分析图的建立去分析和解释其对应的相关功能和技术指标和性能分析是非常重要的，我们不可能只需要建立相关的模块和功能就说自己掌握了所有的相关知识和技术，真正的技术和知识是怎么去分析和解释相关的技术指标和功能参数才是重中之重。就此而言，我坦诚的说自己所掌握的还是十分的有限的，但是老师给我们介绍的相关方法和技巧还是十分有效果的，如果自己真的想在这方面有什么建树对自己以后的要求还是需要更改的要求的，万不可以就此止步不前，自命不凡，我们还需掌握和了解还有许多许多，我们真正所掌握的只是皮毛，要想取得更大的成绩就得不断的去努力学习和汲取相关的知识和技巧。万不可自以为傲，裹足不前，matlab真的是个非常强大和有用的工具我们真正的能把它学懂学透的话还是需要下非常大的功夫和努力的。然而，不是说兴趣才是最大的老师嘛，我也相信，只要你自己有兴趣，即使它再怎么强大和难搞，我们能做的还是非常多的，关键的就只是在于你自己的态度了。我这里想说的事，matlab对于我来说是非常有吸引力的，我不敢说自己多么喜欢它，但是兴趣确实蛮高的，所以我相信在以后的学习和工作当中matlab将成为我非常有用的帮忙工具和好伙伴，也许这要说有点太草率了，但是我觉得对它的评价怎么也不会让每一个接触
浅析Matlab数学实验报告
数学实验报告 姓名: 班级: 学号: 第一次实验任务 过程: a=1+3i; b=2-i; 结果: a+b =3.0000 + 2.0000i a-b =-1.0000 + 4.0000i a*b = 5.0000 + 5.0000i a/b = -0.2000 + 1.4000i 过程: x=-4.5*pi/180; y=7.6*pi/180; 结果: sin(abs(x)+y)/sqrt(cos(abs(x+y))) =0.2098 心得:对于matlab 中的角度计算应转为弧度。 (1)过程: x=0:0.01:2*pi; y1=sin(x); y2=cos(x); y3=exp(x); y4=log(x); plot(x,y1,x,y2,x,y3,x,y4) plot(x,y1,x,y2,x,y3,x,y4) 结果: (2)过程:>> subplot(2,2,1) >> plot(x,y1) >> subplot(2,2,2) >> plot(x,y2) ./,,,,2,311b a b a b a b a i b i a ?-+-=+=计算、设有两个复数 6,7,5.4)
cos()sin(2=-=++y x y x y x ,其中、计算的图形。 下分别绘制)同一页面四个坐标系)同一坐标系下(、在( x y e y x y x y x ln ,,cos ,sin 213==== >> subplot(2,2,3) >> plot(x,y3) >> subplot(2.2.4) >> subplot(2,2,4) >> plot(x,y4) 结果: 心得:在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与hold on进行区别,后者为在同一坐标系中划出多条曲线。 5、随机生成一个3x3矩阵A及3x2矩阵B,计算(1)AB,(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。 过程:A=fix(rand(3,3).*10) ; B=fix(rand(3,3).*10);
MATLAB课程设计实验体会
课程设计实验体会 学生姓名：李祥胜 学生学号：20120704 专业班级：光信息科学与技术 指导老师：miss Chen 学院：信息工程学院 题目: MATLAB学期实验总结
MATLAB概念及介绍 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室）。是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB集成环境主要包括五个部分：MATLAB语言、MATLAB工作环境、句柄图形、MATLAB数学函数库和数学建模、小波分析、MATLAB API(App lication Program Interface)。MATLAB语言是以数组为基本数据单位，包括控制流程语句、函数、数据结构、输人输出及面向对象等特点的高级语言。利用SIMULINK对系统进行仿真与分析，在进入虚拟实验环境后，不需要书写代码，只需使用鼠标拖动库中的功能模块并将它们连接起来，再按照实验要求修改各元器件的参数。通过虚拟实验环境建立实验仿真电路模型，可使一些枯燥的电路变得有趣味，复杂的波形变得形象生动，使得各种复杂的能量转换过程比较直观地呈现。 1.1、MATLAB语言特点及优势 1.1.1、语言特点 MATLAB被称为第四代计算机语言，利用其丰富的函数资源，使编程人员从繁琐的程序代码中解放出来。MATLAB的最突出的特点就是简洁。MATLAB用更直观的、符合人们思维习惯的代码，代替了C和FORTRAN语言的冗长代码。MATLAB 给用户带来的是最直观、最简洁的程序开发环境。以下简单介绍一下MATLAB的主要特点。 （1）语言简洁紧凑，使用方便灵活，库函数极其丰富。MATLAB程序书写形式自由，利用其丰富的库函数避开繁杂的子程序编程任务，压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写，用户不必担心函数的可靠性。 （2)运算符丰富。由于MATLAB是用C语言编写的，MATLAB提供了和C语言几乎一样多的运算符，灵活使用MATLAB的运算符将使程序变得极为简短，具体运算符见附表。 （3）MATLAB既具有结构化的控制语句（如for循环、while循环、break语句和if语句），又有面向对象编程的特性。 （4）语法限制不严格，程序设计自由度大。例如，在MATLAB里，用户无需对矩阵预定义就可使用。 （5）程序的可移植性很好，基本上不做修改就可以在各种型号的计算机和操作系统上运行。
matlab实验心得总结
通过《matlab仿真》实验使我学习掌握了许多知识。首先是对matlab有了一个全新的认识，其次是对matlab的更多操作和命令的使用有了更高的掌握，最重要的事对matlab的处理能力有了一个更高的飞跃尤其是对相关函数的使用及相关问题的处理。 就对matlab相关的命令操作而言，通过这次实验的亲身操作和实践，学习掌握了许多原本不知道的或者不太熟悉的命令。比如说相关m文件的建立，画图用到的标注，配色，坐标控制，同一张图里画几幅不同的图像，相关参数的设置以及相关函数的调用格式等等。就拿建立一个数学方程而言，通过设置不同的参数达到所需要的要求和结果，而且还可以在不同的窗口建立不同的函数而达到相同的效果，比如说可以再命令窗口和m文件中通过不同的命令设置的到相同的所需的效果图。而自己对于矩阵及闭环传递函数的建立原本所掌握的知识几乎为零，而通过这次实验使我彻底的掌握了相关的命令操作和处理的方法，在这里我们不仅可以通过建立函数和参数来达到目标效果，而且还可以通过可视化的编程达到更快更方便，更简洁的效果。就拿可视化编程而言原本根本就只是听说而已罢了，从来就没有亲身去尝试过，然而现在自己却可以和容易的通过搭建不同功能木块来实现相关的函数及功能。这些在原本根本就不敢相信，然而通过《matlab仿真》的学习和实验亲身操作这些原本看似不可能的操作在此就变的轻而易举的事了。 再此我不得不题到的事指导老师教我们怎么去搭建构造相关闭环传递函数的实验，这个实验几乎在我们的这次实验中占据了非常大的比重，在后面的几个大一点的实验中几乎都是涉及这个方面的内容，我现在想说的事怎么去搭建相关的函数和功能模块对我们来说几乎已经不是什么难事了，就拿怎么去对模块功能的实现以及分析确实是个重点和难点。通过对同一个模块分析其对应的不同的参数分析图的建立去分析和解释其对应的相关功能和技术指标和性能分析是非常重要的，我们不可能只需要建立相关的模块和功能就说自己掌握了所有的相关知识和技术，真正的技术和知识是怎么去分析和解释相关的技术指标和功能参数才是重中之重。就此而言，我坦诚的说自己所掌握的还是十分的有限的，但是老师给我们介绍的相关方法和技巧还是十分有效果的，如果自己真的想在这方面有什么建树对自己以后的要求还是需要更改的要求的，万不可以就此止步不前，自命不凡，我们还需掌握和了解还有许多许多，我们真正所掌握的只是皮毛，要想取得更大的成绩就得不断的去努力学习和汲取相关的知识和技巧。万不可自以为傲，裹足不前，matlab真的是个非常强大和有用的工具我们真正的能把它学懂学透的话还是需要下非常大的功夫和努力的。然而，不是说兴趣才是最大的老师嘛，我也相信，只要你自己有兴趣，即使它再怎么强大和难搞，我们能做的还是非常多的，关键的就只是在于你自己的态度了。我这里想说的事，matlab对于我来说是非常有吸引力的，我不敢说自己多么喜欢它，但是兴趣确实蛮高的，所以我相信在以后的学习和工作当中matlab将成为我非常有用的帮忙工具和好伙伴，也许这要说有点太草率了，但是我觉得对它的评价怎么也不会让每一个接触过它的人吝啬是自己赞美之词。它不仅仅可以用来建模分析函数，还可以用来进行图形的建模和仿真，还可以用来分析系统和函数的参数稳定性等等。再次就不一一列举了，我怕三天也不会说完的。
《MATLAB及应用》实验报告2
核科学技术学院 实验报告 实验项目名称MATLAB符号计算 所属课程名称MATLAB及应用 实验类型上机实验 实验日期12月日 指导教师谢芹 班级 学号 姓名 成绩 一、实验名称 MATLAB符号计算 二、实验目的
（1）掌握定义符号对象的方法 （2）掌握符号表达式的运算法则以及符号矩阵运算 （3）掌握求符号函数极限及导数的方法 （4）掌握求符号函数定积分和不定积分的方法 三、实验原理 1. 函数极限及导数的方法 （1）函数极限：limit(F,x,a) 求符号函数f(x)的极限值。即计算当变量x趋近于常数a时，f(x)函数的极限值。 （2）limit(f)：求符号函数f(x)的极限值。符号函数f(x)的变量为函数findsym(f)确定的默认变量；没有指定变量的目标值时，系统默认变量趋近于0，即a=0的情况。 （3）limit(f,x,a,'right')：求符号函数f的极限值。'right'表示变量x从右边趋近于a。 （4）limit(f,x,a,‘left’)：求符号函数f的极限值。‘left’表示变量x从左边趋近于a。 2. 微分： diff(s)：没有指定变量和导数阶数，则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。 diff(s,'v')：以v为自变量，对符号表达式s求一阶导数。 diff(s,n)：按findsym函数指示的默认变量对符号表达式s求n阶导数，n为正整数。 diff(s,'v',n)：以v为自变量，对符号表达式s求n阶导数。 3. 函数定积分和不定积分的方法： int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。 int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分。 int(s,v,a,b)：求定积分运算。a,b分别表示定积分的下限和上限。 梯形法：trapz(x,y)：x为分割点构成的向量，y为被积函数在分割点上的函数值构成的向量； 抛物线法：quad(f,a,b,tol)，f是被积函数，[a,b]是积分区间，tol是精度。
matlab实验报告
Matlab实验报告 ——定积分的近似计算 学生姓名： 学号： 专业：数学与应用数学专业
数学实验报告 实验序号:1001114030 日期：2012年10月20日 班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算 问题背景描述： 利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到， 这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没 有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只 能应用近似方法去计算相应的定积分。 实验目的： 本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。 实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为 相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且
为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数 值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。 实验所用软件及版本：Matlab 7.0.1