matlab中复合中点式程序,《现代数值计算》Matlab程序整理（23页）-原创力文档

第2 章线性方程组的直接解法

1.1 写一个求解线性代数方程x b 的列主元素高斯消去法程序，在程序入口输入n 的

值，其中ij ，bj ln j ，1 i , j  n

i  j 1

解：

function x=gauss(n)

for i=1:n

for j=1:n

A(i,j)=1/(i+j-1);

end

b(i,1)=log(i);

end

%Gauss Elimination

for i=1:n-1

[r,p]=max(abs(A(1:n,i)));

p=p+i-1;

A([p,i],:)=A([i,p],:);

for j=i+1:n

r -A(j,i)/A(i,i);

A(j,:)=A(j,:)+r*A(i,:);

b(j)= b(j)+r*b(i);

end

x=zeros(n,1);

for k=n:-1:1

x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k);

end

1.2 写一个求解线性代数方程A=[1,13,-2,-34;2,6,-7,-10;-10,-1,5,9;-3,-5,0,15];b=[13,-22,14,-36]';

的列主元素高斯消去法程序

解：

function x=gaussel(A,b)

n=length(b);

for i=1:n-1

[r,p]=max(abs(A(i:n,i))); %r 为最大值，p 为所在行

p=p+i-1;

if p>i

t=A(i,:);A(i,:)=A(p,:);A(p,:)=t;

t=b(i);b(i)=b(p);b(p)=t;

end

for j=i+1:n

r -A(j,i)/A(i,i);

A(j,:)=A(j,:)+r*A(i,:);

b(j)= b(j)+r*b(i);

end

x=zeros(n,1);

x(n)=b(n)/A(n,n);

for k=(n-1):-1:1

x(k)=(b(k)-A(k,k+1:n)*x(k+1:n))/A(k,k);

end

>> A=[1,13,-2,-34;2,6,-7,-10;-10,-1,5,9;-3,-5,0,15];

>> b=[13,-22,14,-36]';

>> x=gaussel(A,b)

14.0669

29.9869

17.9950

10.4090

2. LU 分解

function [L,U]=myLU(A) %实现对矩阵A 的LU 分解，L 为下三角矩阵

[n,n]=size(A);

L=zeros(n,n);

U=zeros(n,n);

for i=1:n

L(i,i)=1;

end

for k=1:n

for j=k:n

U(k,j)=A(k,j)-sum(L(k,1:k-1)*U(1:k-1,j));

end

for i=k+1:n

L(i,k)=(A(i,k)-sum(L(i,1:k-1)*U(1:k-1,k)))/U(k,k);

end

matlab中复合中点式程序,《现代数值计算》Matlab程序整理（23页）-原创力文档相关推荐

c语言程序中cost的作用,C语言考试题基础版（21页）-原创力文档
if (x>y)z=x; s=z*z;elsez=y;s=1/(z*z); if (x>y) z=x; s=z*z; else z=y;s=1/(z*z); 7. B. C. PAd, P ...
每个c语言程序写完后都要按照,c语言基础学习小结(习题总结)（5页）-原创力文档...
书山有路勤为径学海无涯苦作舟一.思考题. 1.你如何向别人解释清楚什么是编程.什么是计算机语言? 2.什么是C语言? 二.解答题. 1.用C语言编写程序:求任意两个整数的和.如果是小数的话,要求输 ...
linux编译错误程序中有游离的,操作系统实验报告附思考题（24页）-原创力文档...
课程设计(综合实验)报告 ( 2015 -- 2016 年度第 1 学期) 名称: 操作系统综合实验题目: oslab 综合实验院系: 计算机系班级: 学号: 学生姓名: 指导教师: ...
【MATLAB】基于油猴脚本和MATLAB下载原创力文档
前言最近在油猴找到一个脚本,挺有意思,它能够下载那些允许预览的文档,这个到是很常见的一种脚本,其原理大概就是自动截图得到图片,然后给你合成PDF文档.这也就是为什么只能下载那些能够预览的文档. ...
kiel实现c语言编程,单片机实验1(Kiel-软件的使用和C程序调试)教学文案（18页）-原创力文档...
单片机原理与接口技术实验实验一:Kiel 软件的使用和C程序调试一.实验目的: 1.通过实验学会用C语言进行程序的设计. 2.学会应用keil-uvision集成调试软件对C语言程序进行调试.掌握C ...
买三种文具编程C语言,C语言程序设计实验的个人物品管理系统程序（33页）-原创力文档...
实用标准文案实用标准文案精彩文档精彩文档 <程序设计与算法综合实习I> 个人物品管理程序2015.05.16题目: 班级: 学生学号: 学生姓名: 指导老师: 提交时间: 成绩 ...
利用计算机软件模拟光栅衍射实验,光栅衍射实验地MATLAB仿真（29页）-原创力文档...
实用标准文案届别 2012 届．．学号 200814060106 毕业设计光栅衍射实验的 MATLAB仿真姓名吴帅系别. 专业物理与电子信息工程系应用物理专业导师 ...
计算机在化工设计中的应用论文,计算机软件技术在化工设计的应用-软件技术论文-计算机论文（9页）-原创力文档...
计算机软件技术在化工设计的应用摘要:随着我国社会经济的不断发展,我国的工业生产的效率和技术水平都得到了极大地提升,尤其体现在化工行业中,通过对先进的计算机软件技术的运用,已经基本实现了对化工生产全过 ...
计算机硬件系统中指挥控制计算机工作,2017年广东省中等专业学校教师招聘考试试卷-计算机（6页）-原创力文档...
2017年全国中等专业学校教师招聘考试试卷计算机注意事项: 1.本试卷分为两部分,第一部分为选择题,第二部分为非选择题. 2.应考者必须按试题顺序在答题卡(纸)指定位置上作答,答在试卷上无效. 3 ...

matlab中复合中点式程序,《现代数值计算》Matlab程序整理（23页）-原创力文档

matlab中复合中点式程序,《现代数值计算》Matlab程序整理（23页）-原创力文档相关推荐

最新文章

热门文章