Substituting an expression not giving intended result

Question

John 2019-1-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438653-substituting-an-expression-not-giving-intended-result

评论： Star Strider 2019-1-7

I'm trying to simplify the following expression using the subs function. In the code below, I was expecting Asubs = 3*Phi. Any ideas on what might be wrong in my code?

clear all;
clc;
syms phi(c,cbar,sigma) Phi(c,cbar,sigma) A(csol,cbar,sigma)
assumeAlso(sigma,'real')
assumeAlso(sigma>0)
phi(c,cbar,sigma)= (1/sqrt(2*sym(pi)*sigma^2))*exp(-(c-cbar)^2/(2*sigma^2))
A(csol,cbar,sigma)=3*int(phi(c,cbar,sigma),c,-inf,csol)
Phi(csol,cbar,sigma) = 1/2 + 1/2*erf(sqrt(2)*(csol-cbar)/(2*sigma))
%Phi(x,sigma) = 1/2 + 1/2*erf(sqrt(2)*(x)/(2*sigma))
Asubs = subs(A,Phi)

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-1-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438653-substituting-an-expression-not-giving-intended-result#answer_355375

在 MATLAB Online 中打开

Since ‘Phi’ does not appear as such specifically in ‘A’, the substitution will fail.

Perhaps this is closer to what you want:

As = simplify(A/Phi)

prodcing:

As(csol, cbar, sigma) =
3

2 个评论
显示无隐藏无

John 2019-1-7

I see. The code in the question is a kind of mwe for the following problem. The normal distribution CDF and PDF appear in many of the formulas that I'm working. I wanted to shorten those formulas by using $\Phi$ and $\phi$ to represent them in the expressions. In case you know of a workaround, I'd love to hear.

Star Strider 2019-1-7

I am still at a loss to understand what you are doing.

I would not use the Symbolic Math Toolbox for such calculations. It is primarily designed for derivations and other such. Use the matlabFunction (link) function to create anonymous functions from your symbolic functions. Then use them in your calculations.

请先登录，再进行评论。

Answer 2

madhan ravi 2019-1-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438653-substituting-an-expression-not-giving-intended-result#answer_355370

在 MATLAB Online 中打开

syms c cbar sigma Phi A csol
assumeAlso(sigma,'real')
assumeAlso(sigma>0)
phi(c,cbar,sigma)= (1/sqrt(2*sym(pi)*sigma^2))*exp(-(c-cbar)^2/(2*sigma^2))
A(csol,cbar,sigma)=3*int(phi(c,cbar,sigma),c,-inf,csol)
Phi(csol,cbar,sigma) = 1/2 + 1/2*erf(sqrt(2)*(csol-cbar)/(2*sigma))
%Phi(x,sigma) = 1/2 + 1/2*erf(sqrt(2)*(x)/(2*sigma))
Asubs = vpa(subs(A,Phi)) % cbar csol and sigma values are not known