How to find the expression for theta1?
2 ビュー (過去 30 日間)
古いコメントを表示
theta1=2*atan((2*tan(theta2/2)*tan(gama/2 - theta2/2) + (4*tan(theta2/2)^2 - 4*M^2*tan(theta2/2)^2 - 4*M^2*tan(gama/2 - theta2/2)^2 - M^2*gama^2*l^2 - 8*M^2*tan(theta2/2)*tan(gama/2 - theta2/2) - 4*M^2*gama*l*tan(theta2/2) - 4*M^2*gama*l*tan(gama/2 - theta2/2) - M^2*gama^2*l^2*tan(theta2/2)^2*tan(gama/2 - theta2/2)^2 + 2*M^2*gama^2*l^2*tan(theta2/2)*tan(gama/2 - theta2/2) + 4*M^2*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)^2 + 4*M^2*gama*l*tan(theta2/2)^2*tan(gama/2 - theta2/2) + 4)^(1/2) - 2)/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)) - (2*tan(theta2/2)*tan(gama/2 - theta2/2))/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)))
Can anyone help me to find an expression for theta1 in terms of thetac where
theta2=thetac-theta1
1 件のコメント
Manuela Gräfe
2017 年 4 月 24 日
Hi, umme mumtahina.
I see you are working with the LLC converter and the IEEE document (Optimal design methodology for LLC Resonant Converter... by Zhijian Fang etc.).
I am looking for the same solution at the moment for my bachelor thesis and I was wondering if you could provide me your MATLAB code? So I 'don't have to annoy Walter Roberson with the same issues. Please contact me via private message.
回答 (2 件)
John D'Errico
2017 年 4 月 22 日
Nope. Wanting magic to happen is not sufficient. I am quite confident that there is no way to untangle the implicit function you will get when you eliminate theta2.
Having theta1 inside and out of all of those terms will make it impossible to solve.
If you have values for all of the other variables, then you can use a root finder, but expect multiple solutions to exist.
I'd suggest rethinking your problem.
4 件のコメント
Walter Roberson
2017 年 4 月 23 日
syms theta1 theta2 gama M l
eqn = theta1 == 2*atan((2*tan(theta2/2)*tan(gama/2 - theta2/2) + (4*tan(theta2/2)^2 - 4*M^2*tan(theta2/2)^2 - 4*M^2*tan(gama/2 - theta2/2)^2 - M^2*gama^2*l^2 - 8*M^2*tan(theta2/2)*tan(gama/2 - theta2/2) - 4*M^2*gama*l*tan(theta2/2) - 4*M^2*gama*l*tan(gama/2 - theta2/2) - M^2*gama^2*l^2*tan(theta2/2)^2*tan(gama/2 - theta2/2)^2 + 2*M^2*gama^2*l^2*tan(theta2/2)*tan(gama/2 - theta2/2) + 4*M^2*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)^2 + 4*M^2*gama*l*tan(theta2/2)^2*tan(gama/2 - theta2/2) + 4)^(1/2) - 2)/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)) - (2*tan(theta2/2)*tan(gama/2 - theta2/2))/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)));
F=1.05;
gama=pi/F
thetac=2.98;
theta2=thetac-theta1;
l=0.218;
M=0.95;
seqn = subs(eqn);
sol = vpasolve(seqn);
or
F=1.05;
gama=pi/F
thetac=2.98;
l=0.218;
M=0.95;
eqn = @(theta1) (2*atan((-2*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)+(4*tan(-(1/2)*thetac+(1/2)*theta1)^2-4*M^2*tan(-(1/2)*thetac+(1/2)*theta1)^2-4*M^2*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)^2-M^2*gama^2*l^2+8*M^2*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)+4*M^2*gama*l*tan(-(1/2)*thetac+(1/2)*theta1)-4*M^2*gama*l*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)-M^2*gama^2*l^2*tan(-(1/2)*thetac+(1/2)*theta1)^2*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)^2-2*M^2*gama^2*l^2*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)-4*M^2*gama*l*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)^2+4*M^2*gama*l*tan(-(1/2)*thetac+(1/2)*theta1)^2*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)+4)^(1/2)-2)/(-2*tan(-(1/2)*thetac+(1/2)*theta1)+2*M*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)-2*M*tan(-(1/2)*thetac+(1/2)*theta1)+M*gama*l+M*gama*l*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1))+2*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)/(-2*tan(-(1/2)*thetac+(1/2)*theta1)+2*M*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)-2*M*tan(-(1/2)*thetac+(1/2)*theta1)+M*gama*l+M*gama*l*tan(-(1/2)*thetac+(1/2)*theta1)*tan((1/2)*gama-(1/2)*thetac+(1/2)*theta1)))) - theta1;
sol = fzero(eqn, rand())
Manuela Gräfe
2017 年 4 月 24 日
Hi, umme mumtahina.
I see you are working with the LLC converter and the IEEE document (Optimal design methodology for LLC Resonant Converter... by Zhijian Fang etc.).
I am looking for the same solution at the moment for my bachelor thesis and I was wondering if you could provide me your MATLAB code? So I 'don't have to annoy Walter Roberson with the same issues. Please contact me via private message.
0 件のコメント
参考
カテゴリ
Help Center および File Exchange で Call Python from MATLAB についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!