Finding solution of an equation containing symbolic parameters. I want solution of one symbolic parameter in terms of another

Question

Aiswarya Nair Manikantan Geetha 2022-4-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1703280-finding-solution-of-an-equation-containing-symbolic-parameters-i-want-solution-of-one-symbolic-para

评论： Paul 2022-4-23

syms c_1 c_2 D r_t
%Parameters
sig = 0.2;
p=0.5;
rho_t =1;
w_t =100;
%Function
u = (c_1^sig)/sig ;
v = (c_2^sig)/sig;
u_diff = diff(u);
v_diff = diff(v);
%D Function
du_dd = subs(u_diff,w_t-D);
dv_dd = subs(v_diff,r_t*D);
eqn2 = du_dd == r_t*dv_dd ;
assume(D>0)
assume(D<w_t)
assume(r_t>0)
S2 = solve(eqn2,D,'Real',true)

I am trying to find soltion of D in terms of r_t but it throws an error

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Paul 2022-4-23

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1703280-finding-solution-of-an-equation-containing-symbolic-parameters-i-want-solution-of-one-symbolic-para#answer_949215

在 MATLAB Online 中打开

Not sure why solve doesn't work as is. See below for one approach to a solution.

syms c_1 c_2 D r_t
%Parameters
sig = 0.2;
p=0.5;
rho_t =1;
w_t =100;
%Function
u = (c_1^sig)/sig ;
v = (c_2^sig)/sig;
u_diff = diff(u);
v_diff = diff(v);
%D Function
du_dd = subs(u_diff,w_t-D);
dv_dd = subs(v_diff,r_t*D);
eqn2 = du_dd == r_t*dv_dd ;
assume(D>0)

Need use assumeAlso so as not to wipeout the first assumpton on D

assumeAlso(D<w_t)
assume(r_t>0)

Raise both sides of the equation to the 5/4 power and then try to solve. Not sure if doing so is strictly valid for all possible values of D and r_t given the assumptions, but I don't see why it wouldn't be.

S2 = solve((eqn2.^(5/4)),D,'Real',true,'ReturnConditions',true)
S2 = struct with fields:
             D: x
    parameters: x
    conditions: x^4 + 400*r_t*x^3 + 4000000*r_t*x == r_t*x^4 + 60000*r_t*x^2 + 100000000*r_t & r_t ~= 1 & 0 < x & x < 100

So we get a solution, but there are conditions. But if we simplify and then solve, we get one solution

S2 = solve(simplify(eqn2.^(5/4)),D)

S2 =

Verifty that's a solution

simplify(subs(eqn2,D,S2))
ans = symtrue

2 个评论
显示无隐藏无

Aiswarya Nair Manikantan Geetha 2022-4-23

Wow ! That did work. Can you tell me how u got the idea to raise equation to the 5/4th power.

I just want to improve my way of thinking. Thank you.

Paul 2022-4-23

eqn2 has terms raised the 4/5 power. Raising both sides eqn2 to the 5/4 power just seemed like a logical thing to do. I'm still not sure if there isn't some corner case where doing so wouldn't be valid.

请先登录，再进行评论。