%廣義預測控制_定常GPC算法仿真
%被控模型為y（k）=1.5*y（k-1）-0.7*y（k-2）+u（k-1）+0.5*u（k-2）+e（k），其中e（k）為噪聲
clc;
clear;
%預測模型
A=[0.7-1.51];

B=[0.51];
na=length（A）-1;
nb=length（B）-1;
%預測時域和控制時域
N=5;
Nu=4;
lmd=0.5;
%仿真步數(shù)
tim=100;
%擾動幅度
ampe=1;
%y（k）的最終給定值yr
yr=50;
%柔化因子a
a=0.5;
%一些計算的中間量
Z=[10];%表示z^（-1）
Zj=cell（1N-1）;%Zj中存儲的為z^（-j）次
for?j=1:N-1
????Zj{j}（1）=1;
????Zj{j}（2:j+1）=0;
end
delta=[-11];%△
Adelta=conv（Adelta）;%A△

%迭代求第一個丟番圖方程中的EF
E=cell（1N）;
F=cell（1N）;
E{1}=[1];
F{1}=deconv（polyadd（1-Adelta）Z）;
for?j=1:N-1
????E{j+1}=polyadd（E{j}F{j}（na+1）*Zj{j}）;
???%?F{j+1}=conv（polyadd（F{j}-（F{j}（na+1）*Adelta））Z）;
???F{j+1}=deconv（polyadd（F{j}-（F{j}（na+1）*Adelta））Z）;
end
%求第二個丟番圖方程中的GH
G=cell（1N）;
H=cell（1N）;
for?j=1:N
????EB=conv（E{j}B）;
????H{j}=EB（1:length（EB）-j）;
????G{j}=EB（length（EB）-j+1:length（EB））;
end
%將EFGH左右翻轉(zhuǎn)，因為MATLAB中的多項式系數(shù)是從高次到低次排列，而GPC算法中是相反的，這樣翻轉(zhuǎn)過來使用比較方便
for?j=1:N
????E{j}=fliplr（E{j}）;
????F{j}=fliplr（F{j}）;
????G{j}=f

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????2271??2012-10-28?21:21??GPC\GPC_dingchang.m
?????文件?????????275??2012-09-24?09:21??GPC\polyadd.m
?????目錄???????????0??2012-12-06?18:09??GPC\