數(shù)學(xué)建模第四篇論文研讀報(bào)告

第四篇第五題討論根據(jù)前一段時(shí)間的可觀測(cè)量來(lái)設(shè)計(jì)電流值的控制方法。將離散化，得到：(=1,2,……468)進(jìn)一步推算：(=1,2,……468)使每一對(duì)應(yīng)時(shí)刻與，聯(lián)立上式與，得到：(=1,2,……468)，再?gòu)碾x散化的式子推出瞬時(shí)制動(dòng)扭矩與驅(qū)動(dòng)電流產(chǎn)生的扭矩的關(guān)系式：(=1,2,……,468)由于T(n)與時(shí)瞬時(shí)對(duì)應(yīng)關(guān)系，兩者同時(shí)出現(xiàn)，T(n)只由決定，在本時(shí)刻n時(shí)只能得到，不知道，則：(=1,2,……468)，定義該式子為用前一時(shí)間段的制動(dòng)扭矩預(yù)測(cè)本時(shí)間段的電流扭矩的遞推式，但與T(n)有差異，記e（n）為兩者之差。但e（n）無(wú)法得到，，用上一時(shí)刻的誤差e（n-1）代替。(=1,2,……468)根據(jù)驅(qū)動(dòng)電流與其產(chǎn)生的扭矩值之間的正比關(guān)系，得到試驗(yàn)臺(tái)電流控制方法：(=1,2,……468)評(píng)價(jià)的方法：通過(guò)對(duì)路試與模擬的角速度進(jìn)行比較，程序如下：clearn=load('zhuansu.txt');M=load('niuju.txt');w1=2*pi*514/60;M0(1)=M(1);wz(1)=w1;t=4.67;T=0.01;x=0.01:0.01:ty=0:0.01:tfork=1:1:t/Tw(k)=2*pi*n(k)/60;wz(k+1)=wz(k)-M(k)*T/48;endplot(x,w,'-r')holdonplot(y,wz,'-b')第四篇第六題clear,clc,closeallx=load('data.txt');%shuju.txt為那個(gè)excel數(shù)據(jù)；T0=x(:,1);l=length(T0);fori=1:lT2(i)=T0(i);T4(i+1)=13/48*T0(i);endT=zeros(l,1);%改進(jìn)后的方案用T表示T(1)=200;w1=zeros(l,1);%w1：T產(chǎn)生的轉(zhuǎn)速w1(1)=53.83320667;t=x(:,3);fori=1:lwp(i)=x(i,2)*2*pi/60;endfori=2:lw1(i)=w1(i-1)+(T(i-1)-T0(i-1))*0.01/35;T(i)=T4(i)+(T2(i-1)-T4(i-1));ifw1(i)-wp(i)>0.01T(i)=T(i)-300;elseifwp(i)-w1(i)>0.01T(i)=T(i)+300;ende