c語言矩陣計算

❶ 用c語言編寫一個矩陣運算的程序，高分！

//矩陣三元組之矩陣相加 相乘
#include <iostream>
using namespace std;
typedef int Elemtype;
#define MAXSIZE 12500 //最大非零元素
typedef struct Triple
{
Elemtype value;
int row,col;
}Triple;

typedef struct TSMatrix
{
Triple data[MAXSIZE+1];
int mu,nu,tu;
}TSMatrix;
TSMatrix T;
void InputMatrix(TSMatrix &T) //輸入t個非零元素
{
cout<<"請輸入稀疏矩陣的信息,(行，列，非零元素個數)"<<endl;
cin>>T.mu>>T.nu>>T.tu;
int i;
cout<<"請輸入非零元素的信息(行，列，值),提醒(下標從1開始)"<<endl;
for(i=1;i<=T.tu;++i)
{
cin>>T.data[i].row>>T.data[i].col>>T.data[i].value;
}
}

void Output(TSMatrix T)
{
cout<<"矩陣的三元組表示(ROW=)"<<T.mu<<" COL="<<T.nu<<"非零個數="<<T.tu<<endl;
int i;
for(i=1;i<=T.tu;++i)
{
cout<<"ROW(行):"<<T.data[i].row<<" COL(列):"<<T.data[i].col<<" Value(值)"<<T.data[i].value<<endl;
}
}

void TransposeSMatrix(TSMatrix M,TSMatrix &T) //矩陣的轉置
{
T.mu=M.nu;T.nu=M.mu;T.tu=M.tu;
int i,j,k=1;
for(i=1;i<=M.nu;++i)
{
for(j=1;j<=M.tu;++j)
if(M.data[j].col==i)
{
T.data[k].row=i;
T.data[k].col=M.data[j].row;
T.data[k].value=M.data[j].value;
++k;
}
}
}

void AddMastrix(TSMatrix M,TSMatrix T,TSMatrix &Q) //矩陣相加
{
int index_a,index_b,i=1,j=1,k=1;
Q.mu=M.mu; Q.nu=M.nu;
while (i<=M.tu&&j<=T.tu)
{
index_a=(M.data[i].row)*(M.data[i].col)+M.data[i].col;
index_b=(T.data[j].row)*(T.data[j].col)+T.data[j].col;
if(index_a<index_b)
{
Q.data[k]=M.data[i];
i++;
k++;
}
else if(index_a>index_b)
{
Q.data[k]=T.data[j];
j++;
k++;
}
else if(index_a==index_b)
{
if((M.data[i].value+T.data[j].value)!=0)
{
Q.data[k]=M.data[i];
Q.data[k].value=M.data[i].value+T.data[j].value;
k++;
}
++i;
++j;
}
}
//復制剩餘元素
for(;i<=M.tu;++i)
{
Q.data[k]=M.data[i];
k++;
}
for(;j<=T.tu;++j)
Q.data[k++]=T.data[j];
Q.tu=k-1;
}

void Multiply(TSMatrix M,TSMatrix T,TSMatrix &Q) //相乘
{
if(M.nu!=T.mu)
{
cerr<<"兩矩陣相乘不合法"<<endl;
return ;
}
int *rowSize=new int[T.mu+1]; //存放每行非零元素的個數
int *rowStart=new int[T.mu+2]; //矩陣每行在三元組開始位置
int *temp=new int[T.nu+1]; //存放結果矩陣中每行的計算結果
int i,Current,k,ROWM,COLM,COLB;
for(i=1;i<=T.mu;i++) rowSize[i]=0;
for(i=1;i<=T.tu;++i) rowSize[T.data[i].row]++;
rowStart[1]=1;
for(i=2;i<=T.mu+1;i++)
rowStart[i]=rowStart[i-1]+rowSize[i-1];
Current=1; k=1;
while (Current<=M.tu)
{
ROWM=M.data[Current].row; //當前三元組數據中元素的行號
for(i=1;i<=T.nu;++i) temp[i]=0;
while (Current<=M.tu&&ROWM==M.data[Current].row)
{
COLM=M.data[Current].col; //當前元素的列號，方便與T矩陣的行號相乘
for(i=rowStart[COLM];i<rowStart[COLM+1];i++) //對應T矩陣中每行的個數
{
COLB=T.data[i].col;
temp[COLB]+=(M.data[Current].value)*(T.data[i].value);
}
Current++;
}
for(i=1;i<=T.nu;i++)
{
if(temp[i]!=0)
{
Q.data[k].row=ROWM;
Q.data[k].col=i;
Q.data[k].value=temp[i];
}
k++;
}
}
Q.mu=M.mu;Q.nu=T.nu;
Q.tu=k-1;
}
int main()
{
TSMatrix T,M,Q,S;
InputMatrix(M);
InputMatrix(T);
cout<<"兩矩陣相乘"<<endl;
Multiply(M,T,Q);
Output(Q);
cout<<"兩矩陣相加"<<endl;
AddMastrix(M,M,S);
Output(S);
system("pause");
return 0;
}

❷ 用c語言實現兩個矩陣相乘怎麼做

1、程序運行輸入數據時，第一行為A矩陣的行列數和B矩陣的行列數，接著分別輸入A、B兩個矩陣的值。

❸ c語言矩陣運算

#include<vector>
#include<iostream>
intmain()
{
usingnamespacestd;
intm,n,p;
vector<vector<int>>A,B,C;
cout<<"pleaseinputm,n,p"<<endl;
cin>>m>>n>>p;
vector<int>tmp;
intttp;
for(inti=0;i<m;i++)
{
for(intj=0;j<n;j++)
{
cin>>ttp;
tmp.push_back(ttp);
}
A.push_back(tmp);
tmp.clear();
}
for(inti=0;i<n;i++)
{
for(intj=0;j<p;j++)
{
cin>>ttp;
tmp.push_back(ttp);
}
B.push_back(tmp);
tmp.clear();
}
for(inti=0;i<m;i++)
{
for(intj=0;j<p;j++)
{
ttp=0;
for(intk=0;k<n;k++)
ttp=ttp+A[i][k]*B[k][j];
tmp.push_back(ttp);
}
C.push_back(tmp);
tmp.clear();
}
for(inti=0;i<m;i++)
{
for(intj=0;j<p;j++)
{
cout<<C[i][j]<<'';
}
cout<<endl;
}
system("pause");
}

現在是2017年10月15日，無意中回答了一道2011年12月09日的題，留念。

❹ c語言矩陣運算問題 緊急

#include<stdio.h>
#include<stdlib.h>

typedef struct Matrix{
int *m;
int col;
int row;
};

Matrix* ConsMatrix(int *a,int cn,int rn){
Matrix *ma=(Matrix*)malloc(sizeof(Matrix));
ma->m=a;
ma->col=cn;
ma->row=rn;
return ma;
}

void PrintMatrix(Matrix *ma){
int i,j;
for(i=0;i<ma->row;i++){
for(j=0;j<ma->col;j++){
printf("%d ",ma->m[i*(ma->col)+j]);
}
printf("\n");
}
}

Matrix* MatrixSub(Matrix*ma,Matrix*mb){
Matrix *mc=(Matrix*)malloc(sizeof(Matrix));

if(ma->col==mb->col&&ma->row==mb->row){

int size=ma->row*ma->col;
int i;
int *d=(int*)malloc(sizeof(int)*size);
mc->m=d;
mc->row=ma->row;
mc->col=ma->col;
for(i=0;i<size;i++){
mc->m[i]=ma->m[i]-mb->m[i];
}
return mc;
}
printf("row or col cant match");
return NULL;
}

int main(){
int a[6]={1,2,3,4,5,6};
int b[6]={2,3,4,5,6,7};
Matrix *mb=ConsMatrix(b,2,3);
Matrix *ma=ConsMatrix(a,2,3);
Matrix *mc=MatrixSub(mb,ma);
PrintMatrix(mc);
system("pause");
return 1;
}

❺ C語言中求兩個矩陣的和

C語言中求兩個矩陣的和具體如下：

#include<stdio.h>

#include<stdlib.h>

#include<time.h>

void main()

{

int a[3][3];

int b[3][3];

int i,j,c,x;

printf(" 求兩個矩陣之和 ");

printf("請選擇輸入方式 ");

printf("0為隨機自主輸入，1為人工輸入 ");

scanf("%d",&x);

if(x==1)

{

printf("請輸入兩個3*3矩陣，以回車結束每個數字的輸入 ");

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

printf("a[%d][%d]=",i,j);

scanf("%d",&a[i][j]);

}

}

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

printf("b[%d][%d]=",i,j);

scanf("%d",&b[i][j]);

}

}

}

if(x==0)

{

srand((unsigned int)time(NULL));

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

a[i][j]=rand()%100+1;

}

}

for(i=0;i<=2;i++)

{
for(j=0;j<=2;j++)

{

printf("a[%d][%d]=",i,j);

printf("%d ",a[i][j]);
}
}

srand((unsigned int)time(NULL));

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

b[i][j]=rand()%100+1;

}

}

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

printf("b[%d][%d]=",i,j);

printf("%d ",b[i][j]);

}

}

}

for(i=0;i<=2;i++)

{

for(j=0;j<=2;j++)

{

c=a[i][j]+b[i][j];

printf(" %d ",c);

}

}

}

❻ 如何用C語言編寫程序求矩陣的值

#include<stdio.h>

void main()
{
//定義矩陣1、2為輸入數據矩陣，3為乘積後的矩陣
int cJuZhen1[4][3],cJuZhen2[3][2],cJuZhen3[4][2];
//定義兩個指針*p1,*p2分別指向矩陣1、2
int *p1,*p2,i,j;
p1=&cJuZhen1[0][0];
p2=&cJuZhen2[0][0];
//輸入矩陣1的數據
printf("請輸入4*3的矩陣A: \n");
for(i=0;i<4;i++)
{
for(j=0;j<3;j++)
{
scanf("%d",cJuZhen1[i][j]);
}
}
printf("\n");
//輸入矩陣2的數據
printf("請輸入3*2的矩陣B: \n");
for(i=0;i<3;i++)
{
for(j=0;j<2;j++)
{
scanf("%d",cJuZhen2[i][j]);
}
}
printf("\n");
//輸出矩陣1
printf("矩陣A:"\n);
for(i=0;i<4;i++)
{
for(j=0;j<3;j++)
{
printf("%d ",cJuZhen1[i][j]);
}
if(j==3)
{
printf("\n");
}
}
printf("\n");
//輸出矩陣2
printf("矩陣B: \n");
for(i=0;i<3;i++)
{
for(j=0;j<2;j++)
{
printf("%d ",cJuZhen2[i][j]);
}
if(j==2)
{
printf("\n");
}
}
printf("\n");
//計算矩陣1、2的乘積
printf("兩矩陣乘積為: \n");
for(i=0;i<4;i++)
{
cJuZhen3[i][0]=(*p1)*(*p2)+(*(p1+1))*(*(p2+2))+(*(p1+2))*(*(p2+4));
cJuZhen3[i][1]=(*p1)*(*(p2+1))+(*(p1+1))*(*(p2+3))+(*(p1+2))*(*(p2+5));
p1=p1+3;
}
//輸出矩陣3（即矩陣乘積）數據
printf("矩陣A與B乘積為: \n");
for(i=0;i<4;i++)
{
for(j=0;j<2;j++)
{
printf("%d ",cJuZhen3[i][j]);
}
if(j==2)
{
printf("\n");
}
}
printf("\n");
}

或者這個#include<iostream>
#include<math.h>
using namespace std;
#define n 3 //三階矩陣
#define N 20
#define err 0.0001

void main()
{
int i,j,k;
double A[n][n],X[n],u,y[n],max;

cout<<"Please input the matrix:\n";
for(i=0;i<n;i++)
for(j=0;j<n;j++)
cin>>A[i][j]; //輸入矩陣
cout<<"Please input the initialized vector:\n";
for(i=0;i<n;i++)
cin>>X[i]; //輸入初始向量

k=1;
u=0;
while(1)
{
max=X[0];
for(i=0;i<n;i++)
{
if(max<X[i]) max=X[i]; //選擇最大值
}

for(i=0;i<n;i++)
y[i]=X[i]/max;

for(i=0;i<n;i++)
{
X[i]=0;
for(j=0;j<n;j++)
X[i]+=A[i][j]*y[j]; //矩陣相乘
}

if(fabs(max-u)<err)
{
cout<<"The eignvalue of A is:"<<max<<"\n";
cout<<"The eignvector of A is:";
for(i=0;i<n;i++)
cout<<X[i]<<" ";
break;
}
else
{
if(k<N) {k=k+1;u=max;}
else {cout<<"error\n";break;}
}

}
}