1樓：我行我素
按上面解方程的語句，k,delta都是未知數，一個方程，2個未知數，要求出以其中一個是自變數，另一個是函式，要這樣：
>> syms k delta
>> solve(power((5*sqrt(5)*exp(131/512)/16),k)+power((13*exp(87/512)/16),k)-delta,k)
warning: cannot find explicit solution.
> in solve (line 316)
ans =
empty sym: 0-by-1
上式不能解出，可把k指定一個數，如：
solve(power((5*sqrt(5)*exp(131/512)/16),3)+power((13*exp(87/512)/16),3)-delta)
ans =
3666239720380227/2251799813685248
>> solve(power((5*sqrt(5)*exp(131/512)/16),k)+power((13*exp(87/512)/16),k)-delta,delta)
ans =
(31754343866723/35184372088832)^k + (1084224671755031/1125899906842624)^k
matlab解方程出現sym 20
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
matlab中使用solve命令出現sym問題
3樓：
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的，不知道是否適合其他版本，呵呵
matlab解方程組，答案是sym沒有具體數值
4樓：
修改如下：
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
5樓：宛丘山人
解方程組失效，可將程式改為：
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批量求解方程，並繪製曲面圖，求高手
t 0 30 600 p 1000 2000 40000 l zeros length p length t for ii 1 length t for jj 1 length p a t ii p jj fun l l 2 log l 5 a l jj,ii fsolve fun,1000 end...
這道題用解方程怎麼做，這道題解方程解方程怎麼做
觀察 2，3，6，8，8，4，2，8，6，8，8，4，2，8 發現 第3個數開始，每6個數為相同數。解 n 2 6的餘數為1，2，3，4，5，0時，分別對應6，8，8，4，2，8.n 3的整數 90 2 6 14餘4 答 第90個數是4。 按照數列的產生方法，我們不難寫出數列的前幾項為 2，3，6，...
小學解方程應用題，小學解方程應用題及答案
1 用邊長為9dm和邊長為7dm和正方形的瓷磚輔都正好輔完。說明這塊地板的邊長起碼是9和7的公倍數，即為63dm2 用邊長為9dm的正方形瓷磚鋪地，共需 63 63 9 9 7 7 49塊。3 用邊長為7dm的正方形瓷磚鋪地，共需 63 63 7 7 9 9 81塊。小學解方程應用題及答案 小東 小...