solve seven nonlinear equations with seven variables in matlab
2 views (last 30 days)
Show older comments
hi, I have the following seven equations with seven varibles, I use syms to solve them
syms Eo a b t1 e1 t2 e2
ans = solve('(Eo*0.01)+(a*0.01^2)+(b*0.01^3)+(t1*e1*0.001*(1-exp(-0.1/t1)))+(t2*e2*(1-exp(-0.1/t2)))=30', ...
'(Eo*0.02)+(a*0.02^2)+(b*0.02^3)+(t1*e1*0.001*(1-exp(-0.2/t1)))+(t2*e2*(1-exp(-0.2/t2)))=62', ...
'(Eo*0.03)+(a*0.03^2)+(b*0.03^3)+(t1*e1*0.001*(1-exp(-0.3/t1)))+(t2*e2*(1-exp(-0.3/t2)))=86', ...
'(Eo*0.04)+(a*0.04^2)+(b*0.04^3)+(t1*e1*0.001*(1-exp(-0.4/t1)))+(t2*e2*(1-exp(-0.4/t2)))=101', ...
'(Eo*0.05)+(a*0.05^2)+(b*0.05^3)+(t1*e1*0.001*(1-exp(-0.5/t1)))+(t2*e2*(1-exp(-0.5/t2)))=108', ...
'(Eo*0.06)+(a*0.06^2)+(b*0.06^3)+(t1*e1*0.001*(1-exp(-0.6/t1)))+(t2*e2*(1-exp(-0.6/t2)))=111', ...
'(Eo*0.07)+(a*0.07^2)+(b*0.07^3)+(t1*e1*0.001*(1-exp(-0.7/t1)))+(t2*e2*(1-exp(-0.7/t2)))=111', ...
'Eo','a','b','t1','e1','t2','e2')
I get this answer
Warning: Explicit solution could not be found. > In solve at 83 In Untitled at 5 ans = [ empty sym ]
please help!
0 Comments
Accepted Answer
Walter Roberson
on 12 Apr 2013
There are two solution sets:
Eo = 7308.837897,
a = -111067.8168,
b = 569105.6911,
e1 = 56.07057628+300.3704267*1i,
e2 = -517.5135254,
t1 = -0.5741854156e-2+0.3075914851e-1*1i,
t2 = 0.9863006784e-1
Eo = 7308.837897,
a = -111067.8168,
b = 569105.6911,
e1 = -517513.5254,
e2 = 0.5607057628e-1+.3003704267*1i,
t1 = 0.9863006784e-1,
t2 = -0.5741854156e-2+0.3075914851e-1*1i
3 Comments
Walter Roberson
on 12 Apr 2013
Edited: Walter Roberson
on 12 Apr 2013
In Maple, I used
solve(['0.1e-1*Eo+a*0.1e-1^2+b*0.1e-1^3+0.1e-2*t1*e1*(1-exp(-.1/t1))+t2*e2*(1-exp(-.1/t2)) = 30', '0.2e-1*Eo+a*0.2e-1^2+b*0.2e-1^3+0.1e-2*t1*e1*(1-exp(-.2/t1))+t2*e2*(1-exp(-.2/t2)) = 62', '0.3e-1*Eo+a*0.3e-1^2+b*0.3e-1^3+0.1e-2*t1*e1*(1-exp(-.3/t1))+t2*e2*(1-exp(-.3/t2)) = 86', '0.4e-1*Eo+a*0.4e-1^2+b*0.4e-1^3+0.1e-2*t1*e1*(1-exp(-.4/t1))+t2*e2*(1-exp(-.4/t2)) = 101', '0.5e-1*Eo+a*0.5e-1^2+b*0.5e-1^3+0.1e-2*t1*e1*(1-exp(-.5/t1))+t2*e2*(1-exp(-.5/t2)) = 108', '0.6e-1*Eo+a*0.6e-1^2+b*0.6e-1^3+0.1e-2*t1*e1*(1-exp(-.6/t1))+t2*e2*(1-exp(-.6/t2)) = 111', '0.7e-1*Eo+a*0.7e-1^2+b*0.7e-1^3+0.1e-2*t1*e1*(1-exp(-.7/t1))+t2*e2*(1-exp(-.7/t2)) = 111'])
Note: the '' are not needed in Maple and just happen to not hurt either.
More Answers (0)
See Also
Categories
Find more on Surrogate Optimization in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!