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

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

新規登録して質問してみよう
ただいま回答率
86.12%
多次元配列

1次元配列内にさらに配列を格納している配列を、多次元配列と呼びます。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。

解決済

掃き出し法を用いた逆行列の導出

HYYH
HYYH

総合スコア15

多次元配列

1次元配列内にさらに配列を格納している配列を、多次元配列と呼びます。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。

1回答

0グッド

0クリップ

200閲覧

投稿2022/11/21 15:27

編集2022/11/21 16:27

逆行列を求めるプログアムを作成したが、出力結果が

-nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan -nan

となってしまう。
なぜ非数になるのか、また改善策を教えていただけませんか

c

1#include<stdio.h> 2#include<stdlib.h> 3#define N 8 4int main(void){ 5 double a[N][N]={{0,0,1,0,0,0,0,0},{1,0,1,0,0,0,2,0},{1,1,1,0,0,0,-0.7,-0.7},{0,1,1,0,0,0,0,-0.3},{0,0,0,0,0,1,0,0},{0,0,0,1,0,1,0,0},{0,0,0,1,1,1,-1,-1},{0,0,0,0,1,1,0,-1}}; //入力用の配列 6 double inv_a[N][N]; //ここに逆行列が入る 7 double buf; //一時的なデータを蓄える 8 int i,j,k; //カウンタ 9 10 //単位行列を作る 11 for(i=0;i<N;i++){ 12 for(j=0;j<N;j++){ 13 if(i==j){ 14 inv_a[i][j]=1.0; 15 } 16 else{ 17 inv_a[i][j]=0.0; 18 } 19 } 20 } 21 //掃き出し法 22 for(i=0;i<N;i++){ 23 buf=1/a[i][i]; 24 for(j=0;j<N;j++){ 25 a[i][j]*=buf; 26 inv_a[i][j]*=buf; 27 } 28 for(j=0;j<N;j++){ 29 if(i!=j){ 30 buf=a[j][i]; 31 for(k=0;k<N;k++){ 32 a[j][k]-=a[i][k]*buf; 33 inv_a[j][k]-=inv_a[i][k]*buf; 34 } 35 } 36 } 37 } 38 //逆行列を出力 39 for(i=0;i<N;i++){ 40 for(j=0;j<N;j++){ 41 printf(" %10lf",inv_a[i][j]); 42 } 43 printf("\n"); 44 } 45}

以下のような質問にはグッドを送りましょう

  • 質問内容が明確
  • 自分も答えを知りたい
  • 質問者以外のユーザにも役立つ

グッドが多くついた質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

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

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

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

下記のような質問は推奨されていません。

  • 間違っている
  • 質問になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

適切な質問に修正を依頼しましょう。

matukeso

2022/11/21 16:10

なぜ、以前のガウスジョーダン法の質問では正しかった「単位行列を作る」処理が、この質問ではまちがっているの?
HYYH

2022/11/21 16:26

ご指摘ありがとうございます。 1だけ代入して終わってました、訂正いたします。

回答1

0

ベストアンサー

何が間違っているのかというと、行の交換ステップがないところです。
行列が密ならそれでもいいですが、たとえば1行目は0,0,1,0,0,0,0,0なわけで、行交換しないと1/0とかを計算する羽目になります。
たとえば(k,k)から下に(n,k)まで見て、絶対値の最大のものが第 k 行にくるように行の入れ替えを行うとかですね

投稿2022/11/21 17:34

matukeso

総合スコア1425

良いと思った回答にはグッドを送りましょう。
グッドが多くついた回答ほどページの上位に表示されるので、他の人が素晴らしい回答を見つけやすくなります。

下記のような回答は推奨されていません。

  • 間違っている回答
  • 質問の回答になっていない投稿
  • スパムや攻撃的な表現を用いた投稿

このような回答には修正を依頼しましょう。

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

ただいまの回答率
86.12%

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

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

質問する

関連した質問

同じタグがついた質問を見る

多次元配列

1次元配列内にさらに配列を格納している配列を、多次元配列と呼びます。

配列

配列は、各データの要素(値または変数)が連続的に並べられたデータ構造です。各配列は添え字(INDEX)で識別されています。