Subs does not work in my code.

Question

Valeria Ramirez Ariiola 2024-4-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2112741-subs-does-not-work-in-my-code

评论： Dyuman Joshi 2024-5-28

Hi, my function should replace all the monomials that have negative degree from -M to -(4*M+8) with 1, since 0 will give an error, nevertheless, it doesn't work the substitution function. Every time I try it, it doesn't work. Could you please give me some ideas on how to solve this issue.

function  [sMnew] = Cubic(k,M)
n= 3*k+1;
x= sym('x');
c = sym( 'c', [1 M+3]);
% Define the sum of monomials with negative coefficients
CC = x;
for i = 2:M+3
    CC = CC + c(1,i)*x^(-i+1);
end
sM = expand(CC^(n));
for i = M:1:4*M+8    
    sMnew = subs(sM, x^(-i), 1);
end

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

请先登录，再进行评论。

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

Answer 1

Dyuman Joshi 2024-4-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2112741-subs-does-not-work-in-my-code#answer_1449436

在 MATLAB Online 中打开

You can vectorize operations in the code -

function [sMnew] = Cubic(k,M)
n = 3*k+1;
x = sym('x');
c = sym( 'c', [1 M+3]);
%Indices
k = 2:M+3;
%Vectorized sum
CC = x + sum(c(1,k).*x.^(-k+1));
%Define powers for which to substitute
vec = -1*(M:4*M+8);
%Substitute directly at once
sMnew = subs(z, x.^(vec), ones(size(vec)));
end