Matlab实验[皇泽明]

自动控制原理Matlab仿真实验

姓名：皇泽明

班级：电1

学号：0953305006

完成时间：2011-12-3

实验一 MATLAB 及仿真实验（控制系统的时域分析）

学习利用MATLAB 进行以下实验，要求熟练掌握实验内容中所用到的指令，并按内容要求完成实验。 一、实验目的

学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点

1、系统的典型响应有哪些？

2、如何判断系统稳定性？

3、系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应

1、 阶跃响应：

阶跃响应常用格式：

1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。

2、),(Tn sys step ；表示时间范围0---Tn 。

3、),(T sys step ；表示时间范围向量T 指定。

4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应：

脉冲函数在数学上的精确定义：0

,0)(1)(0?==?∞

t x f dx x f

其拉氏变换为：

)

()()()(1

)(s G s f s G s Y s f ===

所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ②

);

,();

,(T sys impulse Tn sys impulse

③ ),(T sys impulse Y =

（二） 分析系统稳定性

有以下三种方法：

1、利用pzmap 绘制连续系统的零极点图；

2、利用tf2zp 求出系统零极点；

3、利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析

Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.

四、实验内容 (一) 稳定性

1.系统传函为()2

72436

45232

345234+++++++++=s s s s s s s s s s G ，试判断其稳定性 测试程序：den=[1,3,4,2,7,2];p=roots(den) 测试结果： p =

-1.7680 + 1.2673i -1.7680 - 1.2673i 0.4176 + 1.1130i 0.4176 - 1.1130i -0.2991 系统不稳定。

2.用Matlab 求出2

53722)(2342++++++=s s s s s s s G 的极点。 测试程序：den=[1,7,3,5,2];p=roots(den)

测试结果： p =

-6.6553 0.0327 + 0.8555i 0.0327 - 0.8555i -0.4100 系统不稳定。 （二）阶跃响应 1. 二阶系统()10

2102

++=

s s s G 1）键入程序，观察并记录单位阶跃响应曲线

2）计算系统的闭环根、阻尼比、无阻尼振荡频率，并记录

3）记录实际测取的峰值大小、峰值时间及过渡过程时间，并填表：

1测试程序：

num=10;den=[1,2,10];sys=tf(num,den);step(sys);title('G(s)=10/(s^2+2s+10)单位阶跃响应曲线') 测试结果：

012

3

456

0.2

0.4

0.6

0.8

1

1.2

1.4

G(s)=10/(s 2+2s+10)单位阶跃响应曲线

Time (sec)

A m p l i t u d e

2测试程序：num=[10];den=[1,2,10];sys=tf(num,den); damp(sys) 测试结果：

Eigenvalue Damping Freq. (rad/s)

-1.00e+000 + 3.00e+000i 3.16e-001 3.16e+000

4.52%

(00.9)

3.55%

n

s n

t ζωζζω??=??=<

-1.00e+000 - 3.00e+000i 3.16e-001 3.16e+000 3测试程序：

num=10;den=[1,2,10];sys=tf(num,den);step(sys);title('G(s)=10/(s^2+2s+10)单位阶跃响应曲线') 测试结果：

峰值大小：1.35 峰值时间：1.02 过渡过程时间：3.77

G(s)=10/(s 2+2s+10)单位阶跃响应曲线

Tim e (sec)

A m p l i t u d e

0.2

0.4

0.6

0.8

1

1.2

1.4

4 当1=ζ 测试程序：

num=10;den=[1,2*sqrt(10),10];sys=tf(num,den);step(sys);title('G(s)=10/(s^2+2*s*sqrt(10)+10)单位阶跃响应曲线') 测试结果：

2Tim e (sec)

A m p l i t u d e

0.10.20.30.40.50.6

0.70.80.91

当

?=2时

测试程序：

num=10;den=[1,4*sqrt(10),10];sys=tf(num,den);step(sys);title('G(s)=10

/(s^2+4*s*sqrt(10)+10)单位阶跃响应曲线')

测试结果：

2Time (sec)

A m p l i t u d e

0123

456789

0.10.20.30.40.50.6

0.70.80.91

5 当012

1w w n =

测试程序：

num=2.5;den=[1,1,2.5];sys=tf(num,den);step(sys);title('G(s)=2.5/(s^2+1+2.5)单位阶跃响应曲线')

测试结果：

G(s)=2.5/(s 2+1+2.5)单位阶跃响应曲线

Time (sec)

A m p l i t u d e

024

681012

0.2

0.4

0.6

0.8

1

1.2

1.4

当022w w n 时

测试程序：

num=40;den=[1,4,40];sys=tf(num,den);step(sys);title('G(s)=40/(s^2+4+40)单位阶跃响应曲线') 测试结果：

G(s)=40/(s 2+4+40)单位阶跃响应曲线

Time (sec)

A m p l i t u d e

00.51

1.52

2.53

0.2

0.4

0.6

0.8

1

1.2

1.4

2. 作出以下系统的阶跃响应，并与原系统响应曲线进行比较，作出相应的实验分析结果

（1）()10

210

221+++=s s s s G ，有系统零点的情况

测试程序：

G0=tf([10],[1,2,10]);G1=tf([2,10],[1,2,10]);step(G0,G1);title('现系统与原系统的系统响应曲线')

测试结果：

现系统与原系统的系统响应曲线

Time (sec)

012

3456

（2）()10

210

5.0222++++=s s s s s G ，分子、分母多项式阶数相等

测试程序：

G0=tf([10],[1,2,10]);G1=tf([1,0.5,10],[1,2,10]);step(G0,G1);title('现系统与原系统的系统响应曲线')

测试结果：

现系统与原系统的系统响应曲线

Time (sec)

（3）()10

25.0222+++=s s s

s s G ，分子多项式零次项为零

测试程序：

G0=tf([10],[1,2,10]);G1=tf([1,0.5,0],[1,2,10]);step(G0,G1);title('现系统与原系统的系统响应曲线') 测试结果：

-0.4

-0.200.20.40.60.811.21.4现系统与原系统的系统响应曲线

Time (sec)

A m p l i t u d e

（4）()10

22

2++=

s s s

s G ，原响应的微分，微分系数为1/10 测试程序：

G0=tf([10],[1,2,10]);G1=tf([1,0],[1,2,10]);step(G0,G1);title('现系统与原系统的系统响应曲线')

测试结果：

现系统与原系统的系统响应曲线

Time (sec)

A m p l i t u d e

01234567

-0.2

00.20.40.60.8

11.21.4

3. 单位阶跃响应： 25425

)()(2

++=s s s R s C 求该系统单位阶跃响应曲线,并在所得图形上加网格线和标题

测试程序：

num=25;den=[1,4,25];sys=tf(num,den);step(sys);title('G(s)=25/(s^2+4s+25)单位阶跃响应曲线'),grid on

测试结果：

Time (sec)

A m p l i t u d e

00.51

1.52

2.53

0.2

0.4

0.6

0.8

1

1.2

1.4

（三）系统动态特性分析 用Matlab 求二阶系统12012120)(2++=

s s s G 和01

.0002.001.0)(2

++=s s s G 的峰值时间p t 上升时间r t 调整时间s t 超调量%σ。

测试程序：

G0=tf([120],[1,12,120]);step(G0);title('G0=120/(s^2+12s+120)单位阶跃响应曲线)'

测试结果：

G0=120/(s 2+12s+120)单位阶跃响应曲线

Time (sec)

A m p l i t u d e

0.2

0.4

0.6

0.8

1

1.2

1.4

测试程序：

G1=tf([0.01],[1,0.002,0.01]);step(G1);title('G1=0.01/(s^2+0.002s+0.01)响应曲线') 测试结果：

G1=0.01/(s 2+0.002s+0.01)响应曲线

Time (sec)

A m p l i t u d e

10002000

3000400050006000

实验二 MATLAB 及仿真实验（控制系统的根轨迹分析）

一 实验目的

1．利用计算机完成控制系统的根轨迹作图 2．了解控制系统根轨迹图的一般规律 3．利用根轨迹图进行系统分析 二 预习要点

1. 预习什么是系统根轨迹？

2. 闭环系统根轨迹绘制规则。 三 实验方法

（一）方法：当系统中的开环增益k 从0到变化时，闭环特征方程的根在复

平面上的一组曲线为根轨迹。设系统的开环传函为：)

()

()(0s Q s N k

s G =，则系统的闭环特征方程为：0)

()

(1)(10=+=+s Q s N k

s G 根轨迹即是描述上面方程的根，随k 变化在复平面的分布。

（二）MATLAB 画根轨迹的函数常用格式：利用Matlab 绘制控制系统的根轨

迹主要用pzmap ，rlocus ，rlocfind ，sgrid 函数。

1、零极点图绘制

? [p,z]=pzmap(a,b,c,d)：返回状态空间描述系统的极点矢量和零点矢量，而不在屏幕上绘制出零极点图。

? [p,z]=pzmap(num,den)：返回传递函数描述系统的极点矢量和零点矢量，而不在屏幕上绘制出零极点图。

? pzmap(a,b,c,d)或pzmap(num,den)：不带输出参数项，则直接在s 复平面上绘制出系统对应的零极点位置，极点用×表示，零点用o 表示。 ? pzmap(p,z)：根据系统已知的零极点列向量或行向量直接在s 复平面上绘制出对应的零极点位置，极点用×表示，零点用o 表示。 2、根轨迹图绘制

? rlocus(a,b,c,d)或者rlocus(num,den)：根据SISO 开环系统的状态空间描述模型和传递函数模型，直接绘制出系统的根轨迹图。开环增益的值从零到无穷大变化。

? rlocus(a,b,c,d,k)或rlocus(num,den,k)： 通过开环增益k 范围来绘制根轨迹图。

? r=rlocus(num,den,k) 或者[r,k]=rlocus(num,den) ：不在屏幕上直接绘出系统的根轨迹图，而根据开环增益变化矢量k ，返回闭环系统特征方程1＋k*num(s)/den(s)=0的根r ，它有length(k)行，length(den)-1列，每行对应某个k 值时的所有闭环极点。或者同时返回k 与r 。

? 若给出传递函数描述系统的分子项num 为负，则利用rlocus 函数绘制的是系统的零度根轨迹。（正反馈系统或非最小相位系统） 3、rlocfind()函数

? [k,p]=rlocfind(a,b,c,d)或者[k,p]=rlocfind(num,den)

它要求在屏幕上先已经绘制好有关的根轨迹图。然后，此命令将产生一个光标以用来选择希望的闭环极点。命令执行结果：k 为对应选择点处根轨迹开环增益；p 为此点处的系统闭环特征根。

? 不带输出参数项[k,p]时，同样执行，只是此时只将k 的值返回到缺省变量ans 中。 4、sgrid()函数

? sgrid ：在现存的屏幕根轨迹或零极点图上绘制出自然振荡频率wn 、阻尼比矢量z 对应的格线。

? sgrid(‘new’)：是先清屏，再画格线。

? sgrid(z,wn)：则绘制由用户指定的阻尼比矢量z 、自然振荡频率wn 的格线。 四 实验内容 1．()()()

21++=s s s k s G g

要求：

（a ） 记录根轨迹的起点、终点与根轨迹的条数； （b ） 确定根轨迹的分离点与相应的根轨迹增益； （c ） 确定临界稳定时的根轨迹增益gL k a 测试程序：

G=tf([0,1],[1,3,2,0]); figure(1) pzmap(G); figure(2) rlocus(G);

测试结果：

Pole-Zero Map

Real Axis

I m a g i n a r y A x i s

-5-4-3-2-1012345Root Locus

Real Axis

I m a g i n a r y A x i s

b 测试程序：G=tf([0,1],[1,3,2,0]); rlocus(G);

title('根轨迹图形分离点') 测试结果：

根轨迹图形分离点

Real Axis

I m a g i n a r y A x i s

-5-4-3-2-1012

345

c 测试程序：z=[];p=[0,-1,-2];k=1; G=zpk(z,p,k); rlocus(G);

title('临界稳态根轨迹增益Kgl'); [p,k]=rlocfind(G)

测试结果：

Select a point in the graphics window

selected_point = 0.0059 + 1.3509i p =

5.5029 k =

-2.9536 -0.0232 + 1.3648i -0.0232 - 1.3648i

-5-4-3-2-1012

345临界稳态根轨迹增益Kgl

Real Axis

I m a g i n a r y A x i s

2．()()

23)(++=

s s s K s G g 要求：确定系统具有最大超调量时的根轨迹增益； 测试程序：

G0=tf([5.5 16.5],[1 2 0]);

G=feedback(G0,1,-1) step(G) 测试结果：此时Kg=5.5

Step Response

Tim e (sec)

A m p l i t u d e

00.5

1 1.5

0.2

0.4

0.6

0.8

1

1.2

1.4

3．绘制下列各系统根轨迹图。

测试程序：

G=tf([1,2,4],[1,14,65,106,24,0]);sys1=feedback(G,1,-1);rlocus(sys1);t itle('系统根轨迹图')

测试结果：

-15-10

-5

510

15

系统根轨迹图

Real Axis

I m a g i n a r y A x i s

4．绘制下列各系统根轨迹图。开环传递函数: （1）)

6.3()

2.0()()(2

++=

s s s k s H s G ； 测试程序：

G=tf([1,0.2],[1,3.6,0,0]); rlocus(G);

title('根轨迹图') 测试结果：

-1-0.8-0.6-0.4-0.200.20.4

0.60.81根轨迹图

Real Axis

I m a g i n a r y A x i s

（2）)

106.0)(5.0()()(2

+++=

s s s s k

s H s G

测试程序：

G=tf([0,1],[1,1.1,10.3,5,0]); rlocus(G);

title('根轨迹图') 测试结果：

根轨迹图

Real Axis

I m a g i n a r y A x i s

5．试绘制下面系统根轨迹图

测试程序：

G=tf([1,1],[1,3,12,-16,0]); sys1=feedback(G,1,-1); rlocus(sys1);

title('系统根轨迹图')

测试结果：

-20

-15-10

-50510

-20

-15-10-5051015

20系统根轨迹图

Real Axis

I m a g i n a r y A x i s

Matlab上机实验答案

Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -

>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +

Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50

+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =

MATLAB全部实验及答案

MATLAB全部实验及答案 实验一、MATLAB基本操作 实验内容及步骤 4、有关向量、矩阵或数组的一些运算 （1）设A=15；B=20；求C=A+B与c=a+b？ （2）设A=[1 2 3;4 5 6;7 8 9]，B=[9 8 7;6 5 4;3 2 1]；求A*B与 A.*B？ A*B就是线代里面的矩阵相乘 A.*B是对应位置的元素相乘（3）设a=10，b=20；求i=a/b=0.5与j=a\b=2？ （4）设a=[1 -2 3;4 5 -4;5 -6 7] 请设计出程序，分别找出小于0的矩阵元素及其位置（单下标、全 下标的形式），并将其单下标转换成全下标。 clear,clc a=[1 -2 3;4 5 -4;5 -6 7]; [x,y]=find(a<0); c=[]; for i=1:length(x) c(i,1)=a(x(i),y(i)); c(i,2)=x(i); c(i,3)=y(i); c(i,4)=(y(i)-1)*size(a,2)+x(i); end c

（5）在MATLAB命令行窗口运行A=[1,2;3,4]+i*[5,6;7,8]；看结果如何？如果改成运行A=[1,2;3,4]+i[5,6;7,8]，结果又如何？前面那 个是虚数矩阵，后面那个出错 （6）请写出完成下列计算的指令： a=[1 2 3;3 4 2;5 2 3]，求a^2=？,a.^2=？ a^2= 22 16 16 25 26 23 26 24 28 a.^2= 1 4 9 9 16 4 25 4 9 （7）有一段指令如下，请思考并说明运行结果及其原因 clear X=[1 2;8 9;3 6]; X( : ) 转化为列向量 （8）使用三元组方法，创建下列稀疏矩阵 2 0 8 0 0 0 0 1 0 4 0 0 6 0 0 0 方法一： clear,clc

计算方法_全主元消去法_matlab程序

%求四阶线性方程组的MA TLAB程序 clear Ab=[0.001 2 1 5 1; 3 - 4 0.1 -2 2; 2 -1 2 0.01 3; 1.1 6 2.3 9 4];%增广矩阵 num=[1 2 3 4];%未知量x的对应序号 for i=1:3 A=abs(Ab(i:4,i:4));%系数矩阵取绝对值 [r,c]=find(A==max(A(:))); r=r+i-1;%最大值对应行号 c=c+i-1;%最大值对应列号 q=Ab(r,:),Ab(r,:)=Ab(i,:),Ab(i,:)=q;%行变换 w=Ab(:,c),Ab(:,c)=Ab(:,i),Ab(:,i)=w;%列变换 n=num(i),num(i)=num(c),num(c)=n;%列变换引起未知量x次序变化for j=i:3 Ab(j+1,:)=-Ab(j+1,i)*Ab(i,:)/Ab(i,i)+Ab(j+1,:);%消去过程 end end %最后得到系数矩阵为上三角矩阵 %回代算法求解上三角形方程组 x(4)=Ab(4,5)/Ab(4,4); x(3)=(Ab(3,5)-Ab(3,4)*x(4))/Ab(3,3); x(2)=(Ab(2,5)-Ab(2,3)*x(3)-Ab(2,4)*x(4))/Ab(2,2); x(1)=(Ab(1,5)-Ab(1,2)*x(2)-Ab(1,3)*x(3)-Ab(1,4)*x(4))/Ab(1,1); for s=1:4 fprintf('未知量x%g =%g\n',num(s),x(s)) end %验证如下 %A=[0.001 2 1 5 1; 3 -4 0.1 -2 2;2 -1 2 0.01 3; 1.1 6 2.3 9 4]; %b=[1 2 3 4]'; %x=A\b; %x1= 1.0308 %x2= 0.3144 %x3= 0.6267 %x4= -0.0513

matlab实验五答案1

实验五 1、编写程序，该程序在同一窗口中绘制函数在 []0,2π之间的正弦曲线和余弦曲线，步长 为200/π，线宽为2个象素，正弦曲线设置为蓝色实线，余弦曲线颜色设置为红色虚线，两条曲线交点处，用红色星号标记并标注sin(x)=cos(x)；通过函数方式在生成的图形中添加注释，至少应包括：标题，文本注释，图例和坐标轴标注。%x=linspace(0,2*pi,1000);x=[0:pi/200:2*pi]sinx =sin(x);cosx =cos(x); k=find(abs(sinx-cosx)<1e-2);x1=x(k); plot(x,sinx,'LineWidth',2) hold on ,plot(x,cosx,'r:','LineWidth',2)hold on ,plot(x1,sin(x1),'r*')xlabel('x:(0-2\pi)'); ylabel('y:sin(x)/cos(x)'); title('正弦-余弦曲线'); text(x1+0.1,sin(x1),'sin(x)=cos(x)');legend('sin(x)','cos(x)'); 1 2 34 5 6 7 x:(0-2π) y :s i n (x )/c o s (x ) 正弦-余弦曲线 2、绘制图像：双曲抛物面：22 x y z =- ，1616x -<<，44y -<<，并对绘制的双曲抛物面尝试进行视点控制。 [X,Y]=meshgrid(-16:0.4:16,-4:0.1:4);Z =X.^2/16-Y.^2/4; subplot(1,3,1),plot3(X,Y,Z),view(0,180),title('azimuth =0,elevation =180'); subplot(1,3,2),plot3(X,Y,Z),view(-37.5,-30),title('azimuth =

MATLAB)课后实验答案

实验一 MATLAB 运算基础 1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。 (1) 012 2sin 851z e =+ (2) 21ln(2 z x =，其中2 120.45 5i x +??=? ?-?? (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022 a a e e a z a a --+= ++=--L (4) 22 42011 122123t t z t t t t t ?≤

4. 完成下列操作： (1) 求[100,999]之间能被21整除的数的个数。 (2) 建立一个字符串向量，删除其中的大写字母。解：(1) 结果： (2). 建立一个字符串向量例如： ch='ABC123d4e56Fg9';则要求结果是：

实验二 MATLAB 矩阵分析与处理 1. 设有分块矩阵33322322E R A O S ?????? =? ??? ，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22 E R RS A O S +?? =???? 。 解: M 文件如下； 5. 下面是一个线性方程组： 1 231 1 12340.951110.673450.5211145 6x x x ?? ??????????????=??? ?????????????????? ? ch = 123d4e56g9

(1) 求方程的解。 (2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。 (3) 计算系数矩阵A 的条件数并分析结论。 解： M 文件如下： 实验三 选择结构程序设计 1. 求分段函数的值。 2226035605231x x x x y x x x x x x x ?+-<≠-? =-+≤<≠≠??--? 且且及其他 用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 解：M 文件如下：

实验5 Matlab程序设计1

实验5 Matlab 程序设计1 实验目的： 1、 掌握建立和执行M 文件的方法； 2、 掌握实现选择结构的方法； 3、 掌握实现循环结构的方法。 实验内容： 1. 从键盘输入一个4位整数，按如下规则加密后输出。加密规则：每位数字都加上7，然 后用和除以10的余数取代该数字；再把第一位与第三位交换，第二位与第四位交换。 2. 求分段函数的值。 ,x x x x y x x x x x x x ?+- <≠-?=-+ ≤<≠≠??-- ?2226035605231且且及其他 用if 语句实现，分别输出x=-5,-3,0,1,2,2.5,3,5时的y 值。请输入x 的值 x=input('请输入x 的值'); if x<0&x~=-3 y=x*x+x-6; elseif 0<=x&x<5&x~=2&x~=3 y=x*x-5*x+6; else y=x*x-x-1; end y

second 请输入x的值-5 y = 14 second 请输入x的值-3 y = 11 >> second 请输入x的值0 y = 6 >> second 请输入x的值1

y = 2 >> second 请输入x的值2 y = 1 >> second 请输入x的值2.5 y = -0.2500 >> second 请输入x的值3

5 >> second 请输入x的值5 y = 19 >> 3.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E，其中90~100分为A，80~89 分为B，70~79分为C，60~69分为D，60分以下为E。 要求： （1）分别用if语句和swich语句实现。 （2）grade=input('input grade:'); （3）if grade<=100&grade>=90 （4）level='A'; （5）elseif grade<=89&grade>80 （6）level='B'; （7）elseif 70<=grade&grade<=79

MATLAB实验题答案

1、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A

matlab所有实验及答案

实验二 习题 1、 矩阵Y= ???? ? ???? ???3472123100451150425 ，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。取出子矩阵?? ? ? ??21301，并求该矩阵的维数。 解：命令为： Y=[5,2,4;0,15,1;45,100,23;21,47,3] Y(2,3) Y(10) sub2ind([4 3],2,3) [i,j]=ind2sub([4 3],10) find(Y==100) sub2ind([4 3],3,2) B=Y(2:2:4,3:-2:1) 或 B=Y([2 4],[3 1]) [m n]=size(Y) 2、 建立一个数值范围为0—100内4*5的整数随机矩阵。 求出大于50的元素的位置。 解：命令为： G=int8(100*rand(4,5)) find(G>50) 3、 已知矩阵A=[1 0 -1 ；2 4 1； -2 0 5]，B=[0 -1 0；2 1 3；1 1 2] 求2A+B 、A 2-3B 、A*B 、B*A 、A .*B ，A/B 、A\B 解：命令为： A=[1 0 -1 ；2 4 1； -2 0 5] B=[0 -1 0；2 1 3；1 1 2] E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B 4、 利用函数产生3*4阶单位矩阵和全部元素都为8的 4*4阶矩阵，并计算两者的乘积。 解：命令为： A=eye(3,4) B=8*ones(4) C=A*B 5、 创建矩阵a=????? ???????------7023021.5003.120498601 ,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b 构成矩阵d ，计算a*b 、c

MATLAB实验五 函数文件

MATLAB实验报告 学院：光电学院 班级：073-1 姓名：刘颖 学号：200713503117

实验五 函数文件 1.定义一个函数文件，求给定复数的指数、对数、正弦和余弦，并在命令文件中调用该函数文件。 程序设计： function [e ln s c]=num(x) e=exp(x) ln=log(x) s=sin(x) c=cos(x) end 运行结果： >> num(5i) e = 0.2837 - 0.9589i ln = 1.6094 + 1.5708i s = 0 +74.2032i c = 74.2099 ans = 0.2837 - 0.9589i 2.一物理系统可用下列方程组来表示: ??? ? ??? ???????= ?????? ??? ??? ???????????? ??----g g m m N N a a m m m m 2121212 111001cos 0 0sin 00cos 0 sin 0sin cos θ θθ θθθ 从键盘输入 m 1 、 m 2 和θ的值，求 N a a 121、、和 N 2 的值。其中g 取9.8，输入θ时以角度为单位。 程序设计： 函数文件in.m: function [a1,a2,N1,N2]=in(m1,m2,t) g=9.8; A=[m1*cos(t) -m1 -sin(t) 0;m1*sin(t) 0 cos(t) 0;0 m2 -sin(t) 0;0 0 -cos(t) 1]; C=[0;m1*g;0;m2*g]; B=inv(A)*C; a1=B(1); a2=B(2); N1=B(3); N2=B(4); end 调用in.m 的命令文件： >> m1=1;m2=2;t=30*pi/180; >> [a1,a2,N1,N2]=in(m1,m2,t) 运行结果： a1 = 6.5333 a2 = 1.8860 N1 = 7.5440 N2 = 26.1333 4.设 f(x)= 01 .01 1 .01 ) 3() 2(4 2 +++--x x ， 编写一个MATLAB 函数文件fx.m ，使得调用f(x)时，x 可用矩阵代入，得出的f(x)为同阶矩阵。 程序设计： 函数文件fx.m: function A=fx(x) A=1./((x-2).^2+0.1)+1./(((x-3).^4)+0.01) end 调用fx.m 的命令文件： >> A=fx([1 2;2 3;4 3]) 运行结果: A = 0.9716 10.9901 10.9901 100.9091 1.2340 100.9091 5.已知y= ) 20()30() 40(f f f + （1）当f(n)=n+10ln(n 2+5)时，求y 的值。

matlab实验四及其答案

实验四GUI 设计专业 学号姓名成绩电气工程及其自 动化201409140305 杨诚1.创建GUI 绘制方程c bx ax y ++=2图形，需要显示绘图结果的坐标系窗口，还能够输入 a ， b ， c 的值和x 取值范围（最大值和最小值）。 回调函数的编写： a=str2num(get(handles.edit1,'String'));b=str2num(get(handles.edit2,'String'));c=str2num(get(handles.edit3,'String'));xmin=str2num(get(handles.edit4,'String'));xmax=str2num(get(handles.edit5,'String'));x=xmin:0.1:xmax;y=a*x.^2+b*x+c;plot(x,y); 设计的运行界面截图：

2.设计一个GUI，完成画出y=sin(x)、y=cos(x)和y=x 3.的波形图。回调函数的编写： x=-5:0.1:5 plot(x,sin(x)) x=-5:0.1:5 plot(x,cos(x))

x=-5:0.1:5 y=x.^3 plot(x,y)

设计的运行界面截图： 3.创建一个GUI，含有下拉菜单，下拉菜单中有背景颜色选择。回调函数的编写： yanse=get(handles.popupmenu1,'value'); switch yanse case1 set(gcf,'color','r'); case2 set(gcf,'color','y'); case3 set(gcf,'color','g'); end 设计的运行界面截图：

(整理)matlab16常用计算方法.

常用计算方法 1．超越方程的求解 一超越方程为 x (2ln x – 3) -100 = 0 求超越方程的解。 [算法]方法一：用迭代算法。将方程改为 01002ln()3 x x =- 其中x 0是一个初始值，由此计算终值x 。取最大误差为e = 10-4，当| x - x 0| > e 时，就用x 的值换成x 0的值，重新进行计算；否则| x - x 0| < e 为止。 [程序]P1_1abs.m 如下。 %超越方程的迭代算法 clear %清除变量 x0=30; %初始值 xx=[]; %空向量 while 1 %无限循环 x=100/(2*log(x0)-3); %迭代运算 xx=[xx,x]; %连接结果 if length(xx)>1000,break ,end %如果项数太多则退出循环（暗示发散） if abs(x0-x)<1e-4,break ,end %当精度足够高时退出循环 x0=x; %替换初值 end %结束循环 figure %创建图形窗口 plot(xx,'.-','LineWidth',2,'MarkerSize',12)%画迭代线'.-'表示每个点用.来表示，再用线连接 grid on %加网格 fs=16; %字体大小 title('超越方程的迭代折线','fontsize',fs)%标题 xlabel('\itn','fontsize',fs) %x 标签 ylabel('\itx','fontsize',fs) %y 标签 text(length(xx),xx(end),num2str(xx(end)),'fontsize',fs)%显示结果 [图示]用下标作为自变量画迭代的折线。如P0_20_1图所示，当最大误差为10-4时，需要迭代19次才能达到精度，超越方程的解为27.539。 [算法]方法二：用求零函数和求解函数。将方程改为函数 100()2ln()3f x x x =-- MATLAB 求零函数为fzero ，fzero 函数的格式之一是 x = fzero(f,x0) 其中，f 表示求解的函数文件，x0是估计值。fzero 函数的格式之二是 x = fzero(f,[x1,x2])

(完整版)Matlab实验5选择结构程序结构

实验五、选择与循环结构 一、实验目的： 1、 掌握建立和执行M 文件的方法。 2、 掌握利用if 语句实现选择结构的方法。 3、 掌握利用switch 语句实现多分支选择结构的方法。 4、 掌握try 语句的使用。 5、 掌握利用for 语句实现循环结构的方法。 6、 掌握利用while 语句实现循环结构的方法。 7、 熟悉利用向量运算来代替循环的操作方法。 二、实验内容： 1、 列分段函数的值。 ?? ???--≠≠<≤+--≠<-+=其他且且,632,100,6530,6222x x x x x x x x x x x y 要求： (1) 用if 语句实现，分别输出x =-0.5,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。 提示：x 的值从键盘输入，可以是向量。 %homework_5_1_1.m x=input('请输入x 的值：x='); if (x<0 & x~=-3) y= x.*x + x - 6 elseif (x>=0 & x<10 & x~=2 & x~=3) y=x.*x-5.*x+6 else y=x.*x-x-6 end >> homework_5_1 请输入x 的值：x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] y = -5.2500 6.0000 -6.0000 -4.0000 -2.2500 0 14.0000 (2) 用逻辑表达式实现上述函数。 %homework_5_1_2.m x=input('请输入x 的值：x=') y=(x<0 & x~=-3).*(x.*x+x-6)... +(x>=0 & x<10 &x~=2 &x~=3).*(x.*x-5.*x+6)... +(x>=10 | x==-3 | x==3 | x==2).*(x.*x-x-6) >> homework_5_1_2 请输入x=[-0.5 -3.0 1.0 2.0 2.5 3.0 5.0] x = -0.5000 -3.0000 1.0000 2.0000 2.5000 3.0000 5.0000 y = -6.2500 6.0000 2.0000 -4.0000 -0.2500 0 6.0000

MATLAB实验题目及答案

实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0

matlab用于计算方法的源程序

1、Newdon迭代法求解非线性方程 function [x k t]=NewdonToEquation(f,df,x0,eps) %牛顿迭代法解线性方程 %[x k t]=NewdonToEquation(f,df,x0,eps) %x:近似解 %k:迭代次数 %t:运算时间 %f:原函数，定义为内联函数 ?:函数的倒数，定义为内联函数 %x0:初始值 %eps:误差限 % %应用举例： %f=inline('x^3+4*x^2-10'); ?=inline('3*x^2+8*x'); %x=NewdonToEquation(f,df,1,0.5e-6) %[x k]=NewdonToEquation(f,df,1,0.5e-6) %[x k t]=NewdonToEquation(f,df,1,0.5e-6) %函数的最后一个参数也可以不写。默认情况下，eps=0.5e-6 %[x k t]=NewdonToEquation(f,df,1) if nargin==3 eps="0".5e-6; end tic; k=0; while 1 x="x0-f"(x0)./df(x0); k="k"+1; if abs(x-x0) < eps || k >30 break; end x0=x; end t=toc; if k >= 30 disp('迭代次数太多。'); x="0"; t="0"; end

2、Newdon迭代法求解非线性方程组 function y="NewdonF"(x) %牛顿迭代法解非线性方程组的测试函数 %定义是必须定义为列向量 y(1,1)=x(1).^2-10*x(1)+x(2).^2+8; y(2,1)=x(1).*x(2).^2+x(1)-10*x(2)+8; return; function y="NewdonDF"(x) %牛顿迭代法解非线性方程组的测试函数的导数 y(1,1)=2*x(1)-10; y(1,2)=2*x(2); y(2,1)=x(2).^+1; y(2,2)=2*x(1).*x(2)-10; return; 以上两个函数仅供下面程序的测试 function [x k t]=NewdonToEquations(f,df,x0,eps) %牛顿迭代法解非线性方程组 %[x k t]=NewdonToEquations(f,df,x0,eps) %x:近似解 %k:迭代次数 %t:运算时间 %f:方程组（事先定义） ?:方程组的导数（事先定义） %x0:初始值 %eps:误差限 % %说明：由于虚参f和df的类型都是函数，使用前需要事先在当前目录下采用函数M文件定义% 另外在使用此函数求解非线性方程组时，需要在函数名前加符号“@”，如下所示 % %应用举例： %x0=[0,0];eps=0.5e-6; %x=NewdonToEquations(@NewdonF,@NewdonDF,x0,eps) %[x k]=NewdonToEquations(@NewdonF,@NewdonDF,x0,eps) %[x k t]=NewdonToEquations(@NewdonF,@NewdonDF,x0,eps) %函数的最后一个参数也可以不写。默认情况下，eps=0.5e-6 %[x k t]=NewdonToEquations(@NewdonF,@NewdonDF,x0,eps)

实验五MATLAB的基本绘图方法

实验三MATLAB的基本绘图方法 一、实验目的 1．二维平面图形的绘制 2．三维立体图形的绘制 3．隐函数作图 二、实验地点：A404 三、实验日期： 四、实验内容 （一）二维平面图形的绘制 1、Plot的使用方法介绍 plot 是绘制二维图形的最基本函数，它是针对向量或矩阵的列来绘制曲线的。也就是说，使用plot 函数之前，必须首先定义好曲线上每一点的x 及y 坐标，常用格式为：（1）plot(x) 当x 为一向量时，以x 元素的值为纵坐标，x 的序号为横坐标值绘制 曲线。当x 为一实矩阵时，则以其序号为横坐标，按列绘制每列元素值相对于其序号的曲线，当x 为m×n 矩阵时，就由n 条曲线。 （2）plot(x,y) 以x 元素为横坐标值，y 元素为纵坐标值绘制曲线。 （3）plot(x,y1,x,y2,…) 以公共的x 元素为横坐标值，以y1,y2,…元素为纵坐标值绘制多条曲线。 例1：画出一条正弦曲线和一条余弦曲线。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2) 注：在绘制曲线图形时，常常采用多种颜色或线型来区分不同的数据组，MA TLAB 软件专门提供了这方面的参数选项，我们只要在每个坐标后加上相关字符串，就可实现它们的功能。具体参见教材。 2、图形修饰 MATLAB 软件为用户提供了一些特殊的图形函数，用于修饰已经绘制好的图形。 图形修饰函数表如下： 函数含义 grid on (/off) 给当前图形标记添加（取消）网络 xlable(‘string’) 标记横坐标 ylabel(‘string’) 标记纵坐标 title(‘string’) 给图形添加标题 text(x,y,’string’) 在图形的任意位置增加说明性文本信息 gtext(‘string’) 利用鼠标添加说明性文本信息 axis([xmin xmax ymin ymax]) 设置坐标轴的最小最大值 例2、给例1的图形中加入网络和标记。 >> x=0:pi/10:2*pi; >> y1=sin(x); >> y2=cos(x); >> plot(x,y1,x,y2)

Matlab实验指导书(含答案)汇总

实验一：Matlab操作环境熟悉 一、实验目的 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/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x)

Matlab实验

MATLAB实验报告 学校：湖北文理学院 学院：物理与电子工程学院 专业：电子信息工程 学号： 2013128182 姓名：张冲 指导教师：宋立新

实验一 MATLAB环境的熟悉与基本运算 一、实验目的： 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验内容 1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明， 学习使用指令eye（其它不会用的指令，依照此方法类推） 2、学习使用clc、clear，观察command window、command history和workspace 等窗口的变化结果。 3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、 exerc2、exerc3……），学习使用MATLAB的基本运算符。 三、练习 1)help rand，然后随机生成一个2×6的数组，观察command window、 command history和workspace等窗口的变化结果。 2)学习使用clc、clear，了解其功能和作用。 3)用逻辑表达式求下列分段函数的值 4)求[100,999]之间能被21整除的数的个数。（提示：rem,sum的用法） 四、实验结果 1）

2）clc:清除命令窗口所有内容，数值不变；clear:初始化变量的值。3） 4）

实验二 MATLAB数值运算 一、实验目的 1、掌握矩阵的基本运算 2、掌握矩阵的数组运算 二、实验内容 1)输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10) 2)输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，在命令窗 口中执行下列表达式，掌握其含义： A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A 3)二维数组的创建和寻访，创建一个二维数组（4×8）A，查询数组A第2 行、第3列的元素，查询数组A第2行的所有元素，查询数组A第6列的所有 元素。 4)两种运算指令形式和实质内涵的比较。设有3个二维数组A 2×4，B 2×4 ，C 2×2 ， 写出所有由2个数组参与的合法的数组运算和矩阵指令。 5)学习使用表4列的常用函数（通过help方法） 6)学习使用表5数组操作函数。 7)生成一个3行3列的随机矩阵，并逆时针旋转90°，左右翻转，上下翻转。 8)已知a=[1 2 3]，b=[4 5 6]，求a.\b和a./ b 9)用reshape指令生成下列矩阵，并取出方框内的数组元素。 三、实验结果 1）C(i)表示C中的第i个的数值；

实验5 Matlab绘图操作实验报告

Tutorial 5 实验报告 实验名称：Matlab 绘图操作 实验目的： 1、 掌握绘制二维图形的常用函数； 2、 掌握绘制三维图形的常用函数； 3、 掌握绘制图形的辅助操作。 实验内容： 1. 设sin .cos x y x x ?? =+ ??+?? 23051，在x=0~2π区间取101点，绘制函数的曲线。 2. 已知： y x =21，cos()y x =22，y y y =?312，完成下列操作： （1） 在同一坐标系下用不同的颜色和线性绘制三条曲线； （2） 以子图形式绘制三条曲线； （3） 分别用条形图、阶梯图、杆图和填充图绘制三条曲线。 3. 已知：ln(x y x x ≤=??+>??0102 ，在x -≤≤55区间绘制函数曲线。 4. 绘制极坐标曲线sin()a b n ρθ=+，并分析参数a 、b 、n 对曲线形状的影响。 5．在xy 平面内选择区域[][],,-?-8888， 绘制函数z =的三种三维曲面图。 6. 用plot 函数绘制下面分段函数的曲线。 ,(),,x x f x x x x ?+>? ==??+

8. 在同一坐标轴中绘制下列两条曲线。 （1）.y x =-205 （2）sin()cos ,sin()sin x t t t y t t π=?≤≤? =?303 实验结果： 1. 2. （1）

（2）

（3）

计算方法上机实验报告-MATLAB

《计算方法》实验报告 指导教师： 学院： 班级： 团队成员：

一、题目 例2.7应用Newton 迭代法求方程210x x --=在1x =附近的数值解 k x ，并使其满足8110k k x x ---< 原理： 在方程()0f x =解的隔离区间[],a b 上选取合适的迭代初值0x ，过曲线()y f x =的点()() 00x f x ，引切线 ()()()1000:'l y f x f x x x =+- 其与x 轴相交于点：()() 0100 'f x x x f x =-，进一步,过曲线()y f x =的 点()()11x f x , 引切线 ()()()2111: 'l y f x f x x x =+- 其与x 轴相交于点：() () 1211 'f x x x f x =- 如此循环往复，可得一列逼近方程()0f x =精确解*x 的点 01k x x x ,,,,,其一般表达式为： ()() 111 'k k k k f x x x f x ---=- 该公式所表述的求解方法称为Newton 迭代法或切线法。

程序： function y=f(x)%定义原函数 y=x^3-x-1; end function y1=f1(x0)%求导函数在x0点的值 syms x; t=diff(f(x),x); y1=subs(t,x,x0); end function newton_iteration(x0,tol)%输入初始迭代点x0及精度tol x1=x0-f(x0)/f1(x0);k=1;%调用f函数和f1函数 while abs(x1-x0)>=tol x0=x1;x1=x0-f(x0)/f1(x0);k=k+1; end fprintf('满足精度要求的数值为x(%d)=%1.16g\n',k,x1); fprintf('迭代次数为k=%d\n',k); end 结果：