Subscript indices must either be real positive integers or logicals
3 ビュー (過去 30 日間)
古いコメントを表示
Hei All,
I have been getting the error message "Subscript indices must either be real positive integers or logicals. " referring to the calculation of z(i, j) whilst working on the code below. I have spent hours trying to fix the issue, but couldn't figure out what I was missing.
a = 0.2302;
N=4;
theta = linspace(0, 2 * pi * pi, N);
r = linspace(a, 2 * a, N);
sigma_TT = r.^4 .* cos(2 .* theta);
x = zeros(N, N);
y = zeros(N, N);
z = zeros(N, N);
for i = 1:N
for j = 1:N
x(i, j) = (cos(theta(i)))' * r(j);
y(i, j) = (sin(theta(i)))' * r(j);
z(i, j) = real(sigma_TT([theta(i), r(j)]));
end
end
The for loop part is actually what I trying to convert from python to MATLAB. In python it is expressed as:
x = r * np.cos(theta[:,None])
y = r * np.sin(theta[:,None])
z = sigma_TT(SHmax, Shmin, Pp, Pm, a, r, theta[:,None])
Any help on this issue would be much appreciated.
4 件のコメント
Ammar Dodin
2016 年 8 月 16 日
編集済み: Ammar Dodin
2016 年 8 月 16 日
I am not too familiar with numpy, but I think I understand what the first two lines of your python code are doing. It looks like you're creating new vectors that contain the cosine/sin of each of the values of the vector theta after multiplying them with some constant. I don't understand what z is doing? I am assuming sigma_TT is some function you wrote.
Could you possibly provide the Python code for that function?
採用された回答
その他の回答 (1 件)
Image Analyst
2016 年 8 月 16 日
It's pretty well explained in the FAQ: http://matlab.wikia.com/wiki/FAQ#How_do_I_fix_the_error_.22Subscript_indices_must_either_be_real_positive_integers_or_logicals..22.3F
4 件のコメント
Image Analyst
2016 年 8 月 17 日
I thought my comment would explain it more thoroughly than his original 2-line comment. Anyway, hopefully you looked over the FAQ and maybe found some other good stuff in there.
参考
カテゴリ
Help Center および File Exchange で Call Python from MATLAB についてさらに検索
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!