用matlab求解多元一次方程组 求大神指导 急!!急!! 急!! 救命啊!! matlab数学实验中'*什么意思?

\u5728\u6570\u5b66matlab\u600e\u4e48\u5b9a\u4e49\u51fd\u6570

\u9996\u5148\u5efa\u7acbM\u6587\u4ef6\u6216\u76f4\u63a5\u70b9\u51fb\uff08File/New/Function)\u5efa\u7acb\u51fd\u6570\u6587\u4ef6\uff0c\u5176\u4e2d\u51fd\u6570\u6587\u4ef6\u7684\u683c\u5f0f\u662f\uff1a
function [\u8f93\u51fa\u53d8\u91cf] = \u51fd\u6570\u540d\u79f0(\u8f93\u5165\u53d8\u91cf\uff09
% \u6ce8\u91ca
% \u51fd\u6570\u4f53
\u4e3e\u4f8b\u8bf4\u660e\uff1a\u5047\u8bbe\u6709\u4e2a\u51fd\u6570 y=2*x^2+3*x+4\uff0c\u7f16\u5199\u4e00\u4e2a\u5b50\u51fd\u6570\uff08\u81ea\u5b9a\u4e49\u51fd\u6570\uff09
\u5efa\u7acb1\u4e2a*.m\u6587\u4ef6\uff0c\u6587\u4ef6\u5185\u5bb9\u5982\u4e0b
function y=f(x) %x\u4e3a\u5b50\u51fd\u6570\u8f93\u5165\uff0cy\u4e3a\u8f93\u51fa\uff0c\u8f93\u5165\u8f93\u51fa\u53ef\u4ee5\u662f\u591a\u4e2a\uff0c\u4e5f\u53ef\u4ee5\u662f\u6570\u7ec4\u6216\u77e9\u9635
y=2*x.^2+3*x+4

\u4fdd\u5b58\u6b64m\u6587\u4ef6\u540d\u4e3ay
\u8c03\u7528\u65f6
x=2; y=f(x);
\u6216\u8005
x=[1 2 3];y=f(x)

'\u4ee3\u8868\u8f6c\u7f6e,*\u4ee3\u8868\u4e58\u6cd5.
r=0:0.05:1\u8868\u793ar\u662f\u4e00\u4e2a\u884c\u5411\u91cf\uff1b
r'*cos(t)\u8868\u793ar\u8f6c\u7f6e\u540e\uff08\u4e3a\u5217\u5411\u91cf\uff09\u518d\u4e58\u4ee5cos(t)\uff01

1、配料问题一半应该有目标函数吧?最好把原始问题贴出来。

当然,如果没有目标函数,也可以自己随便设定,能够影响到的是,求出的究竟是哪一个解。

 

2、这可以看作一个优化问题来求解,而不是解方程组。要想求出【唯一的一组解】往往是做不到的,因为这取决于你给的限定条件本身是否具有唯一解。具有唯一解的几何意义是N维空间中的一个点,这种可能远远小于有多个解(空间的一个区域)或者无解(约束条件存在冲突)的情况。

 

3、就题主现在的问题而言,就是无解的,这可以用线性规划的方法验证:

Ae=[2 3 5 6; 1 1 1 1];
be=[89; 20];
lb=[5 0 0 0];
ub=[10 5 5 5];
x=linprog([1 0 0 0],[],[],Ae,be,lb,ub)

得到结果是无法找到可行解。

 

上述结论也可以通过Lingo进行验证:

max=x1;
2*x1+3*x2+5*x3+6*x4=89;
x1+x2+x3+x4=20;
5<x1;
x1<10;
0<x2;
x2<5;
0<x3;
x3<5;
0<x4;
x4<5;

用Lingo的好处是,可以通过Debug来判断哪些约束条件造成冲突:

  Constraints and bounds that cause an infeasibility:
  Sufficient Rows:
  (Dropping any sufficient row will make the model feasible.)
   [_3] X1 + X2 + X3 + X4 = 20 ;
   [_2] 2 * X1 + 3 * X2 + 5 * X3 + 6 * X4 = 89 ;
   [_9] X3 <= 5 ;
   [_11] X4 <= 5 ;

  Necessary Variable Bounds:
  (If none of the necessary and sufficient bounds are dropped,
   then the model remains infeasible.)
   X1 >=  0

例如,如果上面的条件去掉x3<5或x4<5这两个约束中的一个,问题就是有解的。



  • matlab鎬庝箞姹傚鍏冧竴娆℃柟绋缁
    绛旓細鍙互鐢╲pasolve姹傝В銆傚疄鐜颁唬鐮侊細for lambda=1:0.1:2 syms x qr=1.449*lambda.*(1-0.1416*lambda.^2).^3.0303;lambda1=vpasolve(1.57744*x*(1-0.1667*x.^2).^2.5==qr)end 杩愯缁撴灉 matlab瑙f柟绋缁刲nx琛ㄧず鎴恖og(x)鑰宭gx琛ㄧず鎴恖og10(x)1-exp(((log(y))/x^0.5)/(x-1))1...
  • matlab姹傝В澶氬厓鏂圭▼缁 鍏厓鏂圭▼缁 姹傚ぇ浣В绛
    绛旓細姹傝В澶氬厓涓娆℃柟绋缁勶紝涓鑸彲浠ョ敤vpasolve锛堬級鍑芥暟鏉ヨВ鍐炽1銆佸x锛宎锛 b锛 c锛 d锛 e鍙橀噺杩涜澹版槑锛屽嵆 syms x a b c d e 2銆佸皢绛夊紡璧嬪肩粰eq1锛宔q2锛宔q3锛宔q4锛宔q5锛宔q6锛屽嵆 eq1=a+c+d==5.7725;eq2=b+2*d+e == 6.2196;eq3=a+2*c+e == 4.1668+0.0889*x;eq4=b*c/(...
  • 濡備綍鐢╩atlab鏉姹傝В澶氬厓涓娆鐜版鏂圭▼
    绛旓細x+2*y=5,3*x-y=7 杩欐槸涓涓簩鍏涓娆绾挎鏂圭▼缁 鎵浠鍦╩atlab涓寜鐓т笅闈㈣緭鍏ワ細>> [x y]=solve('x+2*y=5','3*x-y=7') %琛ㄧず姹倄鍜寉 鍦╯olve鍛戒护閲岃緭鍏ヤ袱涓柟绋嬪嵆鍙 x = 19/7 y = 8/7 涓夊厓绾挎,鍥涘厓绾挎 閮藉彲浠ユ眰,鍙渶瑕佽緭鍏ュ洓涓柟绋嬪氨鏄,姹傝В鐨勫彉閲忔槸x y u v 鍥涗釜鍙橀噺 ...
  • 鐢╩atlab鎬庝箞姹傝В澶氬厓涓娆寰垎鏂圭▼缁
    绛旓細,yN.濡傛灉鏈夊垵濮嬫潯浠讹紝鍙互鎶婃潯浠朵竴璧蜂紶缁欏嚱鏁版潵瀹氳В锛歔y1,,yN]=dsolve(eqns,conds)solvesthesystemofordinarydifferentialequationseqnswiththeinitialorboundaryconditionsconds.缁欏嚭涓涓2涓彉閲忕殑寰垎鏂圭▼缁姹傝В浠g爜锛歴ymsx(t)y(t)z=dsolve(diff(x)==y,...
  • matlab姹傝В澶氬厓鏂圭▼缁
    绛旓細鐢╩atlab姹傝В澶氬厓鏂圭▼缁勶紝鍙互鎶婇棶棰樿浆鎹负鏋佸奸棶棰樸傚嵆 min 3+sum(p.*log2(p))st. sum(p)=1 -sum(p.*log2(p))=3 p1,p2,p3,...p8鈮0 鑰屾瀬鍊奸棶棰樼敤fmincon鏈灏忓煎嚱鏁版潵姹傝В銆傜紪绋嬪悗鍙互寰楀埌濡備笅缁撴灉锛歱1=0.125005874372445锛沺2=0.125007579745882锛沺3=0.124999277594269锛沺4=0.124998762...
  • matlab濡備綍姹傝В涓嬪垪澶氬厓涓娆℃柟绋,绋嬪簭鍙婅繍琛岀粨鏋溿
    绛旓細绋嬪簭澶ц嚧濡備笅銆備絾鏄垜娌℃湁寰楀埌瑙c傜1姝.鐩爣鍑芥暟 鍦∕atlab涓嬭緭鍏ワ細edit锛岀劧鍚庡皢涓嬮潰涓よ鐧惧垎鍙蜂箣闂寸殑鍐呭锛屽鍒惰繘鍘伙紝淇濆瓨 function y=zhidao_bnb20_NGZHUMENG_opt(x)d,e,f,N鍒嗗埆鐢▁(1),x(2),x(3),x(4)浠f浛 y=x(4);绗2姝.绾︽潫鏉′欢 鍦∕atlab涓嬭緭鍏ワ細edit锛岀劧鍚庡皢涓嬮潰涓よ鐧惧垎鍙蜂箣闂寸殑...
  • 濡備綍鍒╃敤matlab姹傝В澶氬厓涓娆℃柟绋缁勭殑鍚勭郴鏁
    绛旓細瀹為檯鍙互solve姹傝В锛屽 Y1='a1*x+b1*y=c1'Y2='a2*x+b2*y=c2'Y3='a3*x+b3*y=c3'[x,y]=solve(Y1,Y2,Y3)
  • 鎬ユ眰,matlab姹傝В澶氬厓涓娆℃柟绋
    绛旓細鍏釜鍙橀噺锛屽叓涓鏂圭▼銆備笅娆″啓鐨勬椂鍊欙紝涓嶈鍑虹幇 琛ㄨ揪寮1=琛ㄨ揪寮2=琛ㄨ揪寮3 鍙﹀鍙橀噺涓庢暟瀛椾箣闂磋鏈*锛屾垜鏀逛簡鑰佸崐澶┿傚懡浠:锛堜綘浠渶濂藉啀妫鏌ヤ竴涓嬫湁娌℃湁鍐欓敊锛塠aa,bb,cc,dd,ee,ff,pp,qq]=solve('a+(p-2)*b+c=280','5580/1.2+(a+0.5)*2400/0.6+1000+40*a/1.2+11200/1.2-(a+0...
  • matlab澶氬厓涓娆℃柟绋缁姹傝В,绋嬪簭鍑洪敊,鍘熷洜
    绛旓細绗竴绉嶄篃灏辨槸浣犳兂杈撳叆鐨勬柟娉曪細syms x y z;%娉ㄦ剰鍦ㄧ敤syms瀹氫箟澶氫釜鍙橀噺鏃朵腑闂翠笉鑳界敤閫楀彿锛岃岃鐢ㄧ┖鏍硷紱[x,y,z]=solve(2*x+3*y+z-450000,x+2*y+z-300000,2*x+2*y+z-350000,x,y,z)绗簩绉嶆柟娉曟槸涓嶇敤瀹氫箟绗﹀彿鍙橀噺锛歝lear x y z;[x,y,z]=solve('2*x+3*y+z=450000','x+...
  • 鎬庝箞鐢∕atlab鏉姹傝В澶氬厓涓娆℃柟绋鍟?
    绛旓細clc;clear;y=sym('a+b*x+c*x^2+d*x^3+e*x^4+f*x^5')%鏄笉鏄痗*x^2鍟 y1=diff(y)y2=diff(y,2)y3=diff(y,3)eq1=subs(y,'x',7)eq2=subs(y1,'x',7)eq3=subs(y2,'x',7)eq4=subs(y,'x',0)-6.5 eq5=subs(y1,'x',0)eq6=subs(y3,'x',0)X=solve(eq1...
  • 扩展阅读:求解方程计算器 ... matlab多元一次方程 ... matlab求解多元方程组 ... matlab求解器在哪 ... matlab x的n次方怎么打 ... matlab求累加和 ... matlab求2的1到m次方之和 ... 和matlab有一样功能的软件 ... matlab求1到100的和 ...

    本站交流只代表网友个人观点,与本站立场无关
    欢迎反馈与建议,请联系电邮
    2024© 车视网