質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

Q&A

解決済

1回答

381閲覧

解をうまく出力できない

Shuu391

総合スコア1

MATLAB

MATLABはMathWorksで開発された数値計算や数値の視覚化のための高水準の対話型プログラミング環境です。

0グッド

0クリップ

投稿2022/04/23 01:20

matlabの4元連立一次方程式です。解を分数で表示できるものはそのまま小数にせず分数で出力できるようにしたいです。あと解がうまく出力できないです。初歩的で申し訳ないです。

%% % disp('4本の4元1次連立方程式して解きます。');

% /* 4×4の係数行列を入力してください。: */
a11=input('a11 = ');
a12=input('a12 = ');
a13=input('a13 = ');
a14=input('a14 = ');
a21=input('a21 = ');
a22=input('a22 = ');
a23=input('a23 = ');
a24=input('a24 = ');
a31=input('a31 = ');
a32=input('a32 = ');
a33=input('a33 = ');
a34=input('a34 = ');
a41=input('a41 = ');
a42=input('a42 = ');
a43=input('a43 = ');
a44=input('a44 = ');

% /*4×1の定数ベクトルBを入力してください。 */
b1=input('b1 = ');
b2=input('b2 = ');
b3=input('b3 = ');
b4=input('b4 = ');

%/クラメルで解く/
w=[a11 a12 a13 a14;a21 a22 a23 a24;a31 a32 a33 a34;a41 a42 a43 a44];
W=det(w);
if W ~=0
A=[b1 b2 b3 b4;a21 a22 a23 a24;a31 a32 a33 a34;a41 a42 a43 a44];
X=det(A);
x1 =X/W;

B=[a11 a12 a13 a14;b1 b2 b3 b4;a31 a32 a33 a34;a41 a42 a43 a44];
X2=det(B);
x2=(X2/W);

C=[a11 a12 a13 a14;a21 a22 a23 a24;b1 b2 b3 b4;a41 a42 a43 a44];
X3=det(C);
x3=X3/W;

D=[a11 a12 a13 a14;a21 a22 a23 a24;a31 a32 a33 a34;b1 b2 b3 b4];
X4=det(D);
x4=(X4/W);

fprintf ('連立方程式の解\n x1=%d\n x2=%d\n x3=%d\n x4=%d\n', x1,x2,x3,x4);

elseif W == 0
fprintf("解なし\n");

end

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

【解が上手く求まらない理由】
Aの各列をbで置きかえるべきところ,Aの各行をbの転置で置きかえているのが誤りの原因と思います.
【解を分数で求める方法】
sym を使えば良いのではないでしょうか.

MATLAB

1% g0festijer8q1j_revised.m 2 3clc, clear; 4 5n = 4; 6A = zeros(n,n); 7b = zeros(n,1); 8 9detAi = zeros(n,1); 10x_cramer = sym(zeros(n,1)); 11 12for i = 1 : n 13 for j = 1 : n 14 % A(i,j) = input(sprintf("A(%d,%d): ",i,j)); 15 A(i,j) = randi(19) - 10; 16 end 17end 18for i = 1 : n 19 % b(i) = input(sprintf("b(%d): ",i)); 20 b(i) = randi(19) - 10; 21end 22 23disp('A = '); 24disp(A); 25disp('b = '); 26disp(b); 27x = sym(A) \ sym(b); 28disp('A \ b ='); 29disp(x); 30 31detA = sym(det(A)); 32fprintf("det(A) = "); 33disp(detA); 34for i = 1 : n 35 Ai = A; 36 Ai(:,i) = b; 37 % fprintf("A%d = \n",i); 38 % disp(Ai); 39 detAi = sym(det(Ai)); 40 fprintf("det(A%d) = ",i); 41 disp(detAi); 42 x_cramer(i) = detAi / detA; 43end 44disp("Cramer's rule: "); 45disp(x_cramer);

投稿2022/05/05 08:57

obeSolitary

総合スコア106

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問