欢迎来到天天文库
浏览记录
ID:14877786
大小:40.00 KB
页数:3页
时间:2018-07-30
《matlab程序求解含有边界条件的线性方程组》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、编写程序,求解含有边界条件(本质)的线性方程组一、编写程序思路:拟编写一MATLAB程序函数,来处理边界条件(),将原线性方程组转化为形如:(边界条件为),然后调用作业(一)中解线性方程组函数x=gauss(A,b),求解该线性方程组。一、变量说明:解含边界条件()线性方程组Ax=b,其中:A-线性方程组系数矩阵;b-列向量n-系数矩阵行数;m-系数矩阵列数;x-未知解向量i-系数矩阵的行变量;k-解向量x的下标变量;三、基于MATLAB软件平台编写程序,如下:functiontr=boundary(A,b,k,a)%定义边
2、界条件处理函数[n,m]=size(A);%获得系数矩阵A的行、列数A(k,:)=0;%系数矩阵第k行化为0fori=1:n%使向量b转化为[b(i)-A(i,k)*a]向量b(i)=b(i)-A(i,k)*a; endA(:,k)=0;%系数矩阵第k列化为0A(k,k)=1;b(k)=a; %使所对应系数矩阵中主对角元素化为1disp(A); %显示边界条件经处理后的A和bdisp(b);x=gauss(A,b)%调用作业一中高斯消元函数解方程组四、程序应用例
3、如:解线性方程组: 其中,边界条件为>>A=[2,-1,5,1,1;4,2,4,8,-1;1,2,0,7,3;3,1,3,3,2;2,4,-3,0,4];>>b=[9;10;12;20;5];>>tr=boundary(A,b,2,3)输出结果为: x=6.52493.00000.34811.9613-4.7514
此文档下载收益归作者所有