matlab解方程出現sym,matlab解方程出現sym

2021-03-17 19:26:50 字數 4623 閱讀 5002

1樓:匿名使用者

修改如下:

clear all; clc

e=2*10^5;

v=0.3;

g=e/(2*(1+v));

lambda=(v*e)/((1+v)*(1-2*v));

epsilon11=1.7833*10^-3;

epsilon22=1.6*10^-3;

syms sigma11 sigma22 epsilon33;

eq1= sigma11==2*g*epsilon11+lambda*(epsilon11+epsilon22+epsilon33);

eq2= sigma22==2*g*epsilon22+lambda*(epsilon11+epsilon22+epsilon33);

eq3= epsilon33==-v*(sigma11+sigma22)/e;

[sigma11,sigma22,epsilon33]=solve(eq1,eq2,eq3);

sigma11=vpa(sigma11,8);

sigma22=vpa(sigma22,8);

epsilon33=vpa(epsilon33,8);

val_sigma11 = double(sigma11)

val_sigma22 = double(sigma22)

val_epsilon33 = double(epsilon33)

matlab2013b下執行通過,版本太低了別怪別人。

val_sigma11 =

-0.0014

val_sigma22 =

497.4286

val_epsilon33 =

469.2286

matlab解方程組,答案是sym沒有具體數值

2樓:

修改如下:

clear all; clc

e=2*10^5;

v=0.3;

g=e/(2*(1+v));

lambda=(v*e)/((1+v)*(1-2*v));

epsilon11=1.7833*10^-3;

epsilon22=1.6*10^-3;

syms sigma11 sigma22 epsilon33;

eq1= sigma11==2*g*epsilon11+lambda*(epsilon11+epsilon22+epsilon33);

eq2= sigma22==2*g*epsilon22+lambda*(epsilon11+epsilon22+epsilon33);

eq3= epsilon33==-v*(sigma11+sigma22)/e;

[sigma11,sigma22,epsilon33]=solve(eq1,eq2,eq3);

sigma11=vpa(sigma11,8);

sigma22=vpa(sigma22,8);

epsilon33=vpa(epsilon33,8);

val_sigma11 = double(sigma11)

val_sigma22 = double(sigma22)

val_epsilon33 = double(epsilon33)

matlab2013b下執行通過,版本太低了別怪別人。

val_sigma11 =

-0.0014

val_sigma22 =

497.4286

val_epsilon33 =

469.2286

3樓:宛丘山人

解方程組失效,可將程式改為:

clear

e=2*10^5;

v=0.3;

g=e/(2*(1+v));

lambda=(v*e)/((1+v)*(1-2*v));

epsilon11=1.7833*10^-3;

epsilon22=1.6*10^-3;

epsilon33=-2*v*(g+lambda)*(sigma11+sigma22)/(2*v+e)

sigma11=2*g*epsilon11+lambda*(epsilon11+epsilon22+epsilon33)

sigma22=2*g*epsilon22+lambda*(epsilon11+epsilon22+epsilon33)

最後執行結果是:

>> epsilon33=-2*v*(g+lambda)*(epsilon11+epsilon22)/(2*v+e)

epsilon33 =

-0.0020

>> sigma11=2*g*epsilon11+lambda*(epsilon11+epsilon22+epsilon33)

sigma11 =

439.5156

>> sigma22=2*g*epsilon22+lambda*(epsilon11+epsilon22+epsilon33)

sigma22 =

411.3156

matlab中使用solve命令出現sym問題

4樓:

function y=func_nonlinear(x)

a=[0.2252 0.2453 -0.1968 -0.1834];

b=[0.1665 -0.1506 -0.1757 0.0798];

xw1=0;xw2=0;xw3=1;xw4=1;

yw1=0.6;yw2=0;yw3=0;yw4=0.6;

%s=solve('xw1*r11+yw1*r12+a(1)*xw1*r31+a(1)*yw1*r32+t1+a(1)*t3=0',...

%'xw1*r21+yw1*r22+b(1)*xw1*r31+b(1)*yw1*r32+t2+b(1)*t3=0',...

%'xw2*r11+yw2*r12+a(2)*xw2*r31+a(2)*yw2*r32+t1+a(2)*t3=0',...

%'xw2*r21+yw2*r22+b(2)*xw2*r31+b(2)*yw2*r32+t2+b(2)*t3=0',...

%'xw3*r11+yw3*r12+a(3)*xw3*r31+a(3)*yw3*r32+t1+a(3)*t3=0',...

%'xw3*r21+yw3*r22+b(3)*xw3*r31+b(3)*yw3*r32+t2+b(3)*t3=0',...

%'xw4*r11+yw4*r12+a(4)*xw4*r31+a(4)*yw4*r32+t1+a(4)*t3=0',...

%'xw4*r21+yw4*r22+b(4)*xw4*r31+b(4)*yw4*r32+t2+b(4)*t3=0',...

%'r11^2+r21^2+r31^2=1',...

%'r21^2+r22^2+r32^2=1',...

%'r11*r12+r21*r22+r31*r32=0');

%r11 r12 r21 r22 r31 r32 t1 t2 t3

%分別為x(1),x(2),x(3),x(4),x(5),x(6),x(7),x(8),x(9)

y(1)=xw1*x(1)+yw1*x(2)+a(1)*xw1*x(5)+a(1)*yw1*x(6)+x(7)+a(1)*x(9);

y(2)=xw1*x(3)+yw1*x(4)+b(1)*xw1*x(5)+b(1)*yw1*x(6)+x(8)+b(1)*x(9);

y(3)=xw2*x(1)+yw2*x(2)+a(2)*xw2*x(5)+a(2)*yw2*x(6)+x(7)+a(2)*x(9);

y(4)=xw2*x(3)+yw2*x(4)+b(2)*xw2*x(5)+b(2)*yw2*x(6)+x(8)+b(2)*x(9);

y(5)=xw3*x(1)+yw3*x(2)+a(3)*xw3*x(5)+a(3)*yw3*x(6)+x(7)+a(3)*x(9);

y(6)=xw3*x(3)+yw3*x(4)+b(3)*xw3*x(5)+b(3)*yw3*x(6)+x(8)+b(3)*x(9);

y(7)=xw4*x(1)+yw4*x(2)+a(4)*xw4*x(5)+a(4)*yw4*x(6)+x(7)+a(4)*x(9);

y(8)=xw4*x(3)+yw4*x(4)+b(4)*xw4*x(5)+b(4)*yw4*x(6)+x(8)+b(4)*x(9);

y(9)=x(1)^2+x(3)^2+x(5)^2-1;

y(10)=x(3)^2+x(4)^2+x(6)^2-1;

y(11)=x(1)*x(2)+x(3)*x(4)+x(5)*x(6);

將該函式儲存成m檔案,檔名為func_nonlinear.m

在命令視窗中輸入:

x0=[0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5];

fsolve(@func_nonlinear,x0)

顯示結果為

ans =

columns 1 through 8

0.8617 -0.0013 0.

1270 -0.9599 0.4908 0.

2509 -0.4245 0.2616

column 9

1.7325

這就是9個變數的迭代數值解。

我的matlab版本是7.1的,不知道是否適合其他版本,呵呵

解方程x305,解方程x

x 3 0.5 12 x 2.5 12 x 12 2.5 x 9.5 x 12 2.5 x 9.5 x 12 3 0.5 9 0.5 9.5 解方程 x 3 0.5 12 5x x 3 35 x 3 0.5 12 x 2.5 12 x 12 2.5 x 9.5 5 x 3 35 x 3 35 5 x...

解方程 要帶過程!解方程,要過程!

3 10 x 1 25 1 4 x 3 10 1 25 1 4 x 30 100 4 100 25 100x 1 100 您好,很高興為您解答,skyhunter002為您答疑解惑。如果本題有什麼不明白可以追問,如果滿意記得。如果有其他問題本題後另發點選向我求助,答題不易,請諒解,謝謝。祝學習進步。...

解方程3234x,解方程3234x

3 4x 11 8 x 11 6 2x 30 46 5 2x 8.9 x 4.45 x 2 3 3 4 x 3 4 2 3 x 1 2 3 4x 33 2 4x 2 11 x 1 22 解方程3 4x 1 8 1 2 3 4x 1 2 1 8 3 4x 3 8 x 3 8 3 4 x 1 2 3 4...