ページ 1 ／ 1
行列とベクトルの積
Posted: 2017年5月31日(水) 23:45
by あい
ベクトルの積を求めたいんですが、うまくいきません
Void MatProから下を教えてほしいです
形はなるべく変えずにお願いします。
#include<stdio.h>
void MatPro(double(*Mat)[3], double *Vec, double *ReVec);
double innerp(int n, double *a, double *b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 }, { 43, 80, -54 }, { 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat, Vec, ReVec);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", *(ReVec + i), *(*(Mat + i)), *(*(Mat + i) + 1), *(*(Mat + i) + 2), *(Vec + i));
return 0;
}
void MatPro(double(*Mat)[3], double *Vec, double *ReVec)
{
int i,j;
for (i = 0; i < 3; i++){
for (j = 0; j < 0; j++){
Mat[j] * Vec;
}
}
}
double innerp(int n, double *a, double *b)
{
double ip = 0;
int i;
for (i = 0; i < n; i++)
ip += *(a + i)**(b + i);
return ip;
}
Re: 行列とベクトルの積
Posted: 2017年6月01日(木) 00:25
by みけCAT
質問者からの返事がなく、あまりいいトピックではありませんが、参考として過去ログに同様の問題に関する質問があります。
ポインタ • C言語交流フォーラム ～ mixC++ ～
ベクトルの積 • C言語交流フォーラム ～ mixC++ ～
Re: 行列とベクトルの積
Posted: 2017年6月01日(木) 09:04
by Math
コード:
#include<stdio.h>
void MatPro(double(*Mat)[3], double* Vec, double* ReVec);
double innerp(int n, double* a, double* b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
//MatPro(Mat, Vec, ReVec);
MatPro(Mat, Vec, ReVec);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", *(ReVec + i), *(*(Mat + i)), *(*(Mat + i) + 1), *(*(Mat + i) + 2), *(Vec + i));
return 0;
}
void MatPro(double(*Mat)[3], double *Vec, double *ReVec)
{
int i;//, j;
for (i = 0; i < 3; i++) {
//for (j = 0; j < 0; j++) {
//Mat[i][j] * Vec[i];
//}
*(ReVec + i) = innerp(3, (Mat)[i], Vec);
}
}
double innerp(int n, double *a, double *b)//d:n次元:a行Vector:b列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
コード:
1>------ すべてのリビルド開始: プロジェクト:ConsoleApplication0601, 構成: Debug Win32 ------
1>c1.c
1>ConsoleApplication0601.vcxproj -> D:\z17\ｃ\ConsoleApplication0601\Debug\ConsoleApplication0601.exe
1>ConsoleApplication0601.vcxproj -> D:\z17\ｃ\ConsoleApplication0601\Debug\ConsoleApplication0601.pdb (Partial PDB)
========== すべてリビルド: 1 正常終了、0 失敗、0 スキップ ==========
コード:
output=Mat[3][3](input)*Vec[3](input)
|4824.00| | 12.00 34.00 82.00| | 2.00|
|3478.00| | 43.00 80.00 -54.00| | 64.00|
|4592.00| | 56.00 31.00 78.00| | 32.00|
続行するには何かキーを押してください . . .
VS2017を使用してます。
EXCELで検証
注：ポインターと2次元行列は複雑なので疑問点があったら質問をしてくださいね。(^^;
[/size]
Re: 行列とベクトルの積
Posted: 2017年6月02日(金) 23:37
by MathM
書き方が雑だったので清書して補足説明をします。
全コード c1.c
コード:
#include<stdio.h>
void MatPro(double(*Mat)[3], double* Vec, double* ReVec);
double innerp(int n, double* a, double* b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat, Vec, ReVec);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", *(ReVec + i), *(*(Mat + i)), *(*(Mat + i) + 1), *(*(Mat + i) + 2), *(Vec + i));
return 0;
}
void MatPro(double(*Mat)[3], double *Vec, double *ReVec)
{
int i;
for (i = 0; i < 3; i++) {
*(ReVec + i) = innerp(3, (Mat)[i], Vec);
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
バッチファイル c.bat
コード:
rem コンパイル後リンク
cl /TC c1.c
rem 実行結果
c1.exe
VS2017 コマンドプロンプトで実行
コード:
D:\z17c\c\0601>c
D:\z17c\c\0601>rem コンパイル後リンク
D:\z17c\c\0601>cl /TC c1.c
Microsoft(R) C/C++ Optimizing Compiler Version 19.10.25019 for x86
Copyright (C) Microsoft Corporation. All rights reserved.
c1.c
Microsoft (R) Incremental Linker Version 14.10.25019.0
Copyright (C) Microsoft Corporation. All rights reserved.
/out:c1.exe
c1.obj
D:\z17c\c\0601>rem 実行結果
D:\z17c\c\0601>c1.exe
output=Mat[3][3](input)*Vec[3](input)
|4824.00| | 12.00 34.00 82.00| | 2.00|
|3478.00| | 43.00 80.00 -54.00| | 64.00|
|4592.00| | 56.00 31.00 78.00| | 32.00|
D:\z17c\c\0601>
注：ポインターと2次元配列は複雑なので以下の通り補足説明をします。
このコードがC言語の正しい書き方なのですが人間にはわかりずらいので”配列”風に書き直して見ることができます。
【C言語では添字演算子[]は配列とは無関係です。C言語には多次元配列は存在しません。多次元配列のように見えるものは「配列の配列」です。この方が使いやすいこともあります。】
シンタックスシュガーを使って書き換えてみます。
以下参照：
http://dixq.net/forum/viewtopic.php?f=3&t=19137
[/size]
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 00:19
by Math
シンタックスシュガーを使って書き換えたコードです。
コード:
#include<stdio.h>
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
double innerp(int n, double *a, double *b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat[3][3], Vec[3], ReVec[3]);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", ReVec[i], Mat[i][0], Mat[i + 1][1],Mat[i + 2][2], Vec[i];
return 0;
}
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
{
int i;
for (i = 0; i < 3; i++) {
ReVec[i] = innerp(3, (Mat)[i], Vec);
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
VS2017で実行
コード:
output=Mat[3][3](input)*Vec[3](input)
|4824.00| | 12.00 34.00 82.00| | 2.00|
|3478.00| | 43.00 80.00 -54.00| | 64.00|
|4592.00| | 56.00 31.00 78.00| | 32.00|
続行するには何かキーを押してください . . .
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 00:29
by Math
// これでも良い
コード:
#include<stdio.h>
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
double innerp(int n, double *a, double *b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat[3][3], Vec[3], ReVec[3]);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", ReVec[i], Mat[i][0], Mat[i + 1][1],Mat[i + 2][2], Vec[i];
return 0;
}
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
{
int i;
for (i = 0; i < 3; i++) {
ReVec[i] = innerp(3, Mat[i][0], Vec);// これでも良い
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 00:41
by box
何か計算はしてますけど、それを「どこにも代入していない」ので、そりゃダメですよね。
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 00:45
by Math
「これでも良い」です。(^^;
コード:
#include<stdio.h>
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
double innerp(int n, double *a, double *b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat[3][3], Vec[3], ReVec[3]);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", ReVec[i], Mat[i][0], Mat[i + 1][1],Mat[i + 2][2], Vec[i];
return 0;
}
void MatPro(double Mat[3][3], Vec[3], ReVec[3]);
{
int i;
for (i = 0; i < 3; i++) {
ReVec[i] = innerp(3, Mat[i][0], Vec[3]);// 「これでも良い」
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 01:14
by かずま
Math さんが書きました:シンタックスシュガーを使って書き換えたコードです。
Math さんが書きました:// これでも良い
Math さんが書きました:「これでも良い」です。(^^;
Mathさん、No.5、No.6、No.8のコードはすべてコンパイルエラーになります。
これ以上、質問者を惑わすような投稿はやめてください。
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 01:36
by Math
あ 寝ぼけて 間違えています。後で訂正します・ごめん。！
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 02:49
by Math
＞Mathさん、No.5、No.6、No.8のコードはすべてコンパイルエラーになります。
＞これ以上、質問者を惑わすような投稿はやめてください。
コマンドプロンプトのエディターにVS2017のIDEを使っていてひどい 勘違いしていました。関数に定義からエラーのはずなのにすみませんでした。No.5、No.6、No.8は初歩的におお間違いです。
------------------------------------------------------------------------------
「訂正版」シンタックスシュガーを使って書き換えたコードです。
コード:
#include<stdio.h>
void MatPro(double Mat[3][3], double Vec[3], double ReVec[3]);
double innerp(int n, double* a, double* b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat, Vec, ReVec);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", ReVec[i], Mat[i][0], Mat[i][1], Mat[i][2], Vec[i]);
return 0;
}
void MatPro(double Mat[3][3], double Vec[3], double ReVec[3])
{
int i;
for (i = 0; i < 3; i++) {
ReVec[i] = innerp(3, Mat[i], Vec);
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
コード:
1>------ すべてのリビルド開始: プロジェクト:ConsoleApplication0601, 構成: Debug Win32 ------
1>c1.c
1>ConsoleApplication0601.vcxproj -> D:\z17\ｃ\ConsoleApplication0601\Debug\ConsoleApplication0601.exe
1>ConsoleApplication0601.vcxproj -> D:\z17\ｃ\ConsoleApplication0601\Debug\ConsoleApplication0601.pdb (Partial PDB)
========== すべてリビルド: 1 正常終了、0 失敗、0 スキップ ==========
コード:
output=Mat[3][3](input)*Vec[3](input)
|4824.00| | 12.00 34.00 82.00| | 2.00|
|3478.00| | 43.00 80.00 -54.00| | 64.00|
|4592.00| | 56.00 31.00 78.00| | 32.00|
続行するには何かキーを押してください . . .
Re: 行列とベクトルの積
Posted: 2017年6月03日(土) 09:03
by Math
シンタックスシュガーであるためこれも正しいプログラムです！
コード:
#include<stdio.h>
void MatPro(double Mat[3][3], double Vec[3], double ReVec[3]);
double innerp(int n, double* a, double* b);
int main(void)
{
double Mat[3][3] = { { 12, 34, 82 },{ 43, 80, -54 },{ 56, 31, 78 } };
double Vec[3] = { 2, 64, 32 }, ReVec[3];
int i;
MatPro(Mat, Vec, ReVec);
printf("output=Mat[3][3](input)*Vec[3](input)\n");
for (i = 0; i < 3; i++)
printf("|%7.2f| |%6.2f %6.2f %6.2f| |%6.2f|\n", i[ReVec], i[Mat][0], i[Mat][1], i[Mat][2], i[Vec]);
return 0;
}
void MatPro(double Mat[3][3], double Vec[3], double ReVec[3])
{
int i;
for (i = 0; i < 3; i++) {
i[ReVec] = innerp(3, i[Mat], Vec);
}
}
double innerp(int n, double *a, double *b)//nは次元:*a は行Vector:*b は列Vector
{
double ip = 0;//内積(inner product)
int i;
for (i = 0; i < n; i++)
ip += *(a + i) * *(b + i);
return ip;//内積(inner product)を返す
}
コード:
1>------ すべてのリビルド開始: プロジェクト:ConsoleApplication2, 構成: Debug Win32 ------
1>c1.c
1>ConsoleApplication2.vcxproj -> D:\z17\ｃ\ConsoleApplication2\Debug\ConsoleApplication2.exe
1>ConsoleApplication2.vcxproj -> D:\z17\ｃ\ConsoleApplication2\Debug\ConsoleApplication2.pdb (Partial PDB)
========== すべてリビルド: 1 正常終了、0 失敗、0 スキップ ==========
コード:
output=Mat[3][3](input)*Vec[3](input)
|4824.00| | 12.00 34.00 82.00| | 2.00|
|3478.00| | 43.00 80.00 -54.00| | 64.00|
|4592.00| | 56.00 31.00 78.00| | 32.00|
続行するには何かキーを押してください . . .