Solve equation and standalone program
显示 更早的评论
Hi ! I created a matlab code that uses the solve function provided by the Symbolic Math Toolbox. The code works fine. Then I want to create a standalone program using the matlab compiler. Apparently the Symbolic Math Toolbox can't be included during the compilation. Do you have an idea, advice... ?
Thanks.
Regards
采纳的回答
Walter Roberson
2013-12-12
0 个投票
If the equations are pre-determined and only coefficients need to be determined, then you can use matlabFunction() to generate the MATLAB Code into a file. Then set up your program to not generate the equations and instead call upon the pre-saved MATLAB to determine the results. This second program can then be compiled.
12 个评论
hamza smair
2017-6-14
Hello Sir, I'm having the same problem. Can you give me an example on how to use the matlabFunction and how to call it in my code? Thank you.
Walter Roberson
2017-6-14
if ~isdeployed()
syms x y
expr = solve(x^3 - 2*x + y, x);
fun = matlabfunction(expr, 'vars', y, 'file', 'abbracordovidae.m');
else
fun = @abbracordovidae;
end
Bryan Law
2021-2-1
Hi Walter, I've been trying to use the matlab function as well as ~isdeployed() as shown from your example above but have not been having luck. Could you assist me with me code as shown below please.
In summary, I'm trying to run the below code within a standalone program
% solve for x
syms x
x1 = solve(slope*x + yint == perpslope*x + perpyint,x);
This is what I have tried with the ~isdeployed and matlabFunction:
if ~isdeployed()
syms x
exp = solve(slope*x + yint == perpslope*x + perpyint,x)
x1 = matlabFunction(exp, x);
end
Thank you so much
Walter Roberson
2021-2-1
if ~isdeployed()
syms x
sol = solve(slope*x + yint == perpslope*x + perapyint, x)
x1fun = matlabFunction(sol, 'vars', {x}, 'File', 'x1_externalized.m', 'optimize', false);
else
%recall the previously saved function. Hope it is up to date with the latest parameters.
x1fun = @x1_externalized;
end
x1 = x1fun(actual_x_value);
However, it is more likely that you would want to generate it with more parameters. Anything that is not constant at that point should be made a symbolic variable and added to the 'vars' list at the time the function is generated.
It is true, though, that in some cases making a value from a numeric parameter into a symbolic parameter can result in solve() not being able to find a solution. For example, if your expression included x^perpexp and perpexp was numerically -4, -3, -2, -1, 0, 1, 2, 3, or 4, then solve() would be able to solve into explicit numeric formulas, but that if you made perpexp symbolic then solve() would not be able to come up with an answer.
João Barbosa
2022-12-21
编辑：João Barbosa
2022-12-21
I have the same problem can you please help me?
syms t80
for i=1:t
eq_FASE= thetai_FASE+(t80/(xd_FASE*E_FASE))-((w_FASE^2*Leq_FASE(i)^2)/(24*xd_FASE*S_FASE^2*t80^2)) == thetak_FASE(i) +(tmax_FASE/(xd_FASE*E_FASE))-((m_FASE(i)^2*w_FASE^2*Leq_FASE(i)^2)/(24*xd_FASE*S_FASE^2*tmax_FASE^2));
eq_CG= thetai_CG+(t80/(xd_CG*E_CG))-((w_CG^2*Leq_CG(i)^2)/(24*xd_CG*S_CG^2*t80^2)) == thetak_CG(i) +(tmax_CG/(xd_CG*E_CG))-((m_CG(i)^2*w_CG^2*Leq_CG(i)^2)/(24*xd_CG*S_CG^2*tmax_CG^2));
xpto_FASE=double(solve(eq_FASE,t80));
xpto_CG=double(solve(eq_CG,t80));
idx_FASE=xpto_FASE>0;
idx_CG=xpto_CG>0;
ti_FASE(i)=xpto_FASE(idx_FASE);
ti_CG(i)=xpto_CG(idx_CG);
T_FASE(i)=double(ti_FASE(i)*S_FASE);
T_CG(i)=double(ti_CG(i)*S_CG);
p_FASE(i)=double(T_FASE(i)/w_FASE);
p_CG(i)=double(T_CG(i)/w_CG);
end
syms t80
syms thetai_FASE xd_FASE E_FASE w_FASE Leq_FASE S_FASE thetak_FASE tmax_FASE m_FASE
syms thetai_CG xd_CG E_CG w_CG Leq_CG S_CG thetak_CG tmax_CG m_CG
eq_FASE= thetai_FASE+(t80/(xd_FASE*E_FASE))-((w_FASE^2*Leq_FASE^2)/(24*xd_FASE*S_FASE^2*t80^2)) == thetak_FASE +(tmax_FASE/(xd_FASE*E_FASE))-((m_FASE^2*w_FASE^2*Leq_FASE^2)/(24*xd_FASE*S_FASE^2*tmax_FASE^2));
eq_CG= thetai_CG+(t80/(xd_CG*E_CG))-((w_CG^2*Leq_CG^2)/(24*xd_CG*S_CG^2*t80^2)) == thetak_CG +(tmax_CG/(xd_CG*E_CG))-((m_CG^2*w_CG^2*Leq_CG^2)/(24*xd_CG*S_CG^2*tmax_CG^2));
t80_FASE = solve(eq_FASE,t80,'MaxDegree',3);
t80_FASE = matlabFunction(t80_FASE)
t80_CG = solve(eq_CG,t80,'MaxDegree',3);
t80_CG = matlabFunction(t80_CG)
%Now call t80_FASE and t80_CG in a loop with the input parameters of your choice.
sol(1:3) = t80_FASE(rand,rand,rand,rand,rand,rand,rand,rand,rand)
%and select the real root(s)
sol = sol(abs(imag(sol))<1e-6)
João Barbosa
2022-12-22
if ~isdeployed()
syms t80
syms thetai_FASE1 xd_FASE1 E_FASE1 w_FASE1 Leq_FASE1 S_FASE1 thetak_FASE1 tmax_FASE1 m_FASE1
syms thetai_CG1 xd_CG1 E_CG1 w_CG1 Leq_CG1 S_CG1 thetak_CG1 tmax_CG1 m_CG1
eq_FASE= thetai_FASE1+(t80/(xd_FASE1*E_FASE1))-((w_FASE1^2*Leq_FASE1^2)/(24*xd_FASE1*S_FASE1^2*t80^2)) == thetak_FASE1 +(tmax_FASE1/(xd_FASE1*E_FASE1))-((m_FASE1^2*w_FASE1^2*Leq_FASE1^2)/(24*xd_FASE1*S_FASE1^2*tmax_FASE1^2));
eq_CG= thetai_CG1+(t80/(xd_CG1*E_CG1))-((w_CG1^2*Leq_CG1^2)/(24*xd_CG1*S_CG1^2*t80^2)) == thetak_CG1 +(tmax_CG1/(xd_CG1*E_CG1))-((m_CG1^2*w_CG1^2*Leq_CG1^2)/(24*xd_CG1*S_CG1^2*tmax_CG1^2));
t80_FASE = solve(eq_FASE,t80,'MaxDegree',3);
t80_FASE = matlabFunction(t80_FASE)
t80_CG = solve(eq_CG,t80,'MaxDegree',3);
t80_CG = matlabFunction(t80_CG);
end
E_FASE1=E_FASE;
S_FASE1=S_FASE;
thetai_FASE1=thetai_FASE;
tmax_FASE1=tmax_FASE;
w_FASE1=w_FASE;
xd_FASE1=xd_FASE;
E_CG1=E_CG;
S_CG1=S_CG;
thetai_CG1=thetai_CG;
tmax_CG1=tmax_CG;
w_CG1=w_CG;
xd_CG1=xd_CG;
for i=1:t
Leq_FASE1=Leq_FASE(i);
m_FASE1=m_FASE(i);
thetak_FASE1=thetak_FASE(i);
Leq_CG1=Leq_CG(i);
m_CG1=m_CG(i);
thetak_CG1=thetak_CG(i);
sol(1:3) = t80_FASE(E_FASE1,Leq_FASE1,S_FASE1,m_FASE1,thetai_FASE1,thetak_FASE1,tmax_FASE1,w_FASE1,xd_FASE1);
idx_FASE=sol>0;
ti_FASE(i)=real(sol(idx_FASE))
T_FASE(i)=double(ti_FASE(i)*S_FASE);
p_FASE(i)=double(T_FASE(i)/w_FASE);
sol(1:3) = t80_CG(E_CG1,Leq_CG1,S_CG1,m_CG1,thetai_CG1,thetak_CG1,tmax_CG1,w_CG1,xd_CG1);
idx_CG=sol>0;
ti_CG(i)=real(sol(idx_CG))
T_CG(i)=double(ti_CG(i)*S_FASE);
p_CG(i)=double(T_CG(i)/w_FASE);
end
This code worked when running on appdesigner, but then when I compilled it didn't work. (I also tried without if ~isdeployed(), same result)
Can you please help ?
if ~isdeployed()
syms t80
syms thetai_FASE1 xd_FASE1 E_FASE1 w_FASE1 Leq_FASE1 S_FASE1 thetak_FASE1 tmax_FASE1 m_FASE1
syms thetai_CG1 xd_CG1 E_CG1 w_CG1 Leq_CG1 S_CG1 thetak_CG1 tmax_CG1 m_CG1
eq_FASE= thetai_FASE1+(t80/(xd_FASE1*E_FASE1))-((w_FASE1^2*Leq_FASE1^2)/(24*xd_FASE1*S_FASE1^2*t80^2)) == thetak_FASE1 +(tmax_FASE1/(xd_FASE1*E_FASE1))-((m_FASE1^2*w_FASE1^2*Leq_FASE1^2)/(24*xd_FASE1*S_FASE1^2*tmax_FASE1^2));
eq_CG= thetai_CG1+(t80/(xd_CG1*E_CG1))-((w_CG1^2*Leq_CG1^2)/(24*xd_CG1*S_CG1^2*t80^2)) == thetak_CG1 +(tmax_CG1/(xd_CG1*E_CG1))-((m_CG1^2*w_CG1^2*Leq_CG1^2)/(24*xd_CG1*S_CG1^2*tmax_CG1^2));
t80_FASE = solve(eq_FASE,t80,'MaxDegree',3);
t80_FASE = matlabFunction(t80_FASE)
t80_CG = solve(eq_CG,t80,'MaxDegree',3);
t80_CG = matlabFunction(t80_CG);
end
E_FASE1=E_FASE;
S_FASE1=S_FASE;
thetai_FASE1=thetai_FASE;
tmax_FASE1=tmax_FASE;
w_FASE1=w_FASE;
xd_FASE1=xd_FASE;
E_CG1=E_CG;
S_CG1=S_CG;
thetai_CG1=thetai_CG;
tmax_CG1=tmax_CG;
w_CG1=w_CG;
xd_CG1=xd_CG;
for i=1:t
Leq_FASE1=Leq_FASE(i);
m_FASE1=m_FASE(i);
thetak_FASE1=thetak_FASE(i);
Leq_CG1=Leq_CG(i);
m_CG1=m_CG(i);
thetak_CG1=thetak_CG(i);
sol = t80_FASE(E_FASE1,Leq_FASE1,S_FASE1,m_FASE1,thetai_FASE1,thetak_FASE1,tmax_FASE1,w_FASE1,xd_FASE1);
sol = sol(abs(imag(sol))<1e-6 & real(sol)>0)
if numel(sol) ~= 1
display('Error')
end
ti_FASE(i)=sol;
T_FASE(i)=ti_FASE(i)*S_FASE1);
p_FASE(i)=T_FASE(i)/w_FASE1;
sol = t80_CG(E_CG1,Leq_CG1,S_CG1,m_CG1,thetai_CG1,thetak_CG1,tmax_CG1,w_CG1,xd_CG1);
sol = sol(abs(imag(sol))<1e-6 & real(sol)>0)
if numel(sol) ~= 1
display('Error')
end
ti_CG(i)=sol;
T_CG(i)=ti_CG(i)*S_CG1;
p_CG(i)=T_CG(i)/w_CG1;
end
João Barbosa
2022-12-22
After the aplication compilation this message still appears and it doesn't work.
"syms" are excluded from packaging for the MATLAB Runtime environment according to the MATLAB Compiler license.
Walter Roberson
2022-12-22
When you use matlabFunction in this situation you need to specify 'file' (and 'optimize', 'false') and the else clause needs to create the function handle from the stored file. Nothing after the ~isdeployed test should refer to any symbolic expression or variable.
João Barbosa
2022-12-23
Thank you guys for the help. I finally understood . Sorry for time wasted.
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Assumptions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
