wrmt.net
当前位置:首页 >> mAtlAB解线性方程组 >>

mAtlAB解线性方程组

方法有很多的,说说高斯列主元消去法解一般线性方程组的做法,以下是liezy.m文件,文件名不要修改就要用这个 function[RA,RB,n,x]=liezy(A,b); B=[A b];n=length(b); RA=rank(A);RB=rank(B); zhicha=RA-RB; if zhicha>0, disp('因为RA~=RB,所以...

解线性方程组属于MATLAB最基本的应用之一。 先把A和B的写出来,然后用后面三种方法求x都可以。 A = [ 5 7 6 5 1 7 10 8 7 2 6 8 10 9 3 5 7 9 10 4 1 2 3 4 5 ]; B = [ 24 96 34 136 36 144 34 140 15 60 ]; x = A\B x = inv(A)*B x = A^-1*B

这个问题不好直接回答,因为N=30000和稀疏性其实基本上不携带任何有价值的信息。 根据需求应该先搞清楚几件事,然后才能选择算法 1.是否具有对称性,如果有的话是否还有正定性 2.需要解一次还是解很多次(相同的A不同的b),或是矩阵还会变(不...

一。用matlab 中的solve函数 >>syms x y; %定义两个符号变量; >>[x ,y]=solve('y=2*x+3','y=3*x-7');%定义一个 2x1 的数组,存放x,y >>x >>x=10.0000 >>y >>y=23.0000 二。用matlab 中的反向斜线运算符(backward slash) 分析: 方程组可化为...

线性方程组,这个很简单。 假设A是系数矩阵,b是右侧向量,x是未知数 那么x=inv(A)*b 比如 1*x(1)+2*x(2)+3*x(3)=4 4*x(1)+1*x(2)+1*x(3)=6 6*x(1)+3*x(3)=8 %那么系数矩阵:A=[1,2,3;4,1,1;6,0,3];%右侧向量:b=[4;6;8];%可以求出x=inv(A)*b x =...

给你一个例子,来说明如何用Matlab求线性方程组的通解。 >> a=[1 -1 1 -1;-1 1 1 -1;2 -2 -1 1]; %线性方程组的系数矩阵 >> b=[1;1;-1]; % 常列向量 >> [rank(a) rank([a,b])] ans = 2 2 %秩相等且小于4,说明有无穷多解 >> rref([a,b]) %简化行...

在求解线性方程组时,会遇到以下几种情形:定解方程组、不定方程组、超定方程组、奇异方程组。作为示例,首先以定解线性方程组为例: 在分析如上方程组时,需要知道,方程中有3个未知数,而方程也有3个,所以可以求出(x,y,z)值,转化为矩阵即为...

%function [RA,RB,n,X]=gaus(A,b) A=[1 2 3; 1 4 9; 1 8 27]; b=[ 1 2 3 ]; B=[A b]; n=length(b); RA=rank(A); %A的秩 RB=rank(B);%B的秩 zhica=RB-RA; if zhica>0 disp('请注意:因为RA~=RB,所以此方程组无解.'); return; end if RA==RB %方...

设A为方程组的系数矩阵,B为方程组的常数项值,X为方程组的解。则 X=pinv(A)*B,或X=A\B。 例如:用matlab求解下列方程组 >>A=[1 2;3 -2;1 -1];B=[1;4;2]; >> pinv(A)*B ans = 1.28378378378378 -0.175675675675675

先写m文件 function [x,y]=line_solution(A,b) [m,n]=size(A); y=[]; if norm(b)>0 if rank(A)==rank([A,b]) if rank(A)==n disp('方程有唯一解x'); x=A\b; else disp('方程有无穷多解,特解为x,其齐次方程组的基础解系为y'); x=A\b; y=null(A,'...

网站首页 | 网站地图
All rights reserved Powered by www.wrmt.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com