matlab矩阵齐次方程求解 matlab中如何求解齐次线性方程组（代数矩阵）的非零解
\u5982\u4f55\u7528matlab\u6c42\u89e3\u9f50\u6b21\u7ebf\u6027\u65b9\u7a0b\u7ec4\u65b9\u6cd5\u6709\u5f88\u591a\u7684\uff0c\u8bf4\u8bf4\u9ad8\u65af\u5217\u4e3b\u5143\u6d88\u53bb\u6cd5\u89e3\u4e00\u822c\u7ebf\u6027\u65b9\u7a0b\u7ec4\u7684\u505a\u6cd5\uff0c\u4ee5\u4e0b\u662fliezy.m\u6587\u4ef6\uff0c\u6587\u4ef6\u540d\u4e0d\u8981\u4fee\u6539\u5c31\u8981\u7528\u8fd9\u4e2a
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('\u56e0\u4e3aRA~=RB\uff0c\u6240\u4ee5\u6b21\u65b9\u7a0b\u65e0\u89e3')
return
end
if RA==RB
if RA==n
disp('\u65b9\u7a0b\u6709\u552f\u4e00\u4e00\u7ec4\u89e3')
x=zeros(n,1); c=zeros(1,n+1);
for p=1:n-1
[Y,j]max(abs(B(p:n,p))),C=B(p,:);
B(p,:)=B(j+p-1,:);
\u7136\u540e\u518d\u547d\u4ee4\u7a97\u53e3\u8f93\u5165\uff1a
\u7cfb\u6570\u77e9\u9635A=[......,.......,.......,.......]
b=[0;0;0;0.........]
[RA,RB,n,x]=liezy(A,b)
\u7136\u540e\u5c31\u6709\u7ed3\u679c\u4e86
MATLAB\u4e2d\u7528AX=0\u7684\u65b9\u6cd5\u6c42X\u53ea\u80fd\u6c42\u51fa0\u89e3
\u8bf4\u660e\u7cfb\u6570\u77e9\u9635\u7684\u884c\u5217\u5f0f\uff0c\u88abMatlab\u8ba1\u7b97\u8ba4\u4e3a\u4e0d\u4e3a0\uff0c
\u9700\u8981\u8865\u5145\u6761\u4ef6
参考资料：矩阵的逆 矩阵的运算 AB=A+2B 移向AB-2B=A 提取B B(A-2)=A计算就好了
>> A=[4.46,-0.54,-0.54;-1.7,3.3,-1.7;-3,-3,2 ]
>> b=[0;0;0]
>> A\b
绛旓細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锛屾墍浠ユ鏂圭▼鏃犺В')return end if RA==RB if RA==n disp('鏂圭▼鏈夊敮涓涓缁勮В')x=zeros(n,1); c=zeros(1,n+1);for p=1:n-1 [Y,j]max(...
绛旓細Matlab璇█鏍煎紡锛歑=inv(A)*b %鍙嗙煩闃 X=A\b X=sym(A)\sym(b)榻愭绾挎ф柟绋嬬粍鐨勯氳В 榻愭绾挎鏂圭▼鐭╅樀褰㈠紡锛欰X=0 Matlab璇█鏍煎紡锛歓=null(A,'r')MATLAB鍒╃敤鐭╅樀姹傜嚎鎬ф柟绋嬬粍 MATLAB鍒╃敤鐭╅樀姹傜嚎鎬ф柟绋嬬粍 闈為綈娆＄嚎鎬ф柟绋嬬粍 闈為綈娆＄嚎鎬ф柟绋嬬粍鐨勪竴鑸舰寮忥細AX=b;瑙ｆ柟绋缁勫涓嬶細MATLAB鍒╃敤鐭╅樀姹傜嚎...
绛旓細鍙互鎶0.1绉诲埌绛夊彿鐨勫彸杈癸紝杩欐牱A'=[4.46,-0.54;-1.7,3.3;-3,-3],B'=[0.054;0.17;-0.2],鐒跺悗鍐嶅仛闄ゆ硶銆侫\B灏卞彲浠ヤ簡銆
绛旓細瑙ｆ硶涓锛欰=[1 1 ;2,1锛3,1锛4,1,锛5,1];B=[9.60,8.85,8.05,7.50,7.15];format rat x1=A\B %姹傚緱闈榻愭鏂圭▼缁凙x=B鐨勪竴涓壒瑙1 Y=null(A,鈥檙鈥) %姹傚緱榻愭鏂圭▼缁凙x=0 鐨勫熀纭瑙ｇ郴Y 瑙ｆ硶浜:formatrat A=[1 1 ;2,1锛3,1锛4,1,锛5,1];B=[9.60,8.85...
绛旓細缁欏嚭鐨勫井鍒嗘柟绋嬶紝鍙互鐢╫de鍑芥暟姹傝В锛屽緱鍒板叾鏁板艰В銆倄=[1,4];y0=1;ode_fun=@(x,y)-(x^2*y*sin(x)+1)/x;[x,y]=ode45(ode_fun,x,y0);plot(x,y);xlabel('x');ylabel('y');grid on 杩愯缁撴灉鐨勫浘褰㈠涓
绛旓細5銆佸浜庨潪澶氶」寮鏂圭▼锛屽彧鑳戒娇鐢╯olve姹傝В銆備緥濡傛眰瑙xp(-x)-x^2+3=0鐨勮В濡備笅銆傞氳繃榛勮壊鐨勮鍛婂彲浠ョ湅鍑猴紝杩欐牱鐨勬柟绋嬫病鏈夎В鏋愯В锛屼笌鎴戜滑宸茬煡鐨勭煡璇嗙浉鍚屻6銆佹渶鍚庝竴绫绘柟绋嬶紝鏄竴鍏冧竴娆℃柟绋嬬粍銆傝繖鏄matlab鏈鎿呴暱鐨勮繍绠楋紝鍙互浣跨敤鐭╅樀杩涜姹傝В銆傚浜榻愭绾挎ф柟绋嬫潵璁诧紝浣跨敤null(A,r)銆傚叾涓璻琛ㄧず浣跨敤绠鍖栭樁姊...
绛旓細MATLAB涓鐢ˋX=0鐨勬柟娉曟眰X鍙兘姹傚嚭0瑙 璇存槑绯绘暟鐭╅樀鐨勮鍒楀紡锛岃Matlab璁＄畻璁や负涓嶄负0锛岄渶瑕佽ˉ鍏呮潯浠
绛旓細鑻=0锛榻愭鏂圭▼锛夛紝鍒欑敤鍏嬭幈濮嗘硶鍒欐眰渚垮彲姹傚叾瑙ｏ紙姹傝鍒楀紡鐢ㄥ懡浠et())锛涜嫢f涓嶄负0锛堥潪榻愭锛夛紝鍒欑敤楂樻柉娑堝厓娉曪紙鍗冲垵绛夊彉鎹㈡硶锛夊彲姹傚叾瑙ｃMatlab涓鎻愪緵鍛戒护lu瀵鐭╅樀杩涜LU鍒嗚В锛屽鏋滄槸绋鐤忕煩闃碉紝鍒欏彲浣跨敤鍛戒护lunic瀵圭煩闃佃繘琛孡U鍒嗚В銆備綘鐨4鍏冧竴娆℃柟绋嬭嫢琛ㄧず鎴愶細AX=f 鍒橝=[a,b,c,d],X=[x,...
绛旓細Phi=[]; % 缃垵濮嬬殑桅(t)鐭╅樀涓虹┖锛屽悗闈緷娆℃坊鍔犲悇鍒 for j=1:n 閫氳繃寰幆渚濇娣诲姞桅(t)鐭╅樀鍚勫垪锛屽嚱鏁癳val绾睘澶氫綑 Phi=[Phi,eval(v(:,j)*exp(lambda(j,j)*t))];end 涓嬮潰杩欎竴鍙ュ張鏄帿鍚嶅叾濡欙紝濡傛灉鎯虫樉绀哄彉閲廝hi锛岀洿鎺ュ啓鍙橀噺鍚嶅嵆鍙紝娌″繀瑕佺粰鑷韩璧嬪 Phi=Phi 閫氳繃澶勪互桅(0)锛屽...
绛旓細1. 璁続x = b锛屾眰x锛岋紙x鍜宐閮芥槸鍚戦噺锛夊垯锛歺 = A\b 灏辫繖涔堢畝鍗曘傛兂涓嶅嚭鏇村鏉傜殑鏂规硶浜嗐2. 浣犱笂闈㈢殑鐪嬭捣鏉ユ槸涓榻愭鏂圭▼锛屼絾鏄洜涓篵1 b2 b3宸茬煡锛屽彲浠ョЩ鍒扮瓑鍙峰彟涓杈癸紱鍐嶆妸V3=V4鏁寸悊鍒鐭╅樀鏂圭▼閲屾垨鑰呭共鑴嗕粠鏂圭▼涓幓鎺塚4锛屽氨鍙互鐢ㄤ笂闈㈡柟娉曡В浜嗐
