久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

 找回密碼
 立即注冊(cè)

QQ登錄

只需一步,快速開始

搜索
查看: 24560|回復(fù): 8
打印 上一主題 下一主題
收起左側(cè)

LMI算法設(shè)計(jì)狀態(tài)反饋H∞控制器的一個(gè)例子

[復(fù)制鏈接]
跳轉(zhuǎn)到指定樓層
樓主
ID:77367 發(fā)表于 2015-4-18 20:32 | 只看該作者 回帖獎(jiǎng)勵(lì) |倒序?yàn)g覽 |閱讀模式

這兩天看到很多關(guān)于線性矩陣不等式(Linear MatrixInequality,LMI)算法方面的東西,但是好像沒一個(gè)實(shí)際的例子可以參考,糾結(jié)了兩天做了一個(gè)單級(jí)倒立擺的例子出來,對(duì)有些人說不定還是有用的,PS:這是我畢業(yè)論文的一部分
    一個(gè)倒立擺系統(tǒng)的數(shù)學(xué)模型轉(zhuǎn)化成標(biāo)準(zhǔn)的H∞控制問題模型,即:

給加權(quán)矩陣C1和D12選擇一個(gè)合適的參數(shù)(通常是通過仿真反復(fù)試驗(yàn)得出來),求解下面的一個(gè)LMI,使γ的值取到最小,從而得到一個(gè)最優(yōu)的
態(tài)反饋H∞控制器







    下面是利用MATLAB LMI工具箱解決這個(gè)問題獨(dú)一無二的例子,網(wǎng)上沒有第二份這樣的代碼了。另外仿真還可以使用MATLAB/Simulink。。。。。


%基于LMI算法的倒立擺狀態(tài)反饋H∞控制器設(shè)計(jì)
%State Feedback H∞ controller design based lmi approach
%by 2013/04/03 hemmingway <hemmingway@163.com>

clc
clear all
%----------------------------------------------%
%定義常數(shù)矩陣
A=[0  1       0       0;
   0 -0.0883  0.6293  0;
   0  0       0       1;
   0 -0.2357  27.8285 0] ;
B1=[0 2.3566  0 104.2027]';
B2=[0 0.8832  0 2.3566]';
C1=[0.064 0   0    0;
    0     1e-3  0    0;
    0     0   0.11 0;
    0     0   0    0.01;
    0     0   0    0];
D12=[0 0 0 0 0.01]';
D11=[0 0 0 0 0]';
C2=[1 0 0 0;
    0 0 1 0];
D21=[0 0 0 0]';
D22=[0 0 0 0]';

%
setlmis([]);       %建立一個(gè)LMI
X=lmivar(1,[4,1]); %定義矩陣變量
W=lmivar(2,[1,4]);
r1=lmivar(1,[1,1]);

%%%%%%%
lmiterm([1 1 1 X],A,1,'s');
lmiterm([1 1 1 W],B2,1,'s');
lmiterm([1 2 1 0],B1');
lmiterm([1 2 2 0],-1);
lmiterm([1 3 1 X],C1,1);
lmiterm([1 3 1 W],D12,1);
lmiterm([1 3 2 0],D11);
lmiterm([1 3 3 r1],-1,1);

%
lmiterm([-2 1 1 X],1,1);

lmisys=getlmis;

%%----------------------------solver---------------------------------------
n = decnbr(lmisys);
c = zeros(n,1);
for j=1:n
     [r1j]=defcx(lmisys,j,r1);
      c(j)=trace(r1j);
end
%c=mat2dec(lmisys,zeros(4,4),zeros(1,4),eye(1))
[copt,xopt]=mincx(lmisys,c, [0 0 0 0 0]);

X=dec2mat(lmisys,xopt,X)
W=dec2mat(lmisys,xopt,W)

K=W*X^(-1);
K=K/100    %%控制器,為什么要除以100? 因?yàn)镈12矩陣哪里是0.001,不是傳統(tǒng)的1

r1=dec2mat(lmisys,xopt,r1);  
gammar=r1^(1/2)               %%gammar

%----------------------------------simu------------------------------------
%
w=0.0;
n=1;
Dt=0.01;
t=-0.8;
t0=t;
x=[-0.2 0 0.3 0]';
for i=1:1500
     if t<0
        %t1=4*pi*t;
        t1=4*pi*t;
        x=[1.1*sin(t1); 1.2*cos(t1); 0.5*sin(t1)+1.0*cos(t1); 0];
    else
      u=K*x;                    %%反饋控制
      Dx=A*x+B1*w+B2*u;
      x=x+Dx*Dt;
    end

    Y(:,n)=x;
    t=t+Dt;
    n=n+1;
end

figure(1)
time = (1:n-1)*Dt+t0;
xpos=Y(1,:);
xangle=Y(3,:);

subplot(2,1,1)
plot((1:n-1)*Dt+t0,xpos,'k')
axis([-0.8 10 -1.5 1.5])
grid on
xlabel('time(s)')
ylabel('Cart positon')

subplot(2,1,2)
plot((1:n-1)*Dt+t0,xangle,'k')
axis([-0.8 10 -1.5 1.5])
grid on
xlabel('time(s)')
ylabel('Pendulum')








評(píng)分

參與人數(shù) 1黑幣 +8 收起 理由
tieq1952 + 8 很給力!

查看全部評(píng)分

分享到:  QQ好友和群QQ好友和群 QQ空間QQ空間 騰訊微博騰訊微博 騰訊朋友騰訊朋友
收藏收藏1 分享淘帖 頂 踩

相關(guān)帖子

回復(fù)

使用道具 舉報(bào)

沙發(fā)
ID:103169 發(fā)表于 2016-1-14 11:14 | 只看該作者
K=K/100    %%控制器,為什么要除以100? 因?yàn)镈12矩陣哪里是0.01,不是傳統(tǒng)的1
請(qǐng)問為什么要這樣處理啊
回復(fù)

使用道具 舉報(bào)

板凳
ID:103169 發(fā)表于 2016-1-14 15:10 | 只看該作者
另外,請(qǐng)問你的simulink模塊是怎么設(shè)計(jì)的呢?
回復(fù)

使用道具 舉報(bào)

地板
ID:125815 發(fā)表于 2016-6-8 11:12 | 只看該作者
你好,可以全部看看你的畢業(yè)論文嗎?在知網(wǎng)上可以下載嗎?題目是什么?麻煩了
回復(fù)

使用道具 舉報(bào)

5#
ID:138094 發(fā)表于 2016-8-31 21:36 | 只看該作者
請(qǐng)問,可以看一下你解出的X矩陣么?用這種方法必須要X是對(duì)角化矩陣么?我的仿真不是倒立擺的,是航天器的,但是接出來的X矩陣對(duì)角元素特別大,求解
回復(fù)

使用道具 舉報(bào)

6#
ID:259946 發(fā)表于 2017-12-10 22:41 | 只看該作者
mincx的例子很棒,非常感謝
回復(fù)

使用道具 舉報(bào)

7#
ID:267412 發(fā)表于 2017-12-26 22:08 | 只看該作者
真的很給力!
回復(fù)

使用道具 舉報(bào)

8#
ID:253767 發(fā)表于 2017-12-27 08:00 | 只看該作者
真是好東西,很給力,值得贊佩
回復(fù)

使用道具 舉報(bào)

9#
ID:430188 發(fā)表于 2018-11-20 20:54 | 只看該作者
樓主,你好,我是東北大學(xué)的一名學(xué)生,想看一下您的設(shè)計(jì)文章,不知道題目是什么,或者您的這些參數(shù)是如何設(shè)計(jì)出來的,我現(xiàn)在很需要,謝謝你
回復(fù)

使用道具 舉報(bào)

本版積分規(guī)則

小黑屋|51黑電子論壇 |51黑電子論壇6群 QQ 管理員QQ:125739409;技術(shù)交流QQ群281945664

Powered by 單片機(jī)教程網(wǎng)

快速回復(fù) 返回頂部 返回列表
主站蜘蛛池模板: 谁有毛片网址 | 一极黄色片 | 国产成人精品亚洲男人的天堂 | 一区| 国产aa视频 | 色婷婷香蕉在线一区二区 | 国产一区二区网站 | 伊人久久久 | 精品久久视频 | 欧美视频免费在线观看 | 三级网站免费 | 亚洲久久久 | 日韩欧美黄色片 | 俺去俺来也在线www色官网 | 精品国产一区二区三区久久久蜜月 | 午夜国产福利 | 成人在线视频网站 | 五月天激情国产综合婷婷婷 | 欧美成人极品 | 国产精品96 | 日日夜夜狠狠操 | 91久久 | 福利视频二区 | 日韩精品久久久久久久酒店 | 久久人人爱 | 欧美日韩一区二区三区四区 | 国产在线一区二区 | 色婷婷av一区二区三区之e本道 | 91在线看片 | 四虎久久久 | 日韩av在线免费 | 国产高清网站 | 免费三片在线播放 | 亚洲欧美另类在线 | 91久久久久 | 啪啪毛片 | 久久综合99| 国产特级黄色片 | 性色av一区二区 | 亚洲欧美在线视频 | 99午夜 |