Download miễn phí Bà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<
được với mô hình không gian trạng thái của đối tượng điều khiển.
Để tải bản Đầy Đủ của tài liệu, xin Trả lời bài viết này, Mods sẽ gửi Link download cho bạn sớm nhất qua hòm tin nhắn.
Ai cần download tài liệu gì mà không tìm thấy ở đây, thì đăng yêu cầu down tại đây nhé:
Nhận download tài liệu miễn phí
Tóm tắt nội dung tài liệu:
DIEN 9 – K53 – DHBKHNBÀ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<
đượ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
DIEN 9 – K53 – DHBKHN
#include
#include
#include
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= new float [n1+1];
for(int j=1;j
}
}
/**********************************************************************/
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[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.p= new float[x.n+1];
for(int j=1;j
cout << "Phan Tu [" << i << "]" << "[" << j << "]" << " = " ;
DIEN 9 – K53 – DHBKHN
is >> x.p[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
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
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();
DIEN 9 – K53 – DHBKHN
}
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[j] += x1.p[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
if(x.p[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[c];
x.p[c]=temp;
}
break;
}
}
else
break;
}
for(i=1; i<=x.m; i++)
for(j=1; j<=x.n; j++)
{
if(x.p[j]!=0&&i==j)
r=i+1;
while(r<=x.m)
{
DIEN 9 – K53 – DHBKHN
float k=x.p[r][j]/x.p[j];
for(c=1; c<=x.n; c++)
x.p[r][c]=x.p[r][c]-k*x.p[c];
r++;
}
}
for(i=1; i<=x.m; i++)
{
for(j=1; j<=x.n; j++)
if(x.p[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
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
for (int j=2;j
mt c=a^(j-1);
d.p[k][j]=(c*b).p[k][1];
}
return d;
}
DIEN 9 – K53 – DHBKHN
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<
#include “mt.h”
#include
#include
#include
#include
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:
DIEN 9 – K53 – DHBKHN
DIEN 9 – K53 – DHBKHN
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 “mt.h”
#include
#include
#include
#include
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" <
cout << "\nHe Khong Dieu Khien Duoc" << "\n";
else
cout << "\nHe Dieu Khien Duoc" << "\n";
}
Sau khi chạy ta có kết quả sau:
DIEN 9 – K53 – DHBKHN
...