求助！一个简单的编程解方程，请大家帮忙。

oxysword

请问，我想通过输入E1，E2，E3，E4，得到W5和W6。应该如何编程。

E1=？;E2=？;E3=？;E4=？;

>> W1=E1;W2=E1+E2;W3=W2+E3;W4=W3+E4;
>> A=[-0.5*(W2+W1),1;-0.5*(W3+W2),1;-0.5*(W4+W3),1];
>> Y=[E2;E3;E4]; A1=A'; P=A1*A; P1=inv(P);

P2=P1*A1*Y;

此时得出P2，P2为一个2*1的矩阵。如P2=[12;44]
然后令b=12，u=44（即P2的两个值分别赋给b和t）

最后计算出
W5=[W1-u/b]*exp(1)^[-b*(5-1)]+u/b
W6=[W1-u/b]*exp(1)^[-b*(6-1)]+u/b

请高手帮忙如何解决啊？

[ 本帖最后由 ChaChing 于 2009-5-6 08:38 编辑 ]

friendchj

help function

oxysword

。。。可以具体说说么？高手。。。

friendchj

% myfun.m
function [w5 w6]=myfun(E1,E2,E3,E4)
W1=E1;W2=E1+E2;W3=W2+E3;W4=W3+E4;
A=[-0.5*(W2+W1),1;-0.5*(W3+W2),1;-0.5*(W4+W3),1];
Y=[E2;E3;E4]; A1=A'; P=A1*A; P1=inv(P);
P2=P1*A1*Y;
W5=[W1-u/b]*exp(1)^[-b*(5-1)]+u/b
W6=[W1-u/b]*exp(1)^[-b*(6-1)]+u/b

cheno9468

原帖由 friendchj 于 2009-5-6 09:09 发表

                               
登录/注册后可看大图

% myfun.m
function [w5 w6]=myfun(E1,E2,E3,E4)
W1=E1;W2=E1+E2;W3=W2+E3;W4=W3+E4;
A=[-0.5*(W2+W1),1;-0.5*(W3+W2),1;-0.5*(W4+W3),1];
Y=[E2;E3;E4]; A1=A'; P=A1*A; P1=inv(P);
P2=P1*A1*Y;
W5=[W1-u/ ...

b和u好象没有被定义