 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2572 | 回復(fù): 9 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
書尋玉鐵桿木蟲 (正式寫手)
|
[求助]
MIMO系統(tǒng)子空間系統(tǒng)辨識(shí)算法求助 已有1人參與
|
|
| 我想問下，大家有沒有做過MIMO系統(tǒng)的子空間系統(tǒng)辨識(shí)啊，最好是雙輸入雙輸出系統(tǒng)的辨識(shí)，有沒有相關(guān)的資料共享下，最好是MATLAB的例程。我按書上的方法編制了MATLAB程序，算出來不對(duì)，其中B、D怎么用最小二乘法求呢，最好有詳細(xì)的資料提供下啊，不勝感激啊。 |
回復(fù)此樓» 猜你喜歡
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
書尋玉
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 5543.2
- 散金: 3
- 紅花: 22
- 帖子: 478
- 在線: 241.1小時(shí)
- 蟲號(hào): 2209924
- 注冊(cè): 2012-12-27
- 性別: GG
- 專業(yè): 制造系統(tǒng)與自動(dòng)化
4樓2014-02-23 18:48:04
southisland
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 88
- 帖子: 10
- 在線: 3.7小時(shí)
- 蟲號(hào): 2987725
- 注冊(cè): 2014-02-23
- 專業(yè): 控制理論與方法
【答案】應(yīng)助回帖
|
我寫的一段程序，可以參考下 /////////////////////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í) |
2樓2014-02-23 16:50:39
southisland
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 88
- 帖子: 10
- 在線: 3.7小時(shí)
- 蟲號(hào): 2987725
- 注冊(cè): 2014-02-23
- 專業(yè): 控制理論與方法
【答案】應(yīng)助回帖
|
找到了當(dāng)時(shí)寫的matlab函數(shù) function [A,B,C,D]=mysubid(Y,U,I,obqr) %基本的線性開環(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 增加可測噪聲模型的辨識(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+12*m+p)*I))...*R(m*I+1 2*m+p)*I,12*m+p)*I)*Q(12*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,; |
3樓2014-02-23 16:55:40
southisland
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 88
- 帖子: 10
- 在線: 3.7小時(shí)
- 蟲號(hào): 2987725
- 注冊(cè): 2014-02-23
- 專業(yè): 控制理論與方法
5樓2014-02-24 16:45:14
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +11 | 手機(jī)用戶 2026-03-17 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 085600材料與化工調(diào)劑 +7 | A-哆啦Z夢(mèng) 2026-03-23 | 12/600 |
|
|
|
[考研] 336求調(diào)劑 +4 | 收到VS 2026-03-20 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +4 | 元夕元 2026-03-20 | 4/200 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運(yùn)來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 266求調(diào)劑 +3 | 哇呼哼呼哼 2026-03-20 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +6 | 想上岸的鯉魚 2026-03-18 | 7/350 |
|
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|