Errors in the MATLAB function added to simulink

2 次查看（过去 30 天）
Awani Khodkumbhe
Awani Khodkumbhe 2019-9-4
评论： Awani Khodkumbhe 2019-9-5
Following is my code for the MATLAB function block in the Simulink.
function gf = fcn(u)
syms x;
y = 0.13*x^4 - 0.9308*x^3 + 2.2267*x^2 + 0.2751*x + 0.0503;
g = finverse(y);
f = 2.19*u;
gf = compose(g,f);
end
I am getting a lot of errors as follows and I do not know how to correct them. I have read the available documentation on the MATLAB function.
Function 'syms' not supported for code generation.
Undefined function or variable 'x'.
Undefined function or variable 'y'. The first assignment to a local variable determines its class.
Undefined function or variable 'g'. The first assignment to a local variable determines its class.
Please help.
Thanks!

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

回答（1 个）

Walter Roberson
Walter Roberson 2019-9-4
编辑：Walter Roberson 2019-9-4
https://www.mathworks.com/help/simulink/ug/calling-matlab-functions.html
Simulink operates in several different modes. In one of the modes, optimization is turned off and everything is simulated at the MATLAB level. That one mode can use the Symbolic Toolbox operations without change if I recall correctly. But it is not the default mode.
The default mode does some optimization producing internal data structures that are faster to simulate but with the MATLAB runtime hanging around in the background. That mode requires that all symbolic toolbox functions be explicitly coded using coder.extrinsic(). You also need to replace the syms x call with x=sym('x') because this mode does not support assignin() that syms uses to create the variable.
Simulink also can operate in a mode where it generates C code and compiles it for faster execution. In this mode, the matlab execution engine is not present and you cannot use the symbolic toolbox at all.
Finally Simulink can also generate code for use on embedded systems or FPGA or raspberry pi or Arduino; the symbolic toolbox is unavailable and most of the targets do not even have any operating system.
So what should you do? Well, you should recode so that you do not need to use the symbolic toolbox in Simulink.
Take your existing fcn and pass it sym('u') at the MATLAB level. You will get out a symbolic expression that can replace the symbolic work you do there. You can then use matlabFunction() with the 'file' option to generate matlab code. Take the generated code and touch it up to do preallocation that Simulink wants and use the resulting pure matlab code instead of the symbolic toolbox calls
  3 个评论
显示 1更早的评论
Walter Roberson
Walter Roberson 2019-9-4
If you do
syms u;
syms x;
y = 0.13*x^4 - 0.9308*x^3 + 2.2267*x^2 + 0.2751*x + 0.0503;
g = finverse(y);
f = 2.19*u;
gf = compose(g,f);
simplify(gf)
then you will get
(39^(1/2)*(272498*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3)*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/2) - 14512125325*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/2) - 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3)*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/2) + 8541000000*u*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/2) + 171145008*39^(1/2)*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/2))^(1/2))/(3900*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/6)*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/4)) - (39^(1/2)*(136249*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/3) - 8541000000*u + 25*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(2/3) + 14512125325)^(1/2))/(3900*(18*13^(1/2)*(9467117712000000000000*u^3 - 46405209395666241120000*u^2 + 63689131782957633368800*u - 1204782748056965920359)^(1/2) - 2792886501600*u + 13803164115967)^(1/6)) + 179/100
If you were to try to use that directly, then you have the problem that it is not vectorized. But
>> matlabFunction(ans)
ans =
function_handle with value:
@(u)sqrt(3.9e+1).*1.0./(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./6.0).*sqrt(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).*(-2.564102564102564e-4)+(sqrt(3.9e+1).*1.0./(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./6.0).*1.0./(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).^(1.0./4.0).*sqrt(sqrt(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).*-1.4512125325e+10+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*sqrt(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).*2.72498e+5-(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*sqrt(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).*2.5e+1+u.*sqrt(u.*-8.541e+9+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(1.0./3.0).*1.36249e+5+(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).^(2.0./3.0).*2.5e+1+1.4512125325e+10).*8.541e+9+sqrt(3.9e+1).*sqrt(u.*-2.7928865016e+12+sqrt(1.3e+1).*sqrt(u.*6.368913178295763e+22-u.^2.*4.640520939566624e+22+u.^3.*9.467117712e+21-1.204782748056966e+21).*1.8e+1+1.3803164115967e+13).*1.71145008e+8))./3.9e+3+1.79e+2./1.0e+2
gives you something that is vectorized. You can replace your fcn with the body of that expression (you might have to break it up into smaller pieces and do some pre-allocation, possibly)
Awani Khodkumbhe
Awani Khodkumbhe 2019-9-5
Thanks a lot for your reply, Walter.
I was able to solve the issue by another way: converting the output format from sym to double.I wrote the MATLAB function and then used the Interprete MATLAB Fcn block instead of the MATLAB Fcn block. It worked with this method.
Thanks again for taking the time and helping me with the MATLAB errors!

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Simulink Functions 的更多信息

标签

产品


版本

R2019a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by