Substituting an expression not giving intended result
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponse acceptée
Star Strider
le 7 Jan 2019
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 commentaires
Star Strider
le 7 Jan 2019
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.
Plus de réponses (1)
madhan ravi
le 7 Jan 2019
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
2 commentaires
Voir également
Catégories
En savoir plus sur Assumptions dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!