Error using sym/subsindex
이전 댓글 표시
Hello!
I got the following mistake:
Error using sym/subsindex (line 836)
Invalid indexing or function definition. Indexing must follow MATLAB indexing. Function arguments must be symbolic variables, and function body must be sym
expression.
Error in examp (line 22)
Q(i,j) = (1/(a(i)^(1/2)))* ((t_max - t_min)./(size(t)-1))*symsum(h(k)*wavelet((tau(k)-ttt(j))/a(i)),k, 1, size(h)-1);
The code is following:
type Ez_dry.csv;
Test = importdata('Ez_dry.csv');
t = Test(:, 1);
h = Test(:, 2);
t_max = max(t);
t_min = min(t);
wav = @(tau)(2/(3.^(1/2)*pi^(1/4))*exp(-tau.^2/2).*(1 - tau.^2));
ttt = t_min:0.1*10^-10:t_max;
a = 0.3*10^-10:0.6*10^-11:10^-9;
for i = 1:length(a)
for j = 1:length(ttt)
syms k
Q(i,j) = (1/(a(i)^(1/2)))*((t_max - t_min)./(size(t)-1))*symsum(h(k)*wav((tau(k)-ttt(j))/a(i)),k, 1, size(h)-1);
end
end
S = real(Q);
figure(2);
contour (ttt, 1./a, S, 500)
Could you, please, help me to find a solution. Thank you in advance!
댓글 수: 2
Stijn Haenen
2020년 8월 7일
편집: Stijn Haenen
2020년 8월 7일
There is no symbolic variable as input for symsum. You do h(k), but what does this mean, k had not a value yet so you cannot take the k th element of h. And you didn't define tau so tau(k) does not exist. And you should use length(h) instead of size, i think.
If you have a function of h expressed in k you can solve your problem, i dont know if you can without a function for h.
this line works, but h is excluded from it:
sym_sum_ans=(1/(a(i)^(1/2)))*((t_max - t_min)./(size(t)-1))*symsum(wav((k-ttt(j))/a(i)), 1, length(h)-1);
num_ans=vpa(sym_sum_ans)
Viktoriia Buliuk
2020년 8월 11일
채택된 답변
추가 답변 (1개)
KSSV
2020년 8월 7일
Initialize the cell first.
Q =cell(length(a),length(ttt)) ; % initialize
USe this in the loop.
Q{i,j} = %% your formula with syms
댓글 수: 3
Viktoriia Buliuk
2020년 8월 7일
KSSV
2020년 8월 7일
Check the spelling.....it should be length
Viktoriia Buliuk
2020년 8월 7일
편집: Viktoriia Buliuk
2020년 8월 7일
카테고리
도움말 센터 및 File Exchange에서 Calculus에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!