这是个关于如何使用继承、重载运算符、虚函数以实现向量和矩阵运算的例子。是我学习时写的，供大家学习C++基本知识使用。

#include <math.h>
#include <iostream.h>

void error(){
    cout<<"error!"<<endl;
    };
//=============================================================
class vm
{
public:
    vm();
    virtual ~vm();
    void set();
    void print();
    
};


class Vec: public vm                            //继承
{
private:
    double *a;
    
public:
    int sv;
    Vec()
    {   a=new double[1];
        sv=1;
    }
    void set()
    {
        cout<<"N=";
        cin>>sv;
        a=new double[sv];
        for(int i=0;i<sv;i++)
        {
            cout<<"a["<<i+1<<"]=";
            cin>>*(a+i);
            cout<<endl;
        };
    };
    void set(double * v)
    {
        for(int i=0;i<sv;i++)
            *(a+i)=*(v+i);
    }

    ~Vec(){
        delete[]a;
    };

    double &operator[](int i)
    {
        if(i<1||i>sv)error();
        return a[i-1];
    };
    void print(){
        for(int i=0;i<sv;i++)
            cout<<"a["<<i+1<<"]="<<a[i]<<endl;
                };

    void operator -(){
        for(int i=0;i<this->sv;i++)
            *(a+i)=(-(*(a+i)));
    };
    Vec &operator=(Vec &m){
        a=new double[m.sv];
        for(int i=0;i<m.sv;i++)
            a[i]=m.a[i];
    };
};

int operator!=(Vec &v1,Vec &v2){
    if(v1.sv!=v2.sv){return 0;}
    else
        for(int i=1;i<=v1.sv;i++)
            if( v1[i] != v2[i] )return 0;
            else
                return 1;
}

Vec operator+(Vec &m1,Vec &m2){            //加号运算符
    if(m1.sv!=m2.sv)error();
    Vec temp;
    for(int i=1;i<=m1.sv;i++)temp[i]=m1[i]-m2[i];
    return temp;
};
Vec operator-(Vec &m1,Vec &m2){            //减号运算符
    if(m1.sv!=m2.sv)error();
    Vec temp;
    for(int i=1;i<=m1.sv;i++)temp[i]=m1[i]-m2[i];
    return temp;
};
Vec operator*(Vec &m1,int &m2){        //乘号运算符
    Vec temp;
    for(int i=1;i<=m1.sv;i++)
        temp[i]=m1[i]*m2;
    return temp;
};
Vec operator*(int &m2,Vec &m1){
    Vec temp;
    for(int i=1;i<m1.sv;i++)
        temp[i]=m1[i]*m2;
    return temp;
};

class Matrix:public vm                                    //继承
{
private:
    double **m;
    
public:
    int si,sj;
    ~Matrix(){
        delete[]m;
    }
    Matrix(int sii=2,int sjj=2);
    Matrix MatrixI(int tt)
    {
        Matrix temp(tt,tt);
        for(int i=0;i<tt;i++)
            temp[i][i]=1;
        return temp;
    }
    Matrix(Vec &v)
    {
        *m=new double[v.sv];
        for(int i=1;i<=v.sv;i++)
            m[0][i-1]=v[i];
            
    }
    void set(double *mm)
    {
        for(int i=0;i<si;i++)
            for(int j=0;j<sj;j++)
                m[i][j]=mm[i*j+j];
    }

    Vec Matrix::operator[](int i){
        if((i<1)&&(i>si))error();
        Vec temp;
        for(int j=0;j<si;j++)
            temp[j+1]=m[i][j];
        return temp;
    }
    Matrix operator=(const Matrix &m1)            //赋值运算符
    {
        Matrix m2(m1.si,m1.sj);
        for(int i=1;i<=m1.si;i++)
            m2[i]=m1[i];
        return m2;
    };
    Matrix operator-()                            //负号运算符
    {
        for(int i=0;i<si;i++)
            for(int j=0;j<sj;j++)
                m[i][j]=(-m[i][j]);
    };
};

Matrix::Matrix(int si,int sj)
{
    for(int j=0;j<sj;j++)
        for(int i=0;i<si;i++)
            m[i][j]=0;
}


    int operator!=(Matrix &m1,Matrix &m2)
    {
        if ((m1.si!=m2.si)||(m1.sj!=m1.sj))
        {return 0;}
        else
            for(int i=1;i<=m1.si;i++)
                if(m1[i]!=m2[i]){
                return 0;}
                else
                    return 1;
    }

    Matrix operator+(Matrix &m1,Matrix &m2)    //加号运算符
    {
        if((m1.si!=m2.si)||(m1.sj!=m2.sj))error();
        Matrix temp;
        for(int i=1;i<=m1.si;i++)
        temp[i]=m1[i]+m2[i];
        return temp;
    }

    Matrix operator-(Matrix &m1,Matrix &m2)        //减号运算符
    {
        if((m1.si!=m2.si)||(m1.sj!=m2.sj))error();
        Matrix temp;
        for(int i=1;i<=m1.si;i++)
        temp[i]=m1[i]-m2[i];
        return temp;
    }
    Matrix operator *(Matrix& m1,Matrix& m2);    //乘号运算符
    {
        if(m1.sj!=m2.si)error();
        Matrix temp(m1.si,m2.sj);
        for(int j=1;j<=m2.sj;j++)
            for(int i=1;i<=m1.si;i++)
                for(int k=1;k<=m1.j;k++)
                    (temp[i])[j]=(temp[i])[j]+(m1[i])[k]*(m2[k])[j];
        return temp;    
    }
        Matrix operator * (Matrix& m1,int m2);
    {
        Matrix temp(m1.si,m1.sj);
        for(int i=1;i<=m1.si;i++)
            for(int j=1;j<=m1.sj;j++)
                (temp[i])[j]=(m1[i])[j]*m2;
            return temp;
    }

    Matrix operator ^ (Matrix &m1,int m2);//幂运算符号
    Matrix operator / (Matrix &m1,int &n);//除号运算符

void main()
{
   Vec v1,v2;
   Matrix m1,m2;
   v1.set();
   v1.print();
   double t[4]={1.234,2.322,1.7665,5.231};
   v2.set(t);
   v2.print();
   v1=v2*3;
   v1.print();

}