版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門(mén)版塊排行榜

返回列表

書(shū)尋玉

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 5543.2
散金: 3
紅花: 22
帖子: 478
在線: 241.1小時(shí)
蟲(chóng)號(hào): 2209924
注冊(cè): 2012-12-27
性別: GG
專業(yè): 制造系統(tǒng)與自動(dòng)化

[求助] MIMO系統(tǒng)子空間系統(tǒng)辨識(shí)算法求助已有1人參與

我想問(wèn)下，大家有沒(méi)有做過(guò)MIMO系統(tǒng)的子空間系統(tǒng)辨識(shí)啊，最好是雙輸入雙輸出系統(tǒng)的辨識(shí)，有沒(méi)有相關(guān)的資料共享下，最好是MATLAB的例程。我按書(shū)上的方法編制了MATLAB程序，算出來(lái)不對(duì)，其中B、D怎么用最小二乘法求呢，最好有詳細(xì)的資料提供下啊，不勝感激啊。

回復(fù)此樓

» 猜你喜歡

085600材料與化工調(diào)劑已經(jīng)有10人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有18人回復(fù)
0703化學(xué)求調(diào)劑已經(jīng)有5人回復(fù)
材料292調(diào)劑已經(jīng)有3人回復(fù)
一志愿陜師大生物學(xué)071000，298分，求調(diào)劑已經(jīng)有3人回復(fù)
306求0703調(diào)劑一志愿華中師范已經(jīng)有7人回復(fù)
333求調(diào)劑已經(jīng)有3人回復(fù)
335求調(diào)劑已經(jīng)有4人回復(fù)
336求調(diào)劑已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有4人回復(fù)

1樓 2014-01-03 10:36:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southisland

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 10
在線: 3.7小時(shí)
蟲(chóng)號(hào): 2987725
注冊(cè): 2014-02-23
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

我寫(xiě)的一段程序，可以參考下
/////////////////////step 5 利用最小二乘計(jì)算B、D///////////////////
/*參考[Peter Van Overschee,96]Subspace Identification for Linear Systems:
Theory Implementation Applications,pp.56*/
Matrix<Type> Mi,Li,MM(I*(p*I-n),m),LL(I*(p*I-n),p*I),IGamd(p*I,p+n),BD;
Mi=trT(U2)*R31*pinv(R11);
Li=trT(U2);
for (int i=1;i<=I;i++)
{
for(int j=1;j<=p*I-n;j++)
{
for(int k=1;k<=m;k++)
MM((i-1)*(p*I-n)+j,k)=Mi(j,(i-1)*m+k);
}
}
for(int i=1;i<=I*(p*I-n);i++)
for(int j=1;j<=p*I;j++)
LL(i,j)=double(0);
for (int i=1;i<=I;i++)
{
for(int j=i;j<=I;j++)
{
for(int k=1;k<=p*I-n;k++)
{
for(int t=1;t<=p;t++)
LL((i-1)*(p*I-n)+k,(j-i)*p+t)=Li(k,(j-1)*p+t);
}
}
}
//設(shè)置矩陣IGamd=[Ip,0;0,Gamd]
for (int i=1;i<=p;i++)
{
for(int j=1;j<=p+n;j++)
{
if (j==i)
IGamd(i,j)=double(1);
else
IGamd(i,j)=double(0);
}
}
for (int i=1;i<=p*(I-1);i++)
{
for(int j=1;j<=p;j++)
IGamd(p+i,j)=double(0);
for(int j=1;j<=n;j++)
IGamd(p+i,p+j)=Gamd(i,j);
}

//QRD<Type> qrBD;qrBD.dec(LL*IGamd);BD=qrBD.solve(MM);
BD=pinv(LL*IGamd)*MM; //最小二乘計(jì)算BD=[D;B]

for(int i=0;i<n+p;i++)
{
if(i<p)
D.setRow(BD.getRow(i),i); //提取D矩陣
else
B.setRow(BD.getRow(i),i-p); //提取B矩陣
}
//計(jì)算算法耗時(shí)

贊一下

回復(fù)此樓

2樓2014-02-23 16:50:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southisland

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 10
在線: 3.7小時(shí)
蟲(chóng)號(hào): 2987725
注冊(cè): 2014-02-23
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

找到了當(dāng)時(shí)寫(xiě)的matlab函數(shù)
function [A,B,C,D]=mysubid(Y,U,I,obqr)
%基本的線性開(kāi)環(huán)子空間辨識(shí)算法
%Y：辨識(shí)輸出信號(hào)， U：辨識(shí)輸入信號(hào)， I：Hankel矩陣行數(shù)
%obqr:如果使用QR分解的斜向投影，則置于非0的任意數(shù)
%obqr==0,使用斜向投影
%計(jì)算時(shí)可選直接投影或者基于QR分解的斜向投影
%extract matrices A & C from matric Gam,then Solve least square:
%M=L*X*[B;D]
% 2012-4-11 rewrite 2012-5-11
%2012-7-11 增加可測(cè)噪聲模型的辨識(shí)

%%%START%%%
if (nargin < 4);obqr = 0;end
[m,N]=size(U);p=size(Y,1);
J=N-2*I+1;
% obqr = 1;
%% STEP 1: 構(gòu)造輸入輸出數(shù)據(jù)hankel矩陣塊% HklU=blkhank(U,I,J );
Up=[];Yp=[];Uf=[];Yf=[];
for i=1:I
Up=[Up;U(:,i:J+i-1)];
Uf=[Uf;U(:,I+i:J+I+i-1)];
Yp=[Yp;Y(:,i:J+i-1)];
Yf=[Yf;Y(:,I+i:J+I+i-1)];
end

%% STEP 2：投影并加權(quán)（直接投影和QR分解斜向投影）
%直接正交投影
Wp=[Yp;Up];
if obqr==0
PI_Uf=eye(J)-Uf'*pinv(Uf*Uf')*Uf;
O=Yf*PI_Uf*Wp'*pinv(Wp*PI_Uf*Wp')*Wp;
else
%采用QR分解計(jì)算斜向投影
YuW=[Uf;Wp;Yf];
[yu_q,yu_r]=qr(YuW');
R= yu_r';% R = R(1:2*I*(m+p),1:2*I*(m+p)); % Truncate
Q=yu_q';
O=R(I*p+2*I*m+1:2*I*(m+p),p*I+1:I*(2*p+m))*inv(R(m*I+1

2*m+p)*I,m*I+1

2*m+p)*I))...
*R(m*I+1

2*m+p)*I,1

2*m+p)*I)*Q(1

2*m+p)*I,

;
end

%% STEP 3：SVD分解，選擇系統(tǒng)狀態(tài)階數(shù)
[Uu,S,V] = svd(O);
ss = diag(S)';
lnss=log(ss);
h = bar(1:I*p,lnss,'b');
grid on;title('Singular Value Histogram');xlabel('The Number of Singular Value ');ylabel('Singular Value');
n = input(' System order ? ');%選擇系統(tǒng)階數(shù)

%% STEP 4：由Gam利用其移不變性計(jì)算C和A
U1=Uu(:,1:n);V1=V(1:n,

;
r=sqrt(S(1:n,1:n));
Gam=U1*r;
Xf=r*V1;
C_h=Gam(1:p,1:n);
A_h=pinv(Gam(1:p*(I-1),1:n))*Gam(p+1:p*I,1:n);

%% STEP 5：計(jì)算B和D: 對(duì)于M=L*X*[B;D]采用LS計(jì)算B,D
tao=Gam';%注意：轉(zhuǎn)置矩陣
PI_Gam=eye(size(tao,2))-tao'*inv(tao*tao')*tao;%Gam轉(zhuǎn)置的正交投影算子，PI_Gam*Gam=0
Ms=PI_Gam*Yf*pinv(Uf);
for k=1:I
M(I*p*(k-1)+1:I*p*k,

=Ms(:,m*(k-1)+1:m*k);%構(gòu)造M的列塊矩陣
L(I*p*(k-1)+1:I*p*k,

=[PI_Gam(:,p*(k-1)+1:p*I),zeros(p*I,p*(k-1))];%構(gòu)造L的塊方陣
end
IG = [eye(p),zeros(p,n);zeros(p*(I-1),p),Gam(1

I-1)*p,

];%構(gòu)造【I,0;0,gamma】矩陣
% Solve least squares
sol_bd = (L*IG)\M;
D_h = sol_bd(1:p,

;B_h = sol_bd(p+1:p+n,

;

贊一下

回復(fù)此樓

3樓2014-02-23 16:55:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

書(shū)尋玉

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 5543.2
散金: 3
紅花: 22
帖子: 478
在線: 241.1小時(shí)
蟲(chóng)號(hào): 2209924
注冊(cè): 2012-12-27
性別: GG
專業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

2樓: Originally posted by southisland at 2014-02-23 16:50:39
我寫(xiě)的一段程序，可以參考下
/////////////////////step 5 利用最小二乘計(jì)算B、D///////////////////
/*參考Subspace Identification for Linear Systems:
Theory Implementation Applications,pp.56*/
M ...

非常感謝你，能分享下其他的step嗎？非常感謝啊

贊一下

回復(fù)此樓

4樓2014-02-23 18:48:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southisland

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 10
在線: 3.7小時(shí)
蟲(chóng)號(hào): 2987725
注冊(cè): 2014-02-23
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

引用回帖:

4樓: Originally posted by 書(shū)尋玉 at 2014-02-23 18:48:04
非常感謝你，能分享下其他的step嗎？非常感謝啊...

沒(méi)有其他步驟了，基本的SIM算法就這幾個(gè)步驟的。

[ 發(fā)自小木蟲(chóng)客戶端 ]

贊一下

回復(fù)此樓

5樓2014-02-24 16:45:14

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

書(shū)尋玉

鐵桿木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 2 (幼兒園)
金幣: 5543.2
散金: 3
紅花: 22
帖子: 478
在線: 241.1小時(shí)
蟲(chóng)號(hào): 2209924
注冊(cè): 2012-12-27
性別: GG
專業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

5樓: Originally posted by southisland at 2014-02-24 16:45:14
沒(méi)有其他步驟了，基本的SIM算法就這幾個(gè)步驟的。
...

嗯，謝謝了，是SIMO的嗎？不是MIMO的嗎

贊一下

回復(fù)此樓

6樓2014-02-24 21:33:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southisland

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 10
在線: 3.7小時(shí)
蟲(chóng)號(hào): 2987725
注冊(cè): 2014-02-23
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

引用回帖:

6樓: Originally posted by 書(shū)尋玉 at 2014-02-24 21:33:50
嗯，謝謝了，是SIMO的嗎？不是MIMO的嗎...

子空間辨識(shí)算法最大優(yōu)點(diǎn)在于不用優(yōu)化求解，而且處理多變量系統(tǒng)的辨識(shí)時(shí)算法無(wú)需修改，只是將辨識(shí)數(shù)據(jù)組成向量即可，是一種真正意義上的多變量系統(tǒng)辨識(shí)算法。

[ 發(fā)自小木蟲(chóng)客戶端 ]

贊一下

回復(fù)此樓

7樓2014-02-24 21:50:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

秦時(shí)小沫

新蟲(chóng) (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 11.2
散金: 1
帖子: 53
在線: 4.6小時(shí)
蟲(chóng)號(hào): 1993421
注冊(cè): 2012-09-11
專業(yè): 控制理論與方法

引用回帖:

您好，我想問(wèn)一下，如何將閉環(huán)子空間辨識(shí)算法應(yīng)用到串級(jí)控制系統(tǒng)啊，此時(shí)我們?nèi)绻孀R(shí)內(nèi)環(huán)的模型，如何辨識(shí)呢，望您指導(dǎo)，感激不盡！

贊一下

回復(fù)此樓

8樓2014-03-17 14:54:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

southisland

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 88
帖子: 10
在線: 3.7小時(shí)
蟲(chóng)號(hào): 2987725
注冊(cè): 2014-02-23
專業(yè): 控制理論與方法

【答案】應(yīng)助回帖

引用回帖:

8樓: Originally posted by 秦時(shí)小沫 at 2014-03-17 14:54:55
您好，我想問(wèn)一下，如何將閉環(huán)子空間辨識(shí)算法應(yīng)用到串級(jí)控制系統(tǒng)啊，此時(shí)我們?nèi)绻孀R(shí)內(nèi)環(huán)的模型，如何辨識(shí)呢，望您指導(dǎo)，感激不盡！...

從算法上看，閉環(huán)SIM辨識(shí)需要三組數(shù)據(jù)，即參考輸入信號(hào)r，系統(tǒng)輸出信號(hào)y和控制器輸出信號(hào)u，基于這些測(cè)量得到的值可以使用閉環(huán)算法得到內(nèi)環(huán)模型。具體的實(shí)驗(yàn)設(shè)計(jì)可以參考相關(guān)論文和書(shū)籍。

[ 發(fā)自小木蟲(chóng)客戶端 ]

贊一下

回復(fù)此樓

9樓2014-03-21 09:21:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

lovskyboy416

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 130.8
散金: 20
帖子: 11
在線: 3.2小時(shí)
蟲(chóng)號(hào): 4696812
注冊(cè): 2016-05-16
專業(yè): 電機(jī)與電器

您好，您的mimo子空間辨識(shí)算法能發(fā)給一下做個(gè)參考嗎

發(fā)自小木蟲(chóng)Android客戶端

贊一下

回復(fù)此樓

10樓2017-05-30 20:32:36

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主書(shū)尋玉的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問(wèn)津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 0854電子信息求調(diào)劑 324 +3	Promise-jyl 2026-03-23	3/150	2026-03-23 13:43 by wangkm
[考研] 工科材料085601 279求調(diào)劑 +8	困于星晨 2026-03-17	10/500	2026-03-23 13:05 by 醉在風(fēng)里
[考研] 291求調(diào)劑 +5	孅華 2026-03-22	5/250	2026-03-23 09:20 by haoshis
[考研] 293求調(diào)劑 +3	濤濤Wjt 2026-03-22	5/250	2026-03-22 22:21 by jiangpengfei
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +11	楓橋ZL 2026-03-18	13/650	2026-03-22 20:26 by edmund7
[考研] 307求調(diào)劑 +11	冷笙123 2026-03-17	11/550	2026-03-22 20:16 by edmund7
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) +6	幸運(yùn)哩哩 2026-03-22	10/500	2026-03-22 20:10 by edmund7
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5	困困困困坤坤 2026-03-20	6/300	2026-03-22 17:41 by hxsm
[考研] 306求調(diào)劑 +5	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 298求調(diào)劑一志愿211 +3	上岸6666@ 2026-03-20	3/150	2026-03-22 15:50 by ColorlessPI
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 考研調(diào)劑 +4	來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-21	4/200	2026-03-22 12:15 by 星空星月
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 一志愿華南師大 070300（化學(xué)）304分求調(diào)劑 +3	0703武芊慧雪304 2026-03-18	3/150	2026-03-21 00:48 by JourneyLucky
[考研] 296求調(diào)劑 +6	www_q 2026-03-18	10/500	2026-03-20 23:56 by JourneyLucky
[考研] 求調(diào)劑 +3	@taotao 2026-03-20	3/150	2026-03-20 19:35 by JourneyLucky
[考研] 一志愿南理工085701環(huán)境302求調(diào)劑院校 +3	葵梓衛(wèi)隊(duì) 2026-03-20	3/150	2026-03-20 19:28 by zhukairuo
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu

24小時(shí)熱門(mén)版塊排行榜

[求助] MIMO系統(tǒng)子空間系統(tǒng)辨識(shí)算法求助 已有1人參與

» 猜你喜歡

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

[求助] MIMO系統(tǒng)子空間系統(tǒng)辨識(shí)算法求助已有1人參與