行列の2乗と3乗を求める方法とは？

前へ 次へ
このQ&Aのポイント
  • 行列Aの2乗と3乗を求める方法について教えてください。
  • 与えられたプログラムでA_2とA_3の計算方法を教えてください。
  • セグメントエラーが出てしまう原因についても教えていただけますか？
回答を見る
  • ベストアンサー

3行3列の行列 A に対して、Aの2乗, Aの3乗 を求めよ。

3行3列の行列 A に対して、Aの2乗, Aの3乗 を求めよ。 という問題で #include <stdio.h> #define N 3 int A[N][N]; int A_2[N][N]; /* 行列Aを2乗したもの */ int A_3[N][N]; /* 行列Aを3乗したもの */ int main() { int i,j,k,l; /* 3行3列の係数行列Aの成分を入力 */ printf("係数行列Aを%d行%d列で入力してください\n", N, N); for( i=0; i<N; i++) { for( j=0; j<N; j++) { printf("A[%d][%d]=", i+1, j+1); scanf("%d", &A[i][j]); } } ここまでがあたえられえているのですが、単純にA_2＝A×Aとしてみたのですが、セグメントエラーが出てしまいうまく出力できませんでした。 A_2とA_3の計算方法とあたえられているプログラムの動作を教えてください。 よろしくお願いします

質問者が選んだベストアンサー

  • ベストアンサー
  • asuncion
  • ベストアンサー率33% (2126/6286)
回答No.3

| A B C | | a b c | | Aa+Bd+Gc Ab+Be+Ch Ac+Bf+Ci | | D E F | | d e f | | Da+Ed+Fc Db+Ee+Fh Dc+Ef+Fi | | G H I | | g h i | | Ga+Hd+Ic Gb+He+Ih Gc+Hf+Ii | 本当ですか？ 積のところで、gがどこにも登場していません。 まあ、それはさておき、 X = A * B のとき、 積Xのi行j列要素は、Aのi行目とBのj行目の積和です。 いきなり3行3列の積を求めようとするより、 まずは2行2列あたりでじゅうぶん考え方をマスターなさるのが よいのではないでしょうか。

dyebe011
質問者

お礼

返信が遅れてしまい申し訳ありません。 　 ＞本当ですか？ ＞積のところで、gがどこにも登場していません。 間違った1行目をコピーしたためｇが出ていませんでした。 2行2列から勉強しようと思います。 こんな私に付き合っていただき、ありがとうございました。

その他の回答 (3)

  • asuncion
  • ベストアンサー率33% (2126/6286)
回答No.4

＞積Xのi行j列要素は、Aのi行目とBのj行目の積和です。 おっと失礼しました。 誤：Aのi行目とBのj行目の積和 正：Aのi行目とBのj列目の積和

  • asuncion
  • ベストアンサー率33% (2126/6286)
回答No.2

＞A_2[i][j] = A[i][j] * A[i][j]; 間違っています。 行列の掛け算の方法を勉強し直してください。

dyebe011
質問者

補足

ご指摘ありがとうございます。 例として、3行3列の行列の積は | A B C | | a b c | | Aa+Bd+Gc Ab+Be+Ch Ac+Bf+Ci | | D E F | | d e f | | Da+Ed+Fc Db+Ee+Fh Dc+Ef+Fi | | G H I | | g h i | | Ga+Hd+Ic Gb+He+Ih Gc+Hf+Ii | となるので、このプログラムでその計算をさせるには for(i=0;i<3;i++) for(j=0;j<3;j++) A_2[i][J] += A[I][j]*A[j][i]; と計算すればよいのでしょうか？

  • asuncion
  • ベストアンサー率33% (2126/6286)
回答No.1

＞単純にA_2＝A×Aとしてみた 考え方は正しいです。 うまくいかない、というコードを貼ってください。

dyebe011
質問者

補足

回答ありがとうございます。 A_2[i][j] = A[i][j] * A[i][j]; としました。 あたえられているプログラムだと int i,j,k,l; となっていますが、k,lはどのようにつかうのでしょうか？

  1. カテゴリ
  2. [技術者向] コンピューター
  3. プログラミング・開発
  4. C・C++・C#

関連するQ&A

  • 行列をべき乗させるプログラム

    2行2列を５乗させるプログラムを作って、一応できたつもりだったんですが結果が合いません・・・ 何かヒントでもいいのでわかる方いらっしゃいましたらよろしくお願いします。 <プログラム> #include <stdio.h> #define N 2 int A[N][N]; int A_NEW[N][N]; int A_5[N][N]; /* 行列Aを５乗したもの */ int main() { int i,j,k,l; /* 2行2列の係数行列Aの成分を入力 */ printf("係数行列Aを%d行%d列で入力してください\n", N, N); for( i=0; i<N; i++) { for( j=0; j<N; j++) { printf("A[%d][%d]=", i+1, j+1); scanf("%d", &A[i][j]); } } for(i=0; i<N; i++) /* A_NEW=A*Aの計算 */ { for(j=0; j<N; j++) { for(k=0; k<N; k++) { A_NEW[i][j] += A[i][k] * A[k][j]; } } } for(l=0; l<3; l++) { for(i=0; i<N; i++) /* A_5の計算 */ { for(j=0; j<N; j++) { for(k=0; k<N; k++) { A_5[i][j] += A_NEW[i][k] * A[k][j]; } } } for(i=0; i<N; i++) { for(j=0; j<N; j++) { A_NEW[i][j] = A_5[i][j]; } } } printf("A_5=\n"); /* 出力 */ for( i=0; i<N; i++) { for( j=0; j<N; j++) { printf("%d ", A_5[i][j]); } printf("\n"); } } <入力例> A= 1 3 2 1 <期待する結果> A= 241 303 202 241 <このプログラムの結果> 406 498 332 406

  • 行列式のプログラムでのループがわかりません。

    行列式を求めるプログラムを考えていて、一応はちゃんと答えがでるのですが、 ループを使って計算させようとしたのですがイマイチわかりません。 #include <stdio.h> main() { 　　　　int a[3][3]={{1,2,2},{-1,5,3},{2,0,1}}; 　　　　int i,j,n; 　　　　for(i=0;i<3;i++) 　　　　{ 　　　　　　　　　for(j=0;j<3;j++) 　　　　　　　　　　　　　printf("%3d",a[i][j]); 　　　　　　　　　　　　　printf("\n"); 　　　　} 　　　　n=a[0][0]*(a[1][1]*a[2][2]-a[1][2]*a[2][1]) 　　　　　-a[1][0]*(a[0][1]*a[2][2]-a[0][2]*a[2][1]) 　　　　　+a[2][0]*(a[0][1]*a[1][2]-a[0][2]*a[1][1]); 　　　　printf("この行列の行列式は \n"); 　　　　printf("%d \n",n); } 1 2 2 -1 5 3 2 0 1 という行列式です。 nの処理をループを使うとどう書いたらいいでしょうか？ これはたまたま３行３列だからひとつひとつ指定して計算できますが これが４行４列とかだったら絶対ループ使わなきゃ面倒ですからね～。 とりあえず３行３列のが理解できたらいいとは思うのですがなかなか複雑で難しいですね。

  • c言語のプログラムで行列の積を計算する

    指定された行・列数（それぞれ10以下とする）の行列 X, Y の積 Z = X × Y を求めるプログラムを作成せよ．行列の要素はすべて整数とする． このプログラムで行と列の成分を入力するときに、1列になってしまうのですがどうしたらちゃんと入力できますか？ int main(void) { double A[10][10]; double B[10][10]; int i,j,m,n,p,k; printf("行列xの行数は？:"); scanf("%d",&m); printf("行列xの列数(行列yの行数)は？:"); scanf("%d",&n); printf("行列yの列数は？:"); scanf("%d",&p); printf("行列xを入力してください。\n"); for (i = 0; i < m; i++) { for (j = 0; j < n; j++) { scanf("%d", &A[i][j]); } } printf("行列yを入力してください。\n"); for (i = 0; i < n; i++) { for (j = 0; j < p; j++) { scanf("%d", &B[i][j]); } }

  • 行列の行の入換え

    | -1 -1 -2 -1 10 | | -1 -2 -1 10 -1 | | -2 -1 10 -1 -1 | | -1 10 -1 -1 -2 | | 10 -1 -1 -2 -1 | この行列を次の順序で並び替えたいのですが 1.第１列にある要素の中で絶対値の最大値がある行を探し，その行と第1行を入れ換える． 2.第２列にある要素の中で絶対値の最大値がある行を第2行以下から探し，その行と第２行を入れ換える． 3.これを左上から右下へ順次移動して，行の入れ換えを行う． #include<stdio.h> #include<math.h> #define N 5 #define M 5 int main(){ float a[N][M]={{-1.0,-1.0,-2.0,-1.0,10.0},{-1.0,-2.0,-1.0,10.0,-1.0},{-2.0,-1.0,10.0,-1.0,-1.0}, {-1.0,10.0,-1.0,-1.0,-2.0},{10.0,-1.0,-1.0,-2.0,-1.0}}; float nmax; int i,j,k=0,s,imax; for(j=0; j<M; j++){ nmax = a[0][j]; imax = 0; for(i=1; i<N; i++){ if(fabs(nmax) < fabs(a[i][j])){ nmax = a[i][j]; imax = i; } } for(j=0; j<M; j++){ s=a[k][j]; a[k][j]=a[imax][j]; a[imax][j]=s; } k++; } for(i=0; i<N; i++){ for(j=0; j<M; j++){ printf("%5.1f ",a[i][j]); } printf("\n"); } return(0); } このように作成したのですが、1行目と5行目しかいれかわりません。 なぜでしょうか？？

  • 行列の積を関数を使って求める・・？

    ２つの行列の行と列を入力し、積を計算するプログラムを関数を使って書きたいのですが、上手く行きません。どこをどのように直したらよいか教えてください！お願いします！！ 以下が私が書いたプログラムです。 #include<stdio.h> #define NUMBER 10 int first(int x1,int x2,int y1,int y2,int i,int j,int k) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; do{ printf("２つの行列の行と列を入力してください\n"); scanf("%d", &x1); scanf("%d", &x2); scanf("%d", &y1); scanf("%d", &y2); if(x1 != y2){ printf("行列の積は計算できません\n"); } }while(x1 != y2); printf("行列Ａの要素を入力してください\n"); for(i=0; i<x1; i++){ for(j=0; j<x2; j++) scanf("%d", &a[i][j]); } printf("行列Ｂの要素を入力してください\n"); for(j=0; j<y1; j++){ for(k=0; k<y2; k++) scanf("%d", &b[j][k]); } } int second(int x1,int x2,int y1,int y2,int i,int j,int k) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; for(i=0; i<x1; i++){ for(k=0; k<y2; k++){ for(j=0; j<x2; j++) c[i][k] = c[i][k] + a[i][j]*b[j][k]; } } for(i=0; i<x2; i++){ for(k=0; k<y2; k++) printf("%3d", c[i][k]); printf("\n"); } } int main(void) { int a[NUMBER][NUMBER] = {0}; int b[NUMBER][NUMBER] = {0}; int c[NUMBER][NUMBER] = {0}; printf("行列の積を計算します\n %d\n", first(x1,x2,y1,y2,i,j,k)); printf("行列Ａと行Ｂの積は\n %3d",second(x1,x2,y1,y2,i,j,k)); }

  • C言語を用いた45×４５の逆行列の表示について・・・・。。

    次のプログラムを動かすと逆行列を求めることができます。 １０×１０の計算までは、スムーズに計算する事ができたのですが、 私が知りたい４５×４５の行列の計算になるとなぜかうまくいきません。 そこでみなさんにお聞きしたいのは、計算がうまくいかない理由がプログラムにあるのか、単にＰＣの限界なのか、それ以外なのか。みなさんからアドバイスをいただけないでしょうか？ ※下のプログラム、見づらくてすいません。 #include <stdio.h> #define N 10 /* 元の数 */ #define N2 2*N /* プロトタイプ宣言 */ void Print(int k); /* 大域変数 */ double a[N][N2]= {/* 行列 Aおよび単位行列 */ { ここに45×45の正則行列＆45×45の単位行列を入れます・・・。。 }; int main(void) { double d; int i,j,k; for (k=0;k<=N-1;k++) { Print(k); for(i=0;i<=N-1;i++) { if(i!=k) { d=a[i][k]/a[k][k]; for (j=k+1;j<=N2-1;j++) a[i][j]=a[i][j]-a[k][j]*d; a[i][k] = 0.0; } } d=a[k][k]; for (j=k;j<=N2-1;j++) a[k][j]=a[k][j]/d; } Print(N); } void Print(int k) { int i,j; printf("\n ステップ %d \n",k); /* 小見出し */ for(i=0;i<=N-1;i++) { /* 第ｉ行を印刷 */ for(j=0;j<=N-1;j++) printf("%6.4f,",a[i][j]); printf(" : "); for(j=N;j<=N2-1;j++) printf("%6.4f,",a[i][j]); printf("\n"); } }

  • 行列の積について

    3行3列の行列AとB、およびその和と積を表示するものです #include <stdio.h> int main (void) { int i,j,k; double a[3][3] = {{2.4, 5.5, -8.5},{0.8, 3.7, 1.1},{3.5, -9.1, 2.6}}; double b[3][3] = {{-5.1, 9.8, 2.3},{-4.1, 0.2, -0.3},{3.3, 6.1, -1.3}}; double c[3][3] = {0}; printf("　行列A\n"); for (i = 0; i < 3; i++){ for (j = 0; j < 3; j++) printf("%5.1f", a[i][j]); printf("\n"); } printf("\n　行列B\n"); for (i = 0; i < 3; i++){ for (j = 0; j < 3; j++) printf("%5.1f", b[i][j]); printf("\n"); } printf("\n　行列A+B\n"); for (i = 0; i < 3; i++){ for (j = 0; j < 3; j++) c[i][j] = a[i][j] + b[i][j]; } for (i = 0; i < 3; i++) { for (j = 0; j < 3; j++) printf("%5.1f", c[i][j]); printf("\n"); } for (i=0; i<3 ; i++){ for (j=0; j<3; j++) for (k=0; k<3; k++) c[i][j] =c[i][j]+ a[i][k]* b[k][j]; } printf("\n　行列AB\n"); for (i = 0; i < 3; i++) { for (j = 0; j < 3; j++) printf("%7.2f", c[i][j]); printf("\n"); } return 0; } これだと積だけが異なった値が出てしまいます。 試しに、積の部分だけで組んでみると正しい答えが出ました。 for文のブロックの組み方がまずいのかなと思っていますが、これ以外だと、行列のレイアウトが崩れてしまいます どのように直したらいいのか、ご教示お願いします

  • 以前も同じ質問したのですがもう一度お願いします

    以前も同じ質問したのですがもう一度お願いします プログラミングで３行３列の行列Ａに対しＡの２乗、Ａの３乗を求めるプログラムを作成したいのですが下記のプログラムでは３行３列まで表示できません どこが間違えているのか修正していただけないでしょうか？ #include <stdio.h> int main(){ int z[3][3]; int w[3][3]; int x[3][3]={1,2,3,4,5,6,7,8,9}; int n,m; int i,j,k; printf("z[3][3] : \n"); for(i=0;i<3;i++) { for(j=0;j<3;j++) {printf("z[%d][%d] = ",i,j); scanf ("%d",&z[i][j]); } } printf("m for [z]^m ="); scanf ("%d",&m); for(n=0;n<m;n++) { for(i=0;i<3;i++) { for(j=0;j<3;j++) { w[i][j] = 0; for(k=0;k<3;k++) { w[i][j] += x[i][k]*z[k][j]; } } } for(i=0;i<3;i++) { for(j=0;j<3;j++) { x[i][j] = w[i][j]; } } printf("[z]^%d\n",n+1); for(i=0;i<3;i++) { for(j=0;j<3;j++) { printf(" %9d",x[i][j]); }printf("\n"); }printf("\n"); } return 0; }

  • c言語で行列の積の値を求める

    行列の成分を入力した後に、 入力された行列は X = 1 2 3 4 5 6 7 8 9 10 11 12 Y = 1 5 2 6 3 7 4 8 のように表示して、(上の数字は適当です。) 行列 X と行列 Y の積を求めて結果を表示するプログラムが作りたいのですが、上手く表示できなくて困っています。 #include <stdio.h> int main(void) { double A[10][10]; double B[10][10]; double C[10][10]; int i,j,m,n,p,k; printf("行列xの行数は？:"); scanf("%d",&m); printf("行列xの列数(行列yの行数)は？:"); scanf("%d",&n); printf("行列yの列数は？:"); scanf("%d",&p); printf("行列xを入力してください。\n"); for (i = 0; i < m; i++) { for (j = 0; j < n; j++) { scanf("%d", &A[i][j]); } } printf("行列yを入力してください。\n"); for (i = 0; i < n; i++) { for (j = 0; j < p; j++) { scanf("%d", &B[i][j]); } } この後にどうすればいいのか教えてください。 よろしくお願いします。

  • C言語　プログラミング　行列演算

    下記のプログラムのおかしい点と解決法を教えてください。 コンパイルは通りますがうまく動きません。。 #include<stdio.h> #define MAX 500 int main(void){ int matrA[MAX][MAX],matrB[MAX][MAX],matrC[MAX][MAX],l,m,n,i,j,k; printf("lとmを入力してください:"); scanf("%d",&l); scanf("%d",&m); printf("行列Aを入力してください"); for(i=0;i<l;i++){ printf("＞"); for(j=0;l<m;j++){ scanf("%d",&matrA[i][j]); } printf("\n"); } printf("nを入力してください(m = %d):",m); scanf("%d",&n); printf("行列Bを入力してください"); for(i=0;i<m;i++){ printf("＞"); for(j=0;j<n;j++){ scanf("%d",&matrB[i][j]); } printf("\n"); } printf("C=\n"); for(i=0;i<l;i++){ for(j=0;j<n;j++){ for(k=0;k<m;k++){ matrC[i][j]+=matrA[i][k]*matrB[k][j]; } printf("%d",matrC[i][j]); } printf("\n"); } }

注目のQ&A

カテゴリ

専門家に質問してみよう

職業から探して質問する

あなたにピッタリな商品が見つかる！ OKWAVEセレクト

質問する