matlab 实验报告

MATLAB 实验报告

姓名：

专业：

班级：

学号：

2009.12.17

一、 实验目的：

熟悉MATLAB 简单的使用办法；学会用MATLAB 解线性方程组。

二、 实验要求：

1、求解线性方程组⎪⎩⎪

⎨⎧=++=++=++14

1022131021210321

321321x x x x x x x x x ；

2、在MA TLAB 命令窗口编写求解线性方程组程序。

三、实验程序

（1）利用直接解法求解

>> A=[10,1,1;2,10,1;2,2,10];

>> b=[12,13,14]'; >>x=A\b

（2）利用LU 求解

>> A=[10,1,1;2,10,1;2,2,10]; >> b=[12,13,14]'; >> [L,U]=lu(A); >> x=U\(L\b)

(3)用QR 分解求解

>> A=[10,1,1;2,10,1;2,2,10]; >> b=[12,13,14]'; >> [Q,R]=qr(A); >> x=R\(Q\b)

四、实验结果

x =

1.0000 1.0000 1.0000

一、实验目的：

熟悉MATLAB简单的使用办法；学会用MATLAB解常微分方程；并画出求解的图形。

二、实验要求：

1、求解常微分方程⎪

⎩

⎪

⎨

⎧

=

-

=

1

)0(

2

'

y

y

x

y

y

；

2、在MATLAB命令窗口编写求解常微分方程程序；

3、求出x在区间

[]20,0上各节点上的数值解并画图。

三、实验程序

(1)建立函数文件f.m。

function f=f(x,y)

f=y-2*x/y

(2) 在命令窗口求解微分方程，结果如下

>> x0=0;xf=20;

>> y0=1;

>> [x,y]=ode23('f',[x0,xf],y0)

四、实验结果

一、 实验目的：

熟悉MATLAB 简单的使用办法；学会用matlab 求解简单极限。

二、实验要求：

1、求极限

x

x x ⎪⎭⎫ ⎝⎛+∞

→11lim

； 2、求极限

x

x x ⎪⎭

⎫ ⎝⎛++→11lim

。 三、实验程序

（1）程序一：

syms x

limit(((1+1/x)^x),x,inf)

（2）程序：

>> syms x

>> f=(1+1/x)^x;

>> limit(f,x,0,'right')

四、实验结果

（1）结果

ans =

exp(1)

（2）结果：

ans =

1

一、实验目的：

熟悉MATLAB简单的使用办法；学会用MATLAB对字符串进行简单的操作。

二、实验要求：

'hello；

1、输入字符串变量a为'

2、将a的每个字符向后移4位，例如''h变为''l；

3、最后逆序排放赋给变量b。

三、实验程序及结果：

实验程序：

（1）>> ch='hello';

>> k=find(ch>='a'&ch<='z');

>> ch(k)=ch(k)+('l'-'h')

实验结果：

ch =

lipps

（2）

>> b=ch(end:-1:1)

实验结果：

b =

sppil

一、 实验目的：

熟悉MATLAB 简单的使用办法；学会用matlab 绘制一些函数图形。

二、 实验要求：

1、在图形窗口中分别绘制2101.0t y =，)2sin(2t y e

t

-=两条曲线，其中[]10,0∈t ；

2、绘制2y 的最大值水平线。

三、 实验程序：

1) 直接输入函数解析式

fplot(‘[0.01*t^2,exp(-t)*sin(2*t)]’,[0 10])

2) 编写M 文件

i. 建立函数文件myfun.m

f unction y=myfu n(t) y(:,1)= 0.01*t^2;

y(:,2)= exp(-t)*sin(2*t);

ii. 命令

fplot （‘myfun ‘, [0 10]）

3）x=0:0.01:10;

>> y1=0.01*x.*x;y2=exp(-x).*sin(2*x); >> max(y2);y=max(y2);plot(x,y1,x,y2,x,y);

012345678910

一、 实验目的：

熟悉MATLAB 简单的使用办法；学会用matlab 进行多项式的一些简单运算。

二、 实验要求：

对与多项式12345)(234++++=x x x x x a 和13)(2+=x x b ，进行如下操作：1）计算：)()()(x b x a x c =；2）计算：)(x c 的根；3）当2=x 时，计算)(x c 的值；4）将)

()(x a x b 进行部分分式展开。

三、 实验程序及结果：

（1）>> A=[5 4 3 2 1];

>>B=[3 0 1]; >>ab=conv(A,B) ab =

15 12 14 10 6 2 1

>> poly2sym(ab)

ans =

15*x^6+12*x^5+14*x^4+10*x^3+6*x^2+2*x+1

（2）计算c(x)的根

>> p=[15 12 14 10 6 2 1]; >> roots(p) ans =

-0.5378 + 0.3583i -0.5378 - 0.3583i 0.1378 + 0.6782i

0.1378 - 0.6782i -0.0000 + 0.5774i -0.0000 - 0.5774i

（3） >> p=[15 12 14 10 6 2 1];

>>poly2sym(p) ans =

15*x^6+12*x^5+14*x^4+10*x^3+6*x^2+2*x+1

>> polyval(p,2)

ans =

1677

（4） >> A=[5 4 3 2 1];

>> B=[3 0 1];

>> [r,p]=residue(B,A)

r =