how to find numerical soloution
1 Ansicht (letzte 30 Tage)
Ältere Kommentare anzeigen
Hi everyone; How can I find the numerical solution of the equation: dldk3=0 in order to get k2? dldk3 is obtained from following codes:
clear all
clc
syms k1 k2 k3 u1 u2
m1=100*10^-9;
m2=100*10^-9;
a=200;
b=1;
s1=0;
s2=0;
x=(k1+k2)/m1+(k2+k3)/m2;
y=(((k1+k2)/m1)^2+((k2+k3)/m2)^2+(2*k2^2-(k1+k3)*k2-2*k1*k3)/(m1*m2))^0.5;
omega1=(2^0.5/2)*(x+y)^0.5;
omega2=(2^0.5/2)*(x-y)^0.5;
delta=omega2-omega1;
eps=(k2)/(y*(m1*m2)^0.5);
L=eps+delta+u1*(a-(2^0.5/2)*(((x+y)^0.5)-((x-y)^0.5))+s1^2)+u2*(b-(k2/(y*(m1*m2)^0.5))+s2^2);
dldk1=jacobian(L,k1);
dldk2=jacobian(L,k2);
dldk3=jacobian(L,k3);
dldu1=jacobian(L,u1);
dldu2=jacobian(L,u2);
B=solve(dldk1,'u1');
u1=B;
dldk2=subs(dldk2);
dldk3=subs(dldk3);
C=solve(dldk2,'u2');
u2=C;
dldu1=subs(dldu1);
dldk3=subs(dldk3);
D=solve(dldu1,'k1');
k1=D(1,1);
dldu2=subs(dldu2);
dldk3=subs(dldk3);
E=solve(dldu2,'k3');
k3=E(1,1);
dldk3=subs(dldk3);
0 Kommentare
Antworten (0)
Siehe auch
Kategorien
Mehr zu Calculus finden Sie in Help Center und File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!