帮忙找个错，关于运算符重载矩阵乘法

Vampierty 2011-12-01 06:48:32

#include "stdafx.h"
#include<iostream.h>
class Matrix
{
public:
Matrix(int mm, int nn)//构造函数初始化data
{
m=mm;n=nn;
int i;
data=new double*[mm];
for(i=0;i<mm;i++)
data[i]=new double[nn];
} //构造M行N列的矩阵
Matrix(const Matrix &src) //重载拷贝构造函数
{
m=src.m;n=src.n;
int i,j;
data=new double*[m];//动态建立二维数组
for(i=0;i<m;i++)
data[i]=new double[n];
for(i=0;i<m;i++)
for(j=0;j<n;j++)
data[i][j]=src.data[i][j];
}
Matrix& operator=(const Matrix &src) //重载"="运算符
{
int i,j;
for(i=0;i<m;i++)
delete []data[i];
delete []data;
m=src.m;n=src.n;
data=new double*[m];//动态建立二维数组
for(i=0;i<m;i++)
data[i]=new double[n];
for(i=0;i<m;i++)
for(j=0;j<n;j++)
data[i][j]=src.data[i][j];
return *this;
}
~Matrix()//
{
for(int i=0;i<m;i++)
delete []data[i];
delete []data;
}
Matrix operator * (Matrix &m2); //矩阵乘法
void display();
void input();
private:
double **data;
int m,n;//矩阵的行数，列数
};//类定义结束
Matrix Matrix::operator *(Matrix &m2)//矩阵乘法的实现
{
Matrix m3(this->m,this->n);
int i,j;
for(i=0;i<m;i++)
for(j=0;j<n;j++)
for(int k=0;k<n;k++)
m3.data[i][j]+=data[i][k]*m2.data[k][j];
return m3;
}
void Matrix::input()
{
for(int i=0;i<m;i++)
for(int j=0;j<n;j++)
cin>>data[i][j];
}
void Matrix::display()
{
int i,j;
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
cout<<data[i][j]<<" ";
}
cout<<endl;
}
}
int main(int argc, char* argv[])
{
int x,y;
cout<<"行数:";
cin>>x;
cout<<"列数:";
cin>>y;
Matrix A(x,y);
cout<<"please input "<<x*y<<"numbers"<<endl;
A.input();
A.display();
cout<<"行数:";
cin>>x;
cout<<"列数:";
cin>>y;
Matrix B(x,y);
cout<<"please input "<<x*y<<"numbers"<<endl;
B.input();
B.display();
Matrix C(y,x);
C=A*B;
C.display();
return 0;
}
运行后得到的是这个样子：
-6.27744e+066 -6.27744e+066 -6.27744e+066
-6.27744e+066 -6.27744e+066 -6.27744e+066
-6.27744e+066 -6.27744e+066 -6.27744e+066
不知道哪里出问题了，这只是我程序的一部分，其他部分如加法、减法、矩阵转置都可以正常运行得到结果的~

...全文

89 4 打赏收藏转发到动态举报

写回复

用AI写文章

4 条回复

切换为时间正序

请发表友善的回复…

发表回复

nuaazdh 2011-12-01

打赏
举报

抱歉，楼主，上面还是有错，改正后：

#include <iostream>



using namespace std;



//Matrix矩阵类

class Matrix

{

public:

    Matrix(int mm, int nn)//构造函数

    {

        m=mm;

        n=nn;

        int i,j;

        data=new double*[mm];

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

            data[i]=new double[nn];

        for(i=0;i<m;i++)//矩阵所有元素清零

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

                data[i][j]=0.0;

    } //构造M行N列的矩阵



    Matrix(const Matrix &src) //拷贝构造函数

    {

        m=src.m;

        n=src.n;

          int i,j;

        data=new double*[m];//动态建立二维数组

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

            data[i]=new double[n];

        for(i=0;i<m;i++)//动态数组赋值

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

                data[i][j]=src.data[i][j];

    }



    ~Matrix()//析构函数

    {

        for(int i=0;i<m;i++)

            delete []data[i];

        delete []data;

    }

    Matrix& operator=(const Matrix &src);//重载"="运算符

    Matrix operator * (const Matrix &m2); //矩阵乘法

    void display();

    void input();

    private:

        double **data;

    int m,n;//矩阵的行数，列数

};//类定义结束



Matrix& Matrix::operator=(const Matrix &src) //重载"="运算符

{

    int i,j;

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

        delete []data[i];

    delete []data;

    m=src.m;n=src.n;

    data=new double*[m];//动态建立二维数组

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

        data[i]=new double[n];

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

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

            data[i][j]=src.data[i][j];

    return *this;

}



//矩阵*运算符重载

Matrix Matrix::operator *(const Matrix &m2)//矩阵乘法的实现

{

    Matrix m3(this->m,m2.n);

    if(this->n!=m2.m)

    {

        cout<<"两矩阵无法进行乘法运算.\n"<<endl;

        exit(0);

    }

    int i,j,k,l;

    for(i=0;i<this->m;i++)

        for(j=0;j<m2.n;j++)

        {

            for(k=0;k<this->n;k++)

             {

                m3.data[i][j]+=this->data[i][k]*m2.data[k][j];

             }



        }

    m3.display();

    return m3;

}



//输入矩阵元素

void Matrix::input()

{

    for(int i=0;i<m;i++)

        for(int j=0;j<n;j++)

            cin>>data[i][j];

}



//显示矩阵元素

void Matrix::display()

{

    int i,j;

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

    {

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

        {

            cout<<data[i][j]<<" ";

        }

        cout<<endl;

    }

}



//主函数

int main(int argc, char* argv[])

{

    int x,y;

    cout<<"矩阵1行数:";

    cin>>x;

    cout<<"矩阵1列数:";

    cin>>y;

    Matrix A(x,y);

    cout<<"请输入矩阵1元素（按行，共"<<x*y<<" 个）"<<endl;

    A.input();

    cout<<"矩阵1："<<endl;

    A.display();

    cout<<"矩阵2行数:";

    cin>>x;

    cout<<"矩阵2列数:";

    cin>>y;

    Matrix B(x,y);

    cout<<"请输入矩阵2元素（按行，共"<<x*y<<" 个）"<<endl;

    B.input();

    cout<<"矩阵2："<<endl;

    B.display();

    Matrix C=A*B;

    cout<<"矩阵1 与 矩阵2 的乘积："<<endl;

    C.display();

    return 0;

}

Vampierty 2011-12-01

打赏
举报

[Quote=引用 2 楼 nuaazdh 的回复:]

C/C++ code
#include <iostream>

using namespace std;

//Matrix矩阵类
class Matrix
{
public:
Matrix(int mm, int nn)//构造函数初始化data
{
m=mm;n=nn;
int i,j;
data=new double……
[/Quote]
正如你所说我确实对矩阵乘法想的太简单了，惭愧啊，以后还得多研究研究。看了你帮忙改的地方对我启发很大，程序也能正常运行得到结果，谢谢啦！

nuaazdh 2011-12-01

打赏
举报

#include <iostream>



using namespace std;



//Matrix矩阵类

class Matrix

{

    public:

    Matrix(int mm, int nn)//构造函数初始化data

    {

        m=mm;n=nn;

        int i,j;

        data=new double*[mm];

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

            data[i]=new double[nn];

        for(i=0;i<m;i++)//矩阵所有元素清零

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

                data[i][j]=0.0;

    } //构造M行N列的矩阵

    Matrix(const Matrix &src) //重载拷贝构造函数

    {

        m=src.m;n=src.n;

          int i,j;

        data=new double*[m];//动态建立二维数组

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

            data[i]=new double[n];

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

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

                data[i][j]=src.data[i][j];

    }

    Matrix& operator=(const Matrix &src) //重载"="运算符

    {

        int i,j;

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

            delete []data[i];

        delete []data;

        m=src.m;n=src.n;

        data=new double*[m];//动态建立二维数组

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

            data[i]=new double[n];

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

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

                data[i][j]=src.data[i][j];

        return *this;

    }

    ~Matrix()//

    {

        for(int i=0;i<m;i++)

            delete []data[i];

        delete []data;

    }

    Matrix operator * (Matrix &m2); //矩阵乘法

    void display();

    void input();

    private:

        double **data;

    int m,n;//矩阵的行数，列数

};//类定义结束



//矩阵*运算符重载

Matrix Matrix::operator *(Matrix &m2)//矩阵乘法的实现

{

    Matrix m3(this->m,m2.n);

    if(this->n!=m2.m)

    {

        cout<<"两矩阵无法进行乘法运算.\n"<<endl;

        exit(0);

    }

    int i,j,k,l;

    for(i=0;i<this->m;i++)

        for(j=0;j<m2.n;j++)

        {

            for(k=0;k<this->n;k++)

                for(l=0;l<m2.m;l++)

                    m3.data[i][j]+=this->data[i][k]*m2.data[l][j];

        }



    return m3;

}



//输入矩阵元素

void Matrix::input()

{

    for(int i=0;i<m;i++)

        for(int j=0;j<n;j++)

            cin>>data[i][j];

}



//显示矩阵元素

void Matrix::display()

{

    int i,j;

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

    {

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

        {

            cout<<data[i][j]<<" ";

        }

        cout<<endl;

    }

}



//主函数

int main(int argc, char* argv[])

{

    int x,y;

    cout<<"矩阵1行数:";

    cin>>x;

    cout<<"矩阵1列数:";

    cin>>y;

    Matrix A(x,y);

    cout<<"please input "<<x*y<<" numbers"<<endl;

    A.input();

    A.display();

    cout<<"矩阵2行数:";

    cin>>x;

    cout<<"矩阵2列数:";

    cin>>y;

    Matrix B(x,y);

    cout<<"please input "<<x*y<<" numbers"<<endl;

    B.input();

    B.display();

      Matrix C(y,x);

    C=A*B;

    C.display();

    return 0;

}

楼主矩阵乘法规则没有弄清除吧

nuaazdh 2011-12-01

打赏
举报

#include <iostream>

using namespace std;

//Matrix矩阵类
class Matrix
{
public:
Matrix(int mm, int nn)//构造函数初始化data
{
m=mm;n=nn;
int i,j;
data=new double*[mm];
for(i=0;i<mm;i++)
data[i]=new double[nn];
for(i=0;i<m;i++)//矩阵所有元素清零
for(j=0;j<n;j++)
data[i][j]=0.0;
} //构造M行N列的矩阵
Matrix(const Matrix &src) //重载拷贝构造函数
{
m=src.m;n=src.n;
int i,j;
data=new double*[m];//动态建立二维数组
for(i=0;i<m;i++)
data[i]=new double[n];
for(i=0;i<m;i++)
for(j=0;j<n;j++)
data[i][j]=src.data[i][j];
}
Matrix& operator=(const Matrix &src) //重载"="运算符
{
int i,j;
for(i=0;i<m;i++)
delete []data[i];
delete []data;
m=src.m;n=src.n;
data=new double*[m];//动态建立二维数组
for(i=0;i<m;i++)
data[i]=new double[n];
for(i=0;i<m;i++)
for(j=0;j<n;j++)
data[i][j]=src.data[i][j];
return *this;
}
~Matrix()//
{
for(int i=0;i<m;i++)
delete []data[i];
delete []data;
}
Matrix operator * (Matrix &m2); //矩阵乘法
void display();
void input();
private:
double **data;
int m,n;//矩阵的行数，列数
};//类定义结束

//矩阵*运算符重载
Matrix Matrix::operator *(Matrix &m2)//矩阵乘法的实现
{
Matrix m3(this->m,m2.n);
if(this->n!=m2.m)
{
cout<<"两矩阵无法进行乘法运算.\n"<<endl;
exit(0);
}
int i,j,k,l;
for(i=0;i<this->m;i++)
for(j=0;j<m2.n;j++)
{
for(k=0;k<this->n;k++)
for(l=0;l<m2.m;l++)
m3.data[i][j]+=this->data[i][k]*m2.data[l][j];
}

return m3;
}

//输入矩阵元素
void Matrix::input()
{
for(int i=0;i<m;i++)
for(int j=0;j<n;j++)
cin>>data[i][j];
}

//显示矩阵元素
void Matrix::display()
{
int i,j;
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
cout<<data[i][j]<<" ";
}
cout<<endl;
}
}

//主函数
int main(int argc, char* argv[])
{
int x,y;
cout<<"矩阵1行数:";
cin>>x;
cout<<"矩阵1列数:";
cin>>y;
Matrix A(x,y);
cout<<"please input "<<x*y<<" numbers"<<endl;
A.input();
A.display();
cout<<"矩阵2行数:";
cin>>x;
cout<<"矩阵2列数:";
cin>>y;
Matrix B(x,y);
cout<<"please input "<<x*y<<" numbers"<<endl;
B.input();
B.display();
Matrix C(y,x);
C=A*B;
C.display();
return 0;
}

《跟朱老师从C高级到C++》属于系列课程《朱有鹏老师热门编程语言系列》的第一部分，主要针对C++语言完全学习。本课程可以看做是嵌入式linux核心课程中《C语言高级专题》课程的延续，在学了C语言前提下继续深度学习...

文章目录类的特殊属性与魔术方法特殊的属性查看属性特殊函数魔术方法实例化可视化可哈希(hash)与等等运算符(==)bool运算符重载1.比较运算符重载2.算术运算符重载3.第二参数运算符重载4.赋值运算符的重载运算符重载的...

方法2.1特殊内建函数与查看属性2.2实例化2.3可视化2.4可哈希(hash)与等等运算符(==)2.5bool3运算符重载3.1比较运算符重载3.2算术运算符重载3.3反向运算符重载3.4赋值运算符的重载4.容器相关方法5.可调用对象6.上下文...

通过重载运算符实现加法运算通过函数实现转置运算 #define _CRT_SECURE_NO_WARNINGS #define SWAP(a,b){a=a+b;b=a-b;a=a-b;} #include<stdio.h> #include<stdlib.h> typedef struct Matrix { int ...

没发生任何改变，原来是使用了Mat4的transformVector 方法，然后又发现Mat4的矩阵有重载 Vec3的 *= 运算符，又用该运算符计算结果还是没有任何改变，有点怀疑人生了。后来又搜索Mat4中所有和Vec3相关的函数，...