www.ypnh.net > 如何用mAtlAB求解线性方程组

如何用mAtlAB求解线性方程组

一。用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) 分析: 方程组可化为...

方法有很多的,说说高斯列主元消去法解一般线性方程组的做法,以下是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=[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]) %简化行...

给你举个例子,已知两个方程表达式,可以利用solve函数求解 [x,y] = solve('-0.43*x -0.67*y -90=0','-0.32*x-0.38*y-280=0')

function [A,det]=GaJo_inv(A) % Gauss-Jordan 列主元消去法求方阵逆 P178 % [A,det]=GaJo_inv(A) % A 要求逆矩阵 % det 按需求返回 A 的行列式 % A 返回逆矩阵,放在 A 中 n=size(A); if n(1)~=n(2) error('不是方阵!'); end n=n(1);det=1; flag...

你的问题提的范围太大了。有具体的线性方程组吗? 一般求解线性方程组的最常用的方法有,solve函数,x=A/b(A是线性方程组的系数,b是等式值),x=inv(A)*b(inv(A)是A的逆矩阵)等等。

先写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,'...

x=fsolve(@fun,x0)求解fun(x)=0的解,x0是初值,fun是函数,x就是解 因为fsolve使用迭代法求解方程的,所以总要有个迭代的初值吧,这个初值就是你给的x0。 比如解方程组 x(1).^2+x(2).^2=1 x(1)=2*x(2) 可以写成 x=fsolve(f,[1 1]) 这里[1 1]就...

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

直接用矩阵算吧?给你一个 你试试看看是不是答案。 A=[20,10,10,15; 5,5,10,15; 5,15,5,10; 8,10,10,20,]; % A是你的非其次方程的系数矩阵 B=[70,35,35,50]; C=A/B % C=A/B 的意思是A除B实际上矩阵是不能除的,这是matlab的一个简单的表达,其实...

网站地图

All rights reserved Powered by www.ypnh.net

copyright ©right 2010-2021。
www.ypnh.net内容来自网络,如有侵犯请联系客服。zhit325@qq.com