Implementation of infinite series in MATLAB.
조회 수: 31 (최근 30일)
이전 댓글 표시
How I can impliment this type of infinite series in MATLAB given a function as
where
, here l,\gamma are constant. when I use
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/516292/image.png)
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/516297/image.png)
syms k;
exp = symsum(H(k + 1,gamma_av),k,1,Inf);
I used to get error like this "The following error occurred converting from sym to double:
Unable to convert expression into double array."
Please help me out.
댓글 수: 5
Walter Roberson
2021년 2월 11일
The error is not in the part of the code that you posted. We need to see the rest of the code.
채택된 답변
Walter Roberson
2021년 2월 13일
The error message you get is because the symbolic engine was not able to find a convergent value for the infinite series when you use double(). And depending on the number of digits you use for vpa() it might or might not find a solution.
Interestingly, if you use a finite series such as 75 terms, then double() does well at converting the value, whereas you need on the order of 150+ digits in order for vpa to be able to resolve it: if you use a high finite Limit and too few digits then vpa() will tend to return exactly 0.
format long g
syms k;
snr = 2;
l = 2;
Limit = Inf;
expression2 = symsum(H(k + 1,snr,l),k,1,Limit);
for D = 120:129
try
D
ev = vpa(expression2, D)
ed = double(ev)
catch ME
fprintf('failed working at %d digits\n', D);
end
end
double(expression2)
function [out] = H(y,z,l)
Pi = sym(pi);
%gam, y,z,l are user inputs
m = -1/log2(cos(Pi/6)); % in paper its used as small gamma.(\gamma)
out1 = (2/z)^y;
out2 = gamma(y) - igamma(y,(l^(-2*(m + 2))/2)*z);
out = out1 * out2;
end
댓글 수: 0
추가 답변 (1개)
Sahil
2022년 9월 18일
편집: Walter Roberson
2022년 9월 18일
function [out] = H(y,z,l)
Pi = sym(pi);
%gam, y,z,l are user inputs
m = -1/log2(cos(Pi/6));
% in paper its used as small gamma.(\gamma)
out1 = (2/z)^y;
out2 = gamma(y) - igamma(y,(l^(-2*(m + 2))/2)*z);
out = out1 * out2;
end
댓글 수: 1
Walter Roberson
2022년 9월 18일
What is the difference between that and what I posted last year at https://www.mathworks.com/matlabcentral/answers/742277-implementation-of-infinite-series-in-matlab#answer_622427 ?
참고 항목
카테고리
Help Center 및 File Exchange에서 Calculus에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!