没有合适的资源？快使用搜索试试~ 我知道了~

文库首页行业研究行业报告数值分析Matlab作业.docx

数值分析Matlab作业.docx

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

0 下载量 201 浏览量 2023-02-20 18:40:02 上传评论收藏 600KB DOCX 举报

温馨提示

试读

18页

资源推荐

资源详情

资源评论

数值分析编程作业

数值分析

2012 年 12 月

第二章

14.考虑梯形电阻电路的设计，电路如下：

电路中的各个电流{i1,i2,…,i8}须满足下列线性方程组：

2i  2i

 V / R

 0

2i  5i  2i

 0

2i  5i  2i

 0

2i  5i  2i

 0

2i  5i  2i

 0

2i  5i  2i  0

2i  5i  0

这是一个三对角方程组。设 V=220V，R=27 ，运用追赶法，求各段电路的电流量。



Matlab 程序如下：

function chase ()

%追赶法求梯形电路中各段的电流量

a=input('请输入下主对角线向量a=');

b=input('请输入主对角线向量b=');

c=input('请输入上主对角线向量c=');

d=input('请输入右端向量d=');

n=input('请输入系数矩阵维数n=');

u(1)=b(1);

for i=2:n

l(i)=a(i)/u(i-1);

u(i)=b(i)-c(i-1)*l(i);

end

y(1)=d(1);

for i=2:n

y(i)=d(i)-l(i)*y(i-1);

end

x(n)=y(n)/u(n);

i=n-1;

while i>0

x(i)=(y(i)-c(i)*x(i+1))/u(i);

i=i-1;

end

输入如下:

2（共 18 页）

数值分析

>> chase

请输入下主对角线向量a=[0,-2,-2,-2,-2,-2,-2,-2];

请输入主对角线向量b=[2,5,5,5,5,5,5,5];

请输入上主对角线向量c=[-2,-2,-2,-2,-2,-2,-2,0];

请输入方程组右端向量d=[220/27,0,0,0,0,0,0,0];

请输入系数矩阵阶数n=8

运行结果如下：

x = 8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194

0.0477

第三章

14.试分别用(1)Jacobi 迭代法；(2)Gauss-Seidel 迭代法解线性方程组



10 1 2 3 4

   





   

27



1 9 1 2 3

   



 2 1 7 3 5     14 



3 2 3 12 1

 



   

17





   

 x  



4 3 5 1 15





(0)

 (0,0,0,0,0)

迭代初始向量

。

（1）雅可比迭代法程序如下：

function jacobi()

%Jacobi迭代法

a=input('请输入系数矩阵a=');

b=input('请输入右端向量b=');

x0=input('请输入初始向量x0=');

n=input('请输入系数矩阵阶数n=');

er=input('请输入允许误差er=');

N=input('请输入最大迭代次数N=');

for i=1:n

for j=1:n

if i==j

d(i,j)=a(i,j);

else

d(i,j)=0;

end

m=eye(5)-d\a;

g=d\b;

%迭代矩阵

x=m*x0+g;

k=1;

while k<=N

for i=1:5

%进行迭代

if max(abs(x(i)-x0(i))) >er

x=m*x+g;

k=k+1;

3（共 18 页）

数值分析

else

return

end

continue

end

程序执行如下：

>>jacobi

请输入系数矩阵a=[10 1 2 3 4;1 9 -1 2 -3;2 -1 7 3 -5;3 2 3 12 -1;4 -3 -5 -1 15]

请输入右端向量b=[12 -27 14 -17 12]'

请输入初始向量x0=[0 0 0 0 0]'

请输入系数矩阵阶数n=5

请输入允许误差er=1.0e-6

请输入最大容许迭代次数N=60

x =

1.0000

-2.0000

3.0000

-2.0000

1.0000

（2）高斯－赛德尔迭代法程序如下：

function gs_sdl()

%gauss-seiddel迭代法

a=input('请输入系数矩阵a=');

b=input('请输入右端向量b=');

x0=input('请输入初始向量x0=');

n=input('请输入系数矩阵阶数n=');

er=input('请输入允许误差er=');

N=input('请输入最大迭代次数N=');

for i=1:n

for j=1:n

if i<=j

l(i,j)=0;

else

l(i,j)=-a(i,j);

end

for i=1:n

for j=1:n

if i<j

u(i,j)=-a(i,j);

else

4（共 18 页）

剩余17页未读，继续阅读

评论收藏

内容反馈

1.该资源内容由用户上传，如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款（资源遇到问题，请及时私信上传者）

版权申诉

资源评论

资源反馈

评论星级较低，若资源使用遇到问题可联系上传者，3个工作日内问题未解决可申请退款~

xxpr_ybgg

粉丝: 6520
资源: 3万+

上传资源快速赚钱

我的内容管理展开

我的资源快来上传第一个资源

我的收益

登录查看自己的收益

我的积分登录查看自己的积分

我的C币登录后查看C币余额

我的收藏

我的下载

下载帮助

前往需求广场，查看用户热搜

数值分析Matlab作业.docx

数值分析matlab代码.docx

工程数值计算matlab实验报告——辛普森数值微分数值积分.docx

数值分析第四章外推法计算数值微分MATLAB计算实验报告.docx.docx

(完整word版)数值分析的MATLAB程序.docx

数值分析matlab程序.docx

数值分析matlab实验.docx

数值分析作业MATLAB.docx

数值分析matlab完整版实验报告.docx

电子科大数值分析作业MATLAB.docx

常微分方程的数值解法及其matlab实现.docx

基于matlab的数值分析(矩阵的LU分解)实验报告.docx

matlab数值分析处理线性拟合.docx

数值分析的MATLAB程序.docx

【老生谈算法】数值分析-MATLAB相关算法.docx

数值分析原子弹爆炸作业.docx

数值计算大作业代码.docx各种插值公式matlab

数值分析数值积分实验报告.docx

Cobalt Strike下载

北京邮电大学计算机考研复试笔试资料

计算机系统-笔记-HUN2021级

cs1.6老版本供下载

合成孔径雷达的经典成像算法cs(matlab)仿真代码（吐血整理，内容全，注释全）

港大CS（MSC）面试整理

合成孔径雷达RD CS OmegaK算法点目标仿真.rar

计算机科学导论原书第二版答案.zip

Cobalt-Strike-4.5

cobaltstrike4.3.zip

在dataGridView的列中出现日历选择控件的类型

541118 深入理解计算机系统（原书第三版）课程ppt

最新资源