How to handle mupadmex error?

Question

Niloufar 2022-12-27

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1884917-how-to-handle-mupadmex-error

回答： Walter Roberson 2023-1-2

在 MATLAB Online 中打开

I have a wave equation that is like this:

and here is the initial conditions:

I should create the gif for a variety of c but I get mupadmex error Which I don't know how to fix and here is my code:

clear;clc;close all;
a = 1;
t0 = 1;
alpha = 1;
syms x;
%define lambda function
f = @(x) ((-a/t0)*x-a)*(x>-2*t0-0.1 & x<-t0) + ((a/t0)*x+a)*(x>-t0 & x<=0) + ((-a/t0)*x+a)*(x>0 & x<t0) + ((a/t0)*x-a)*(x>t0 & x<2*t0+0.1);
g = @(x) alpha.*f(x);
for c = 0:0.05:1
    h(x) = int(g(x));
    G = matlabFunction(h);
    desired_fun = G(x+c) - G(x-c);
    my_plot = fplot(((f(x-c) + f(x+c))/2) + desired_fun,[-t0 t0]);
    frame = getframe(1);
    im = frame2im(frame); %pause(0.01);
    [imind,cm] = rgb2ind(im,256);
    imwrite(imind,cm,"C:\Users\ernika\wave_equation.gif",'gif','WriteMode','append'); %saved as mygif1
    drawnow;
    if(c~=1)
        delete(my_plot)
    end
end

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

请先登录，再进行评论。

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

Answer 1

Ajay Gajulapally 2023-1-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1884917-how-to-handle-mupadmex-error#answer_1139447

在 MATLAB Online 中打开

The mupadmex error is caused because of your initial values of a, t0, alpha.

Consider Changing those values which resolves your mupadmex error.

a = 10;
t0 = 100;
alpha = 1000;

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

请先登录，再进行评论。

Answer 2

Walter Roberson 2023-1-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1884917-how-to-handle-mupadmex-error#answer_1139472

在 MATLAB Online 中打开

Do not define the functions that way. The symbolic toolbox does not define a result for multiplying a symbolic expression by a logical value.

You need to define the function symbolically using piecewise. Then integrate it with a definite integral specifying the variable of integration and the upper and lower bound. You can make a bound the variable of integration if you need to.

Then use matlabFunction with the 'file' option and with 'optimize' off. You must use 'file' for this purpose. 'file' option is able to generate code for piecewise.

Now, the function that results is not vectorized, so if you try to use it in a context such as integral() or fplot() then that will fail; you would need an arrayfun wrapper in such cases.

You appear to be passing symbolic x back into the result of matlabFunction, but that is going to fail.

I would suggest that you reconsider using matlabFunction and use

G = h;

which will avoid the issues about code generation and will allow a proper desired_function to be defined symbolically.