0
私には4つの未知数(2x2)を持つ関数があります。 私は解決策を得ました。しかし、私がやろうとしていることは、いくつかのパラメータを変更して、元のソリューションがどのように変化するかを確認することです。fsolveでループを使用し、結果を関数にプラグインする
しかし、matlabはA(I)= B、 'I'と要素Bの数が同じでなければならないと言っています。
コードの例は次のようなものです(psi01 psi02 psi03 psi04はさまざまなパラメータです)。
R=902;
psi01=0.9:-0.1:0.1;
psi11=0.9:-0.1:0.1;
psi02=0.1:0.1:0.9;
psi12=0.1:0.1:0.9;
E0=[R/4 R/4; R/4 R/4];
for i=1:9
[email protected](E) [
g/E(2,1)-(th1+psi12(i)/psi11(i)*th2)-((psi01(i)/psi11(i))+a*b)/b*(g/E(1,1)-(rho*b* (psi11(i)/psi01(i))/(psi01(i)/psi11(i)+a*b))*(th1+psi12(i)/psi11(i)*th2)*(th1+psi02(i) /psi11(i)*th2) );
g/E(2,2)-(psi11(i)/psi12(i)*th1+th2)-((psi02(i)/psi12(i))+a*b)/b*(g/E(1,2)-(rho*b*(psi12(i)/psi02(i))/(psi02(i)/psi12(i)+a*b))*(psi11(i)/psi12(i)*th1+th2)*(psi01(i)/psi02(i)*th1+th2) );
R-(E(1,1)+E(1,2)+E(1,1)+E(2,2));
E(1,2)- c(E)*E(1,1);
E(2,2)- d(E)*E(2,1)
];
Ep2(i)=fsolve(vv, E0);
end