两矩阵相加纹理存储器

① C语言中怎么用指针实现两矩阵的相加（矩阵的行数与列数不相等 比如说是两行三列）

最简单的 约定矩阵的行列 然后用二维数组存储矩阵
在加法的时候 转换成一维数组 然后按顺序相加就可以了
比如
int a[2][3], b[2][3];
//input values
add((int *)a, (int *)b, 2*3);

....
void add(int * a, int *b, int size)
{
int i;

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

a[i] += b[i];

}
这个是把a+b的结果存在a中的 当然也可以小小修改一下 加个参数c 然后把结果存在c中

② Matlab中不同维数矩阵相加问题

这样：

>>A=[1,2,3;4,5,6;7,8,9],B=[1,2;3,4]

A=

123

456

789

B=

12

34

>>[rA,cA]=size(A);

>>[rB,cB]=size(B);

>>C=zeros(max(rA,rB),max(cA,cB));

>>C(1:rA,1:cA)=A;

>>C(1:rB,1:cB)=C(1:rB,1:cB)+B;

>>C

C=

243

796

789

>>[2,4,3;8,10,6;7,8,9]

ans=

243

8106

789

(2)两矩阵相加纹理存储器扩展阅读：

注意事项

两个大小不同的矩阵是不能相加的，要实现对应位置相加，需要将纬度较小的扩展为纬度较大的矩阵，没有数据的位置填上零，形成两个大小相同的矩阵，然后再相加。

matlab最大的好处就是不用事先定义变量，按列存储数据，只要每一列的数据个数是不变的，就可以直接使用，直接使用的变量默认为矩阵，只要存储了第一列数据，其行数就确定了，但列数是可变的，不会报错，放心使用，但前提是你的行数一旦确定就不能超过第一次使用的行数，否则会报错，需要另行编制程序。

③ 请C语言高手帮忙编写两个稀疏矩阵相加的程序,急!!!

这个程序能实现矩阵的加减乘。
从中删除你不需要的部分你会吧。
#include<iostream.h>
#include<malloc.h>
#define TRUE 1
#define ERROR 0
#define OK 1
#define MAXISIZE 100
typedef int Elemtype ;
typedef int Status ;
struct Triple
{
int i,j; //行下标，列下标
Elemtype e; //非零元数的值
};
struct TSMatrix
{
Triple data[MAXISIZE+1];
int mu,nu,tu; //矩阵的行数，列数，非零元数
};
Status CreateSMatrix(TSMatrix &M)
{
int i,m,n;
Elemtype e;
Status k;
cout<<"输入矩阵的行数，列数，非零元数:\n";
cin>>M.mu>>M.nu>>M.tu;
M.data[0].i=0;
for(i=1;i<=M.tu;i++)
{
do
{
cout<<"输入第"<<i<<"个非零元素所在的行(1~"<<M.mu<<"),列(1~"<<M.nu<<"),元素值:\n";
cin>>m>>n>>e;
k=0;
if(m<1||m>M.mu||n<1||n>M.nu)
{ k=1;
cout<<"元素的行列数输入错误!\n";}
if(m<M.data[i-1].i||m==M.data[i-1].i&&n<=M.data[i-1].j)
k=1;
}while(k);
M.data[i].i=m;
M.data[i].j=n;
M.data[i].e=e;
}
return OK;
}
void DestroySMatrix(TSMatrix &M)
{
M.mu=0;
M.nu=0;
M.tu=0;
}

void PrintSMatrix(TSMatrix M)
{
int i;
cout<<M.mu<<"行"<<M.nu<<"列"<<M.tu<<"个非零元数\n";
cout<<"行 列 元素值\n";
for(i=1;i<=M.tu;i++)
cout<<M.data[i].i<<" "<<M.data[i].j<<" "<<M.data[i].e<<endl;
}
Status Cmp(int a,int b)
{ if(a<b) return 1;
else if(a==b) return 0;
else return -1;
}
Status AddSMatrix(TSMatrix M,TSMatrix N,TSMatrix &Q)
{
Triple *Mp,*Me,*Np,*Ne,*Qh,*Qe;
if(M.mu!=N.mu||M.nu!=N.nu)
return ERROR;
Q.mu=M.mu;
Q.nu=M.nu;
Mp=&M.data[1];
Np=&N.data[1];
Me=&M.data[M.tu];
Ne=&N.data[N.tu];
Qh=Qe=Q.data;
while(Mp<=Me&&Np<=Ne)
{ Qe++;
switch(Cmp(Mp->i,Np->i))
{
case 1: *Qe=*Mp;
Mp++;
break;
case 0: switch(Cmp(Mp->j,Np->j))
{
case 1: *Qe=*Mp;
Mp++;
break;
case 0: *Qe=*Mp;
Qe->e+=Np->e;
if(!Qe->e)
Qe--;
Mp++;
Np++;
break;
case -1: *Qe=*Np;
Np++;
}
break;
case -1: *Qe=*Np;
Np++;
}
}
if(Mp>Me)
while(Np<=Ne)
{
Qe++;
*Qe=*Np;
Np++;
}
if(Np>Ne)
while(Mp<=Me)
{
Qe++;
*Qe=*Mp;
Mp++;
}
Q.tu=Qe-Qh;
return OK;
}
Status SubSMatrix(TSMatrix M,TSMatrix N,TSMatrix &Q)
{
int i,l;
for(i=1;i<=N.tu;i++)
N.data[i].e*=-1;
l=AddSMatrix(M,N,Q);
return l;
}
Status MultSMatrix(TSMatrix M,TSMatrix N,TSMatrix &Q)
{
int i,j,h=M.mu,l=N.nu,Qn=0;
Elemtype *Qe;
if(M.nu!=N.mu)
return ERROR;
Q.mu=M.mu;
Q.nu=N.nu;
Qe=(Elemtype *)malloc(h*l*sizeof(Elemtype));
for(i=0;i<h*l;i++)
*(Qe+i)=0;
for(i=1;i<=M.tu;i++)
for(j=1;j<=N.tu;j++)
if(M.data[i].j==N.data[j].i)
*(Qe+(M.data[i].i-1)*l+N.data[j].j-1)+=M.data[i].e*N.data[j].e;
for(i=1;i<=M.mu;i++)
for(j=1;j<=N.nu;j++)
if(*(Qe+(i-1)*l+j-1)!=0)
{
Qn++;
Q.data[Qn].e=*(Qe+(i-1)*l+j-1);
Q.data[Qn].i=i;
Q.data[Qn].j=j;
}
free(Qe);
Q.tu=Qn;

return OK;
}
void main()
{char a;
int h;
TSMatrix A,B,C;
cout<<"创建矩阵A:\n";
CreateSMatrix(A);
PrintSMatrix(A);
cout<<"创建矩阵B:\n";
CreateSMatrix(B);
PrintSMatrix(B);
cout<<"选择操作:\n"<<"A.两稀疏矩阵的和\n"<<"B.两稀疏矩阵的差\n"<<"C.两稀疏矩阵的积\n"<<"Q.退出\n";
cin>>a;
while(a!='Q')
{
switch(a)
{
case 'A':
h=AddSMatrix(A,B,C);
if(h==1)
{cout<<"两稀疏矩阵的和为:\n";
PrintSMatrix(C);
}
else cout<<"该两稀疏矩阵不能求和!\n";
break;
case 'B':
h=SubSMatrix(A,B,C);
if(h==1)
{ cout<<"两稀疏矩阵的差为:\n";
PrintSMatrix(C);
}
else cout<<"该两稀疏矩阵不能求差!\n";
break;
case 'C':
h=MultSMatrix(A,B,C);
if(h==1)
{cout<<"两稀疏矩阵的积为:\n";
PrintSMatrix(C);
}
else cout<<"两稀疏矩阵不能求积!\n";
break;
default: cout<<"输入错误!请重新输入\n";
}cin>>a;
}DestroySMatrix(A);
DestroySMatrix(B);
DestroySMatrix(C);
}

④ 用C语言编写程序（两个矩阵相加）

用C语言编写程序（两个矩阵相加）代码如下：

/**

MatrixAddition.c

实现两个矩阵相加

*/

#include<stdio.h>

const int COLS=3;

void printMatirx(int *pArray,int rows,int cols);

void printMatirx2(int (*pArray)[COLS],int rows);

int main()

{

int A[3][4]={{15,10,9,12},

{18,14,8,7},

{16,13,6,11}};

printf("矩阵A= ");

printMatirx(A,3,4);

int B[3][4]={{4,3,5,2},

{0,9,6,1},

{5,7,2,6}};

printf("矩阵B= ");

printMatirx(B,3,4);

int C[3][4];

int i,j;

//矩阵相加：两个矩阵必须行数和列数一样才能相加，

//和矩阵的每个元素分别是两个矩阵对应的元素的和

printf("矩阵A+矩阵B= ");

for(i=0;i<3;i++)//控制行

{

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

{

C[i][j]=A[i][j]+B[i][j];

//printf("%3d",C[i][j]);//输出结果

}

// printf(" ");

}

int (*p)[COLS]=C;

printMatirx2(p,3);

/*结果应为：

19 13 14 14

18 23 14 8

21 20 8 17

*/

printf("矩阵C-矩阵A= ");

for(i=0;i<3;i++)//控制行

{

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

{

printf("%3d",C[i][j]-A[i][j]);//输出结果

}

printf(" ");

}

/*结果应为：

4 3 5 2

0 9 6 1

5 7 2 6

*/

//矩阵数乘：

int D[3][4];

printf("矩阵D: ");

for(i=0;i<3;i++)//控制行

{

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

{

D[i][j]=1;

printf("%d ", D[i][j]);//输出结果

}

printf(" ");

}

int mul;

printf("矩阵D数乘以：");

scanf("%d",&mul);

printf("矩阵D数乘以%d= ",mul);

for(i=0;i<3;i++)//控制行

{

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

{

D[i][j]=D[i][j]*mul;

printf("%d ", D[i][j]);//输出结果

}

printf(" ");

}

return 0;

}

(4)两矩阵相加纹理存储器扩展阅读：

C语言矩阵加减法函数：

void matrix_a(double **a_matrix, const double **b_matrix, const double **c_matrix,int krow, int kline, int ktrl)

////////////////////////////////////////////////////////////////////////////

// a_matrix=b_matrix+c_matrix

// krow :行数

// kline :列数

// ktrl :大于0: 加法 不大于0:减法

////////////////////////////////////////////////////////////////////////////

{

int k, k2;

for (k = 0; k < krow; k++)

{

for(k2 = 0; k2 < kline; k2++)

{

a_matrix[k][k2] = b_matrix[k][k2]

+ ((ktrl > 0) ? c_matrix[k][k2] : -c_matrix[k][k2]);

}

}

}