左辺のサイズが 1x1 で右辺のサイズが 0x1 であるため、代入を実行できません。

28 次查看（过去 30 天）
hiroaki miyauchi
hiroaki miyauchi 2019-11-22
评论： hiroaki miyauchi 2019-11-26
左辺と右辺のサイズが異なるため、計算数値を行列に代入できない。
サイズを合わせようとreshape関数を使ったが、うまくいきません。
また、右辺のサイズが 0x1 であるというエラーメッセージの意味がよくわかりません。
よろしくお願いします。
N = 3;
syms r
x = zeros(N+1,N+1);
for i = 0:N
for j = 0:N
f = besselj(i,r);
x(i+1,j+1) = vpasolve(f,r,[j*pi (j+1)*pi]);
end
end
disp(x)
>> eignvalue
左辺のサイズが 1x1 で右辺のサイズが 0x1 であるため、代入を実行できません。
エラー: eignvalue (line 7)
x(i+1,j+1) = vpasolve(f,r,[j*pi (j+1)*pi]);

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

采纳的回答

Hiroki Okawa
Hiroki Okawa 2019-11-23
ツールボックスを持っていないため未検証ですが、エラーメッセージによるとvpasolve関数で解が得られていないようです。
  1 个评论
hiroaki miyauchi
hiroaki miyauchi 2019-11-26
おかげさまで原因が分かりました。ご指摘ありがとうございます。

请先登录，再进行评论。

更多回答（1 个）

Shota Kato
Shota Kato 2019-11-25
右辺のサイズが 0x1 であるというエラーメッセージの意味は，vpasolveによって解が求まっていない，ということです．
vpasolveのドキュメント (https://jp.mathworks.com/help/symbolic/vpasolve.html) が役に立つのではないのでしょうか．
vpasolveの返す値は空のオブジェクトなので，行列の要素として代入することはできません．
これでいいかはわかりませんが，解が求まらないときはx(i+1, j+1)に他の値を代入することでプログラムは回るようになります．
N = 3;
syms r
x = zeros(N+1,N+1);
for i = 0:N
for j = 0:N
f = besselj(i,r);
x_ = vpasolve(f,r,[j*pi (j+1)*pi]);
if isempty(x_)
x(i+1, j+1) = nan;
else
x(i+1,j+1) = x_;
end
end
end
disp(x)
>> eignvalue
2.4048 5.5201 8.6537 11.7915
0 3.8317 7.0156 10.1735
0 5.1356 8.4172 11.6198
0 NaN 6.3802 9.7610
  1 个评论
hiroaki miyauchi
hiroaki miyauchi 2019-11-26
方程式の解を求めることについては お教えしていただいた方法で求められました。
vpasolve 関数は空のプロジェクトを返すことについては初めて知りました。
ご指摘ありがとうございます。

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Symbolic Math Toolbox 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!