Chương trình 2: Viết chương trình minh họa việc sử dụng thư viện ma trận để khảo sát tính điều khiển được với mô hình không gian trạng thái của đối tượng điều khiển Main2.cpp #include “
Trang 1BÀI TẬP LỚN KỸ THUẬT LẬP TRÌNH
Đề bài:
Xây dựng một thư viện về ma trận bằng C++ cho phép cung cấp các chức năng để sử dụng viết chương trình sau:
cin>>m;
mkq=(m1+m2)*(m3-m4); //m1,m2 là các ma trận đối tượng
cout<<mkq;
Viết chương trình minh họa việc sử dụng thư viện ma trận để khảo sát tính điều khiển được với mô hình không gian trạng thái của đối tượng điều khiển
Bài làm:
Xây dựng thư viện các thuật toán ma trận cần thiết:
mt.h
class mt
{
private:
int m,n;
float** p; /* Cap phat bo nho dong cho mang 2 chieu (con tro tro toi con tro)*/ public:
mt();
mt(int m1,int n1);
friend ostream& operator<< (ostream& os,const mt& x); /*Ham in ma tran*/ friend istream& operator>> (istream& is,mt &x); /*Ham nhap ma tan*/
friend mt operator+(const mt& x1,const mt& x2); /*Nap chong toan tu cong
ma tran*/
friend mt operator-(const mt& x1,const mt& x2); /*Nap chong toan tu tru ma tran*/
friend mt operator*(const mt& x1,const mt& x2); /*Nap chong toan tu nhan
ma tran*/
friend mt operator^(const mt& x,const int& n); /*Nap chong toan tu luy thua
ma tran*/
friend int Rank(const mt& x); /*Ham tinh hang ma tran*/
friend mt kalman(const mt& a,const mt& b); /*Ham tao ma tran Kalman de khao sat tinh dieu khien duoc cua he thong*/
};
mt.cpp (Định nghĩa các hàm)
#include “mt.h”
#include<iostream.h>
Trang 2#include<conio.h>
#include<math.h>
#include<iomanip.h>
mt::mt()
{
this->m=0;this->n=0;this->p=NULL;
}
/**********************************************************************/ mt::mt(int m1,int n1)
{
this->m=m1;this->n=n1;
this->p = new float* [m1+1];
for(int i=1;i<=m1+1;i++) {
this->p[i]= new float [n1+1];
for(int j=1;j<n1+1;j++)
this->p[i][j]=0;
}
}
/**********************************************************************/ ostream& operator<< (ostream& os,const mt& x)
{
for (int i=1 ; i<= x.m ; ++i)
{
os << "\n" ;
for (int j=1; j<=x.n; ++j)
os << setw(6) << x.p[i][j] ;
}
os << "\n" ;
return os;
}
/**********************************************************************/ istream& operator>> (istream& is, mt& x)
{
cout << "So Hang: " ;
is >> x.m;
cout << "So Cot: " ;
is >> x.n;
x.p = new float* [x.m+1];
for (int i=1;i<x.m+1;i++)
{
x.p[i]= new float[x.n+1];
for(int j=1;j<x.n+1;j++) {
cout << "Phan Tu [" << i << "]" << "[" << j << "]" << " = " ;
Trang 3is >> x.p[i][j];
} }
return is;
}
/**********************************************************************/
mt operator+(const mt& x1,const mt& x2)
{
if (x1.m != x2.m || x1.n !=x2.n)
{ cout << "\n Khong thuc hien duoc do hai ma tran khong cung cap";
getch();
}
else
{
mt d(x1.m,x1.n);
for(int i=1; i<x1.m+1; i++)
for(int j=1; j<x1.n+1; j++)
d.p[i][j] = x1.p[i][j] + x2.p[i][j];
return d;
}
}
/**********************************************************************/
mt operator-(const mt& x1,const mt& x2)
{
if (x1.m != x2.m || x1.n !=x2.n)
{ cout << "\n Khong thuc hien duoc do hai ma tran khong cung cap";
getch();
}
else
{
mt d(x1.m,x1.n);
for(int i=1; i<x1.m+1; i++)
for(int j=1; j<x1.n+1; j++)
d.p[i][j] = x1.p[i][j] - x2.p[i][j];
return d;
}
}
/**********************************************************************/
mt operator*(const mt& x1,const mt& x2)
{
if (x1.n != x2.m)
{ cout << "\n Khong thuc hien duoc do hai ma tran khong cung cap";
getch();
Trang 4}
else
{
mt d(x1.m,x2.n);
for (int i=1; i < x1.m+1; i++)
for (int j=1; j < x2.n+1; j++)
{
for (int k=1; k < x1.n+1; k++)
d.p[i][j] += x1.p[i][k] * x2.p[k][j];
}
return d;
}
}
/**********************************************************************/ int Rank(const mt& x)
{
int rankmt,i,j,r,c;
rankmt=x.m;
for(i=1; i<=x.m; i++)
for(j=1; j<x.n; j++)
{
if(x.p[i][j]==0&&i==j) {
for(r=i+1; r<=x.m; r++) if(x.p[r][j]!=0)
{
for(c=1; c<=x.m; c++) {
float temp= x.p[r][c];
x.p[r][c]=x.p[i][c];
x.p[i][c]=temp;
} break;
} }
else
break;
}
for(i=1; i<=x.m; i++)
for(j=1; j<=x.n; j++) {
if(x.p[i][j]!=0&&i==j)
r=i+1;
while(r<=x.m) {
Trang 5float k=x.p[r][j]/x.p[i][j];
for(c=1; c<=x.n; c++)
x.p[r][c]=x.p[r][c]-k*x.p[i][c];
r++;
} }
for(i=1; i<=x.m; i++)
{
for(j=1; j<=x.n; j++)
if(x.p[i][j]!=0)
break;
else
if(j==x.n)
rankmt ;
}
return rankmt;
}
/**********************************************************************/
mt operator^(const mt& x,const int& n)
{
mt d(x.m,x.n);
for (int i=1; i <x.m+1; i++)
for (int j=1;j<x.n+1;j++)
d.p[i][j]=x.p[i][j];
for (int k=2; k <=n; k++)
d=d*d;
return d;
}
/**********************************************************************/
mt kalman(const mt& a,const mt& b)
{
mt d(b.m,b.m);
for (int i=1;i<b.m+1;i++)
d.p[i][1]=b.p[i][1];
for (int j=2;j<b.m+1;j++)
for(int k=1;k<b.m+1;k++) {
mt c=a^(j-1);
d.p[k][j]=(c*b).p[k][1];
} return d;
}
Trang 6
Chương trình sử dụng thư viện trên:
Chương trình 1:
cin>>m;
mkq=(m1+m2)*(m3-m4); //m1,m2 là các ma trận đối tượng
cout<<mkq;
Main1.cpp
#include “mt.h”
#include<iostream.h>
#include<conio.h>
#include<math.h>
#include<iomanip.h>
void main()
{
mt m1,m2,m3,m4,mkq;
cout << "\n Nhap ma tran M1 " <<"\n";
cin >> m1;
cout << "\n Nhap ma tran M2 " <<"\n";
cin >> m2;
cout << "\n Nhap ma tran M3 " <<"\n";
cin >> m3;
cout << "\n Nhap ma tran M4 " <<"\n";
cin >> m4;
mkq=(m1+m2)*(m3-m4);
cout << "\n Ma tran ket qua MKQ=(M1+M2)*(M3-M4): " <<"\n"; cout << mkq;
}
Sau khi chạy ta có kết quả sau:
Trang 8Chương trình 2:
Viết chương trình minh họa việc sử dụng thư viện ma trận để khảo sát tính điều khiển được với mô hình không gian trạng thái của đối tượng điều khiển
Main2.cpp
#include “mt.h”
#include<iostream.h>
#include<conio.h>
#include<math.h>
#include<iomanip.h>
void main()
{
mt a,b,p;
int n;
cout << "\n Nhap ma tran A " <<"\n";
cin >> a;
cout << "\n Nhap ma tran B " <<"\n";
cin >> b;
cout << "\n Nhap bac cua he " <<"\n";
cin >> n;
p=kalman(a,b);
cout << "\n Ma tran Kalman: " <<"\n" <<p;
if (Rank(p)!=n)
cout << "\nHe Khong Dieu Khien Duoc" << "\n";
else
cout << "\nHe Dieu Khien Duoc" << "\n";
}
Sau khi chạy ta có kết quả sau: