数值分析作业

发布 2022-08-27 20:16:28 阅读 4584

一﹑求齐次方程组的解。

求解程序:> a=[1 2 2 1;2 1 -2 -2;1 -1 -4 -3];r=rank(a)r =

> c=null(a)c =

> c'*cans =

> cl=null(sym(a))cl =

于是解向量为[x y z w] =k+k=k [2 ﹣2 1 0] +k [﹣0 1],即x=k+2k,y=﹣2k-k,z=k,w=k,其中k,k为任意常数,为方程组的一个基础解系。

二、求方阵a3=的行列式、特征多项式、特征值和特征向量。

求解程序:> a3=[-2 1 1;0 2 0;-4 1 3];det(a3)ans =

> p=poly(a3)p =

> p1=poly2str(p,'y')

p1 =y^3 - 3 y^2 + 4

> roots(p)

ans =2.0000 + 0.0000i

2.0000 - 0.0000i

> eig(a3)ans =

> [x r]=eig(a3)x =

r => [x r]=eig(sym(a3))x =

r =三、用quadl 指令求积分,结果取14位定点小数。

求解程序。打开编辑调试窗用该函数的积分,可编成m﹣文件。

function y=li8_12(x)

y=sin(x.^2);

将上述文件以“li8_12”为名存盘,退出编辑调试窗。

在指令窗中键入:

> format long,s=quadl('li8_12',exp(-5),5)

回车得出。s =

这样可以得出理想的数值积分结果。

四、求解+2xy(x) =1,y(﹣2.5) =0

求解程序:1)现在编辑调试窗口中编出下列m﹣函数文件:

function y1=li9_3(x,y)

y1=1-2*x*y

以li9_3存盘,退回指令窗口中。

2)键入:> [x,y]=ode23(@ li9_3,[-2.5 3],[0])x =

y =1.0e+002 *

3)若键入。

ode23(@ li9_3,[-2.5 3],[0]),grid

回车得微分方程的图示解。

有上图可看出使用ode23指令时写于不写输出变量的差异。常微分方程的初值问题的数值解,就是求出界函数的列表法或图示法表示,输出数值是函数的列表法表示,输出曲线是函数的图示法表示。

数值分析作业

12.求在 0,1 上的一次最佳平方逼近多项式与二次最佳平方逼近多项式。函数 function s zjpfbj n,a,b 创建一个函数,里面填入次数,和区间范围。base inline x j 1 x j 定义多项式。quan inline 1 x 权函数。a zeros n 1 y zeros...

数值分析作业

实验2.2算法设计与比较。实验目的 编制不同方法的matlab程序,这些方法的计算效果和特点。问题提出 非线性方程的数值解法很多,不同的方法效果如何,要靠计算的实践来分析 比较。实验内容 考虑下列算法 1 牛顿法 2 弦割法 3 抛物线法。分别编写它们的matlab程序。牛顿法程序 function...

数值分析作业

数值分析上机作业。zhh mushan 班级 土木工程一班。姓名 赵惠惠。班级 2010510359 二零一一年六月 舍入误差与有效数。一 问题说明。1 设,其精确值为。1 编制按从小到大的顺序计算通用程序。2 编制按从小到大的顺序计算通用程序。3 按两种顺序分别计算并指出有效位数。编制程序用单精度...