Problem running symbolic toolbox in a loop

Question

Prithvi Raj Murthy 2022-9-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1793005-problem-running-symbolic-toolbox-in-a-loop

编辑： Walter Roberson 2023-10-5

hi,i am working on a model identification techinque, so i have used symbolic tool box for some mathematics. I am facing issues when i run the mathematics in a for loop. The program is taking too long to execute and sometimes it even hangs up. I would like to know if it is the symbolic tool box is so heavy on the compuation side if i run it in a loop.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Prithvi Raj Murthy 2022-9-1

编辑：Walter Roberson 2023-10-5

在 MATLAB Online 中打开

syms a b c d f g h j
p_z=[a;b;c;d;f;g;h;j];
theta(1) = (2*f*(h^2+j^2)+2*h*(f^2+g^2))/((f^2)*(h^2)+(f^2)*(j^2)+(g^2)*(h^2)+(g^2)*(j^2));
theta(2) =  (h^2+j^2+f^2+g^2+4*f*h)/((f^2)*(h^2)+(f^2)*(j^2)+(g^2)*(h^2)+(g^2)*(j^2));
theta(3) = (2*h+2*f)/((f^2)*(h^2)+(f^2)*(j^2)+(g^2)*(h^2)+(g^2)*(j^2));
theta(4) = 1/((f^2)*(h^2)+(f^2)*(j^2)+(g^2)*(h^2)+(g^2)*(j^2));
theta(5) = ((a^2)*(c^2)+(a^2)*(d^2)+(b^2)*(c^2)+(b^2)*(d^2));
theta(6) = (2*a*(c^2+d^2)+2*c*(a^2+b^2));
theta(7) = (c^2+d^2+a^2+b^2+4*a*c);
theta(8) = (2*c+2*a);
theta(9) = 1;
%% Simulation error method
for k = 1:N
    if k == 1
        Ysim(k,1) = theta(5).*U(k)
    elseif k == 2
        Ysim(k,1) = theta(1).*Ysim(k-1) + theta(5).*U(k) + theta(6).*U(k-1)
    elseif k == 3
        Ysim(k,1) = theta(1).*Ysim(k-1) + theta(2).*Ysim(k-2) + theta(5).*U(k) + theta(6).*U(k-1) + theta(7).*U(k-2)
    elseif k == 4
        Ysim(k,1) = theta(1).*Ysim(k-1) + theta(2).*Ysim(k-2) + theta(3).*Ysim(k-3) ...
            + theta(5).*U(k) + theta(6).*U(k-1) + theta(7).*U(k-2) + theta(8).*U(k-3)
    else
        Ysim(k,1) = theta(1).*Ysim(k-1) + theta(2).*Ysim(k-2) + theta(3).*Ysim(k-3) + theta(4).*Ysim(k-4)...
            + theta(5).*U(k) + theta(6).*U(k-1) + theta(7).*U(k-2) + theta(8).*U(k-3) + theta(9).*U(k-4)
    end
end
Unrecognized function or variable 'N'.

This is the part of the code where i have used symbolic toolbox and i am using it in the loop

Torsten 2022-9-1

Imagine how big the expressions for Ysim become for N=9 - this will make calculating Ysim symbolically extremely slow.

In the code above, I don't see any reason to use the symbolic toolbox. Try to avoid it for N large.

Prithvi Raj Murthy 2022-9-1

Actually i am using a transfer function model to predict my system from data. so i used symbolic toolbox for the variables to calculate the different parameters of the model. Is there anyother alternate that i can work on

请先登录，再进行评论。

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

Answer 1

Vidip Jain 2023-10-5

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1793005-problem-running-symbolic-toolbox-in-a-loop#answer_1326439

I understand you are using symbolic toolbox for some mathematical computation in a for loop and the program is taking too long to execute.

Yes, the Symbolic Toolbox in MATLAB can be heavy on the computation side, especially when running in a for loop. This is because the Symbolic Toolbox has to perform a lot of symbolic operations, such as differentiation, integration, and simplification.

If you are using the Symbolic Toolbox to calculate the parameters of a transfer function model, then there are a few things you can do to improve the performance of your code:

Avoid using the Symbolic Toolbox in a for loop. If possible, try to calculate the parameters of the transfer function model outside of the for loop. For example, you could calculate the parameters in a separate function and then pass them to the for loop.
Use the Symbolic Toolbox function ‘simplify()’ to simplify the symbolic expressions. This can help to reduce the size of the symbolic expressions and improve the performance of the Symbolic Toolbox.
Use the Symbolic Toolbox function ‘assume()’ to specify the assumptions about the symbolic variables. This can help to improve the performance of the Symbolic Toolbox by reducing the number of symbolic operations that need to be performed.

For further information, refer to the documentation links below: