发布时间 : 星期四 文章MATLAB实验指导书更新完毕开始阅读d8deca86240c844768eaee14
f(x1)*f(xm)<0，则解存在的区间缩小为[x1,xm]，否则解存在的区间缩小为[xm,x2]。重复这样的步骤，直到区间的长度小于一个可以接受的小数（比如1e-10）,则认为中点即是原方程的解。
实验四 函数的编写及调试
一、实验目的：掌握MATLAB函数的编写及调试方法。 二、实验内容：
1、编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。
function [y]=myfun1(x)
?sinx,?y(x)??x,??x?6,?x?00?x?3x?3
选择一些数据测试你编写的函数。
2、编写一个函数求向量x中元素的平均值、最大值、最小值、均方根值。
function [m_x,max_x,min_x,rms_x]=myfun2(x)方均根值(Root
Mean Square)的计算公式为：
rms?1N2x?ii?1N
用下面数据测试你写的函数：
（1）x=sin(0:0.01:6*pi)
（2）x=rand(1,200)，得到的x为200个（0，1）之间均匀分布的随机数。
3、编写一个函数，给出一个向量x?[x1,x2,?xn]，生成如下范德蒙矩阵。
1?1?xx2?1function [v]=myvander(x)?x12x22????n?1n?1x2??x1例如：>>v=myvander([2 3 4 5]) 得v=
?????1?xn??2xn? ???n?1xn??1??111?234?5?? ?491625???82764125??生成一些数据测试你写的函数。
三、思考题
4
编写程序，用如下迭代公式求
a，a的值分别为：3，17，113。迭代的终止条件为
100次。分别对迭代结果和迭
代
次
数
。
xn?1?xn?10?5，迭代初值x0?1.0,迭代次数不超过
准
2确值进
2行
2比较，并统计
xn?1ax?a?2x?x?a?xn?1??22xn?1
实验五 MATLAB的绘图
1、在同一坐标系下绘制下面三个函数在t?[0，4?]的图象。
y1?ty2?t
y3?4?e?0.1tsin(t)2、编写程序，选择合适的步距，绘制下面函数在区间[-6，6]中的图象。
?sinx,?y(x)??x,??x?6,?x?00?x?3x?3
3、用compass函数画下面相量图
ua = 1 ; ub = cos(-2*pi/3)+sin(-2*pi/3)*i ; uc=cos(2*pi/3)+sin(2*pi/3)*i; compass([ua,ub,uc,ua-ub,ub-uc,uc-ua]) 90120 1.5150 1 0.530 2601800210330240270300 4、三维空间曲线绘制z=0:0.1:4*pi;
x=cos(z); y=sin(z);
5
plot3(x,y,z)
5、用mesh或surf函数，绘制下面方程所表示的三维空间曲面，x和y的取值范围设为[-3，3]。
x2y2z???
1010三、思考题
在同一坐标系下，用不同颜色和线型绘制以下两个函数在t?[-2?，2?]范围内的图象。
y1?2
0.5t
y2?2e?0.2t
实验六 MATLAB数值运算
一、实验目的：掌握MATLAB常用的数值运算函数。 二、实验内容：
1、求代数方程3x?4x?7x?2x?9x?12?0的5个根，并将其用星号（*）标记在复平面图上。（用roots和plot函数）。
2、求代数方程x?1?0的5个根，并将其用星号（*）标记在复平面图上。（用roots和plot函数）。
3、求下面函数在[0.5,4]区间内的过零点。（用fzero函）
55432f(x)?x3?2x2sin(x)?5xcos(x)?4、
1x
已知R=50欧姆，U=4V，二极管D正向电流与电压的关系为：
Id?Ise?Udq???KT?1????
其中：
Ud为二极管正向电压
Is为反向饱合电流，取10-12A K为玻尔茨曼常数，1.38*10-23
T为绝对温度，取300开尔文（27摄氏度） q为电子电荷1.6*10-19C
6
求此电路中的电流Id和二极管正向电压Ud（要求用fsolve函数求解）
5、实验数据处理：已知某压力传感器的测试数据如下表 p 0.0 1.1 2.1 2.8 4.2 5.0 6.1 6.9 u 10 11 13 14 17 18 22 24 p为压力值，u为电压值，试用多项式u(p)8.1 29 9.0 34 9.9 39 ?ap3?bp2?cp?d来拟合其特
性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上。
实验七 MATLAB应用
1、以原点为奇对称中心的方波y(wt)，可以用相应频率的基波及其奇次谐波合成。
?4?111sinwt?sin3wt?sin5wt???sin(2n?1)wt? ??35(2n?1)??n?1,2,3,?y(wt)?取的阶数越多，越接近方波，但总消除不了边缘上的尖峰，这称为吉布斯效应。设方波频率为50Hz,时间t取0~0.04秒(f=50Hz,w=2*pi*f,h=1e-5,tf=40e-3,t = 0:h:tf)，编写程序，画出如下用1次谐波、1,3次谐波、1,3,5,7,9次谐波，1,3,5,…,19次谐波合成的近似方波。(产生方波的函数为:square)
2、用Simulink求解下图所示电路0~100微秒内的响应。已知R=6*10-4欧，C=1700微法，L=6*10-9享，uc(0)=15kV。
RGain21s1sUcUciScopeUL1/LGain1ii-1/CGainIntegratorIntegrator1Scope1
? 模块参数设置:
Integrator1的Initial condition:15kV
在命令窗口为R,L,C赋值。 ? 仿真参数设置如下：
Start time:0
Stop time:100e-6
7