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...