matlab程序求解含有边界条件的线性方程组

《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