性无码一区二区三区在线观看,少妇被爽到高潮在线观看,午夜精品一区二区三区,无码中文字幕人妻在线一区二区三区,无码精品国产一区二区三区免费

  • 回復
  • 收藏
  • 點贊
  • 分享
  • 發(fā)新帖

矢量變頻器軟件算法

大家好:
     沒有想到在電源論壇里還有這么多的變頻器愛好者和一些變頻器高手,不過好象本人看到做硬件的多,做軟件的幾乎沒有??
     不知道是沒有人做軟件還是因為沒有開頭打開這個話題,如果是后者的話,那我來開個頭吧,不過也希望更多的從事硬件的朋友多來參與,畢竟沒有硬件軟體也是法實現(xiàn)的.
全部回復(16)
正序查看
倒序查看
2006-05-16 11:21
以下是采用TMS320LF2407A的主體程序算法,希望大家能來支持.
0
回復
2006-05-16 11:26
;--------------------------------------主循環(huán)程序(SET_F由外部輸入)---------------------------------
MAIN_LOOP  
            LDP  #6
LT    SET_F                  ; 將頻率調(diào)節(jié)比(Q15格式)轉換成角頻率
MPY  F_OMEGA               ; Q20
PAC                          ; 積送ACC, Q4格式
SACH  OMEGA,1               ; 保存角頻率,Q5格式
LACC  OMEGA
            LT    SET_F                   ;將頻率調(diào)節(jié)比(Q15格式)轉換成參考電壓
MPY  MAX_V                 ; Q15*Q14=Q29
PAC
SACH  SET_V,1                ; 保存參考電壓幅值,Q14格式
B     MAIN_LOOP              ; 循環(huán)
;--------------------------------------假中斷處理----------------------------------------------------------
PHANTOM  
            CLRC     INTM
            RET
;--------------------------------------T1下溢中斷處理子程序-----------------------------------------
_C_INT2    
            SST   #ST0,ST0_SAVE          ; 保存現(xiàn)場ST0
SST   #ST1,ST1_SAVE          ; 保存ST1
LDP   #5
SACH  ACCH
SACL  ACCL                 ; 保存ACC
SPH    P_HI
SPL   P_LO                    ; 保存 P
MPY   #1                      ; P<=T
SPL   T_SAVE                 ; 保存 T
SAR   AR0,AR0_SAVE          ; 保存AR0
LDP   #224  
LACC  PIVR                  ;讀偏移地址
SUB   #029H                  ; T1 下溢中斷?
CC    T1UF_ISR,EQ           ;是下溢中斷跳T1UF_ISR
REST      
            LDP   #5                     ;否則恢復現(xiàn)場
LAR   AR0, AR0_SAVE          ; 恢復AR0
LT    P_LO                     ;恢復P
MPY   #1  
LPH   P_HI  
LT    T_SAVE                  ;恢復T
LACC  ACCH,16
ADDS  ACCL                  ;恢復ACC
LDP   #0                      ;指向B2
LST   #ST1,ST1_SAVE           ; 恢復ST1
LST   #ST0,ST0_SAVE           ; 恢復ST0
CLRC  INTM                   ; 開中斷
RET                            ; 返回
T1UF_ISR  
            LDP   #232
SPLK   #0200H,EVIFRA          ; 清中斷標志
LDP    #6                      ; 計算轉角增量
LT     OMEGA                  ; Q5
MPY    T_SAMPLE               ;Q5*Q24
PAC                            ;積存ACC, Q13
SFR                            ; 右移一位成Q12格式
ADD     THETA_H,16            ; Q12
ADDS    THETA_L               ;計算絕對位置
SACH    THETA_H               ;保存
SACL    THETA_L
BCND   CHK_UPLIM,GEQ         ;如果大于0檢查上限
ADD     THETA_360,16           ;否則與2π比較,Q12
SACH     THETA_H               ; 保存
B        RND_THETA    
CHK_UPLIM  
            SUB     THETA_360,16           ;與2π比較,Q12
BCND   REST_THETA,LEQ        ; 如果在范圍內(nèi)恢復THETA_H
SACH    THETA_H    
B        RND_THETA
REST_THETA  
            ADD    THETA_360,16         ; 恢復THETA_H
RND_THETA  
            ADD     #1,15                  ; 圓整并保留高字
SACH     THETA_R               ;保存
LACC     #1                      ;假設THETA_H在第一象限
SACL     SS                     ; SIN符號=1
SACL     SC                     ; COS符號=1
LACC     THETA_R
SACL     THETA_M               ; 存入THETA_M
SUB      THETA_90 ;
BCND    E_Q,LEQ              ;在第一象限則跳轉
SPLK     #-1,SC                ; COS符號=-1
LACC     THETA_180
SUB      THETA_R              ; 180-θ
SACL    THETA_M              ;存入THETA_M
BCND    E_Q,GEQ             ;在第二象限則跳轉
SPLK     #-1,SS                ;SIN符號=-1
LACC    THETA_R
SUB      THETA_180           ;THETA-180
SACL    THETA_M             ;存入THETA_M
LACC    THETA_270
SUB      THETA_R
BCND    E_Q,GEQ             ;在第三象限則跳轉
SPLK    #1,SC                 ; COS符號=1
LACC    THETA_360
SUB      THETA_R
SACL    THETA_M             ; 在第四象限
E_Q        
            LT       THETA_M             ; Q12.計算查表索引
MPY     THETA_I              ; Q12*Q9
PAC                           ; Q5
SACH    SIN_INDX            ; Q5
LACC    SIN_INDX,11          ;左移11位
SACH    SIN_INDX            ;相當于右移5位變成Q0格式,即整數(shù)
LACC    SIN_ENTRY           ;查SIN表
ADD     SIN_INDX
TBLR    SIN_THETA
LACC    SIN_END
SUB      SIN_INDX ;
TBLR    COS_THETA ;
LT       SS                    ;查COS表
MPY     SIN_THETA            ;修改符號,Q14
PAC
SACL    SIN_THETA            ; 左移16位保存,Q14
LT       SC
MPY     COS_THETA           ; 修改符號, Q14
PAC
SACL    COS_THETA           ;左移16位保存,Q14
LT       SET_V                ;開始計算UA , UB
MPY     COS_THETA            ; Q14*Q14
PAC                           ;存入ACC, Q12
SACH    UA,1                 ; UA ,Q13格式
MPY     SIN_THETA           ; Q14*Q14
PAC                           ; 存入ACC, Q12
SACH    UB,1                  ; UB, Q13格式
LT      THETA_R               ; Q12格式.確定扇區(qū)
MPY     THETA_S              ; Q12*Q15
PAC
SACH    SECTOR
LACC    SECTOR,5
SACH    SECTOR              ; 相當于右移11位變成Q0格式(整數(shù))
LACC    #DEC_MS            ;逆陣數(shù)據(jù)首地址
ADD     SECTOR,2
SACL    TEMP               ;產(chǎn)生地址指針
LAR     AR0,TEMP           ; 指向逆陣表
LT       UA                  ; Q13格式.計算UA*M(1,1)+UB*M(1,2)
MPY     *+                   ; M(1,1) UA, Q13*Q14
PAC                          ; Q11格式
LT       UB                 ; Q11
MPY      *+                 ; M(1,2) UB, Q13*Q14
APAC                         ; 0.5*C1, Q11
BCND    CMP1BIG0,GEQ      ; 如果大于0繼續(xù)
LACC     #0                  ; 否則0
CMP1BIG0  
            SACH    TEMP                ; 0.5*C1,Q11格式
LT       TEMP                ; Q11格式
MPY     T1_PERIODS          ; Q11*Q5
PAC                          ; Q0格式
ADD      #1,16               ; 防止C1=0
SACH     CMP_1              ; 0.5*C1*TP,Q0格式
LT       UA                  ;計算UA*M(2,1)+UB*M(2,2)
MPY      *+                  ; M(2,1) UA,Q13*Q14
PAC                           ; Q11格式
LT       UB                   ; Q13格式
MPY     *+                   ; M(2,2) UB: Q13*Q14
APAC                          ; 0.5*C2,Q11
BCND   CMP2BIG0,GEQ        ; 如果大于0繼續(xù)
LACC   #0                     ;否則0
CMP2BIG0  
            SACH   TEMP                  ; 0.5*C2,Q11格式
LT       TEMP                 ; Q11格式
MPY     T1_PERIODS          ; Q11*Q5
PAC                           ; Q0格式
ADD     #1,16                 ; 防止C2 = 0
SACH    CMP_2                ; 0.5*C2*TP,Q0格式
LACC    #CCKWISE_          ; ACTRA寄存器值表首地址
ADD      SECTOR            ; 查表地址
TBLR     SVPAT    
LAR      AR0,#ACTRA        ;指向ACTRA
LACC     *                ; 讀ACTRA
AND      #0FFFH          ; 清高四位
OR       SVPAT             ; 高四位賦新值
SACL     *                ; 更新ACTRA
LAR      AR0,#CMPR1     ; 指向CMPR1
LACC     CMP_1
SACL     *+               ; 更新CMPR1, 指向CMPR2
ADD      CMP_2
SACL     *                ; 更新CMPR2
SUB      #500              ;減T1周期值,CMPR2是否超限
BCND    IN_LMT,LEQ       ;沒超跳轉
SPLK     #500,*            ;否則賦T1周期值
IN_LMT    
            B         REST              ; 返回
.DATA                      
;--------------------------------------數(shù)據(jù)段------------------------------------------------------------------
ANGLES_  
            .WORD 01922H                ; π/2, Q12格式
.WORD 03244H                ; π, Q12格式
.WORD 04B66H                ; 3π/2, Q12格式
.WORD 06488H                ; 2π, Q12格式
.WORD 20066           ; 矩陣A的逆陣數(shù)據(jù),每一個逆陣有4個數(shù)據(jù),Q14格式
.WORD –11585               ; 按參考電壓所在的扇區(qū)索引
.WORD 0
.WORD 23170
.WORD 20066
.WORD 11585
.WORD -20066
.WORD 11585
.WORD 0
.WORD 23170
.WORD -20066
.WORD -11585
.WORD -20066
.WORD 11585
.WORD 0
.WORD -23170
.WORD -20066
.WORD -11585
.WORD 20066
.WORD -11585
.WORD 0
.WORD -23170
.WORD 20066
.WORD 11585
CCKWISE_  
             .WORD 0001000000000000B      ;逆時針旋轉時ACTR高4位值,按扇區(qū)索引
.WORD 0011000000000000B
             .WORD 0010000000000000B
             .WORD 0110000000000000B
             .WORD 0100000000000000B
             .WORD 0101000000000000B
SIN_ENTRY_
             .WORD 0                       ; 0-90度SIN值表,Q14格式
             .WORD 286,572,857,1143,1428
             .WORD 1713,1997,2280,2563,2845
           .WORD 3126,3406,3686,3964,4240
           .WORD 4516,4790,5063,5334,5604
           .WORD 5872,6138,6402,6664,6924
           .WORD 7182,7438,7692,7943,8192
           .WORD 8438,8682,8923,9162,9397
           .WORD 9630,9860,10087,10311,10531
           .WORD 10749,10963,11174,11381,11585
           .WORD 11786,11982,12176,12365,12551
           .WORD 12733,12911,13085,13255,13421
           .WORD 13583,13741,13894,14044,14189
           .WORD 14330,14466,14598,14726,14849
           .WORD 14968,15082,15191,15296,15396
           .WORD 15491,15582,15668,15749,15826
           .WORD 15897,15964,16026,16083,16135
           .WORD 16182,16225,16262,16294,16322
           .WORD 16344,16362,16374,16382,16384
.END
0
回復
feifeiyu
LV.3
4
2006-05-22 18:34
樓主,即然是交流,就該整理個flow chart,或者C的也行啊,發(fā)個ASM,大家沒有那么多時間看啊,何況用這種DSP的有多少啊?
0
回復
yrr_2002
LV.1
5
2006-05-23 14:45
@vector_inveter
;--------------------------------------主循環(huán)程序(SET_F由外部輸入)---------------------------------MAIN_LOOP              LDP  #6LT    SET_F                  ;將頻率調(diào)節(jié)比(Q15格式)轉換成角頻率MPY  F_OMEGA              ;Q20PAC                          ;積送ACC,Q4格式SACH  OMEGA,1              ;保存角頻率,Q5格式LACC  OMEGA            LT    SET_F                  ;將頻率調(diào)節(jié)比(Q15格式)轉換成參考電壓MPY  MAX_V                ;Q15*Q14=Q29PACSACH  SET_V,1                ;保存參考電壓幅值,Q14格式B    MAIN_LOOP              ;循環(huán);--------------------------------------假中斷處理----------------------------------------------------------PHANTOM              CLRC    INTM            RET;--------------------------------------T1下溢中斷處理子程序-----------------------------------------_C_INT2                SST  #ST0,ST0_SAVE          ;保存現(xiàn)場ST0SST  #ST1,ST1_SAVE          ;保存ST1LDP  #5SACH  ACCHSACL  ACCL                ;保存ACCSPH    P_HISPL  P_LO                    ;保存PMPY  #1                      ;P
先謝謝了,樓主何不再把算法給大家講一下,比如這個程序是如何工作的,里面的SIN表是如何得來的,基本原理是什么.
0
回復
2006-05-26 16:25
@feifeiyu
樓主,即然是交流,就該整理個flowchart,或者C的也行啊,發(fā)個ASM,大家沒有那么多時間看啊,何況用這種DSP的有多少啊?
人氣太低了,也沒有人搞, 沒意思
0
回復
2006-05-31 23:49
@vector_inveter
人氣太低了,也沒有人搞,沒意思
**此帖已被管理員刪除**
0
回復
topman
LV.2
8
2006-06-22 17:10
@vector_inveter
以下是采用TMS320LF2407A的主體程序算法,希望大家能來支持.
我們用TMS240F2812
0
回復
yaron88
LV.3
9
2006-06-22 22:27
@hejianmonk
**此帖已被管理員刪除**
這個程序不是矢量控制,也不完整.只是個SVPWM實驗程序而已,至多也只能用在變頻空調(diào)上,不過沒有PFC控制恐怕以后空調(diào)上也不能用了
0
回復
jsmmp
LV.3
10
2006-06-28 23:06
@vector_inveter
人氣太低了,也沒有人搞,沒意思
幾年前搞過,現(xiàn)在看這些程序還是一頭霧水.

感覺還是用C語言寫可讀性強.

是否有試過用C寫的矢量控制或直接轉矩控制程序嗎?

不知道是否有支持2407的C語言開發(fā)系統(tǒng)?
0
回復
nihaohong
LV.4
11
2006-06-29 08:13
@jsmmp
幾年前搞過,現(xiàn)在看這些程序還是一頭霧水.感覺還是用C語言寫可讀性強.是否有試過用C寫的矢量控制或直接轉矩控制程序嗎?不知道是否有支持2407的C語言開發(fā)系統(tǒng)?
當然有.TI的CCS支持C編程.
0
回復
yaron88
LV.3
12
2006-06-29 08:35
@nihaohong
當然有.TI的CCS支持C編程.
TI 公司提供的電機控制程序很多都是匯編程序
這個程序是最簡單的了!
0
回復
nihaohong
LV.4
13
2006-06-29 18:23
@yaron88
TI公司提供的電機控制程序很多都是匯編程序這個程序是最簡單的了!
自己用C編程就是了,為什么一定要用TI提供的.
0
回復
yaron88
LV.3
14
2006-06-30 09:33
@nihaohong
自己用C編程就是了,為什么一定要用TI提供的.
老兄,這是2407,不是2812.你用過2407沒有?
0
回復
jsmmp
LV.3
15
2006-07-01 00:17
@yaron88
TI公司提供的電機控制程序很多都是匯編程序這個程序是最簡單的了!
很多人都是靠它起步的啊.

但是時間隔久了,真的容易忘的.

哪里有TI的電機控制C樣本程序下載嗎? 要能讀讀肯定長見識.
0
回復
yaron88
LV.3
16
2006-07-01 11:30
@jsmmp
很多人都是靠它起步的啊.但是時間隔久了,真的容易忘的.哪里有TI的電機控制C樣本程序下載嗎?要能讀讀肯定長見識.
/*regulator reference*/
int w_kp,w_ki;
int i_kp,i_ki;
int i_kp,i_ki;
int if_kp,if_ki;  
             
/*fundmental variables*/
int ua,ub,uc;
int ud,uq;
int ia,ib,ic;

int theta=0;  
int pos=0,pos2;

int i=0,j=0;       //用于設置緩沖時間
int ii=0,jj=0;          //用于速度給定
int iii=0,jjj=0;   //用于求出速度
int iiii=0,jjjj=0;   //用于求出id,iq

int bFirstTime=1;  
int t=5;                    /*標示程序是否在運行*/    
int bStop=0;

long int temp,temp2,temp3;  /*臨時的長整數(shù)*/
int n;                /*臨時的整數(shù)*/      

long int MaxU,MaxI;
int UU;
int error=0;

  
int K[2][5]={{100,200,300,500,1000},  
               {10,15,20,30,50}     /* i_kp*/
              };  
            
int a=0,b=0,c=0,d=0;  

int bAuto=0;
int direct=1;
int kk=8;

void  initial()
{
   asm(" clrc      SXM"); // 符號位擴展
   asm(" clrc     OVM"); // 累加器中結果正常溢出
   asm(" clrc      CNF"); // B0被配置為數(shù)據(jù)存儲空間
   *SCSR1=0x83CE; // CLKIN=6M,CLKOUT=4*CLKIN=24M
   *WDCR=0x0E8;
   *IMR=0x0001; //  允許INT1中斷
   *IFR=0x0FFFF; //  清除全部中斷標志,"寫1清0"
}

void  ADINIT()
{
   *T4CNT=0; // T4計數(shù)器清0
   *T4CON=0X150C; // T4為連續(xù)增計數(shù)模式,32分頻,且選用內(nèi)部時鐘源
   *T4PER=1830; // 設置T4的周期寄存器
   *GPTCONB=0X400; // T4周期中斷標志觸發(fā)AD轉換
   *EVBIMRB=0X01; // 清除EVB中斷標志,寫"1"清0
   *EVBIFRB=0X0FFFF; // 清除EVB中斷標志,寫"1"清0
   *ADCTRL1=0X990; // 采樣時間窗口預定標位ACQ PS3-ACQ PS0為0,
// 轉換時間預定標位CPS為9,AD為啟動/停止模式,排
// 序器為級連工作方式,且禁止特殊的兩種工作模式
   *ADCTRL2=0X8404; // 可以用EVB的一個事件信號觸發(fā)AD轉換,
// 且用中斷模式1
   *MAXCONV=3; // 4通道
   *CHSELSEQ1=0XA720;
}
// 啟動AD轉換子程序(通過啟動定時器4的方式間接啟動)
void  ADSOC()
{
   *T4CON=*T4CON|0X40;// 啟動定時器4
}
// 若是其它中斷則直接返回子程序
void interrupt nothing()
{
   return;
}

void  ValueInit()
{
    
    
    if_ki=0;    /*5411.8->42*//*周期為1/1024秒時,21*/
    w_ki=0;/*2032.9->798*/ /*周期為1/1024秒時,399*/
    i_ki=0;/*654.7*->82*/ /*周期為1/1024秒時,41*/
    
    if_kp=500; /*779.853->3119*/  
    w_kp=200;/*34.56->6949*/
    i_kp=100;/*37.62->2408*/
    
    UU=970;
    w_gd=192;
    id_gd=0;                          
    if_gd=216;  /* 約1.69A */
    MaxU=6925216;
    MaxI=10240; /*20480=10A*/
    
    bAuto=0;
    kk=4;  
}


void  interrupt  adint()
{
          
            int flag;  
           int a_cos,a_sin,b_cos,b_sin,c_cos,c_sin;
           static int id_buf[M],iq_buf[M],if_buf[M],w_buf[16];  
           static long int  sumId=0,sumIq=0,sumIf=0;
            *T4CNT=0X0000;
            *MAINSTART=0;
            *RECSTART=0;
              
           flag=*EVBIFRB&0X01;
           if(flag!=0x01)
           {
             asm(" clrc INTM"); // 允許總中斷
             return; // 如果不是定時器1周期中斷,則直接返回
           }
           
           if(bAuto)
           {
           ii++;
            if(ii==1500)
            {
            a++;
            if(a==5)
            {
              a=0;
              b++;
            }
            if(b==5)
              b=5;
            //w_kp=K[0][a];
            //i_kp=K[1][b];  
              
            ii=0;
            if(jj            {
              w_gd=W[jj]*direct;
              if_gd=IFGD[jj];
              MaxU=MAXU[jj];
              i_kp=IKP[jj];
              kk=KK[jj];
              jj++;
            }
            else
            {
              if(jj==NUM)
              {
                jj=0;
                direct*=-1;
              }      
            }
            }
            }
            *POSITION=0xff00;
           pos=*POSITION;
          pos2=pos;
          
/********************AD采樣值讀取***********************************/          
     
            temp=(*RESULT0)>>6;
            temp=(temp*325)>>10;
            if_fk=temp;
              
            temp=(*RESULT1)>>6;
            temp=((temp-449)*609)>>10;
            ia=temp;  
              
              
            temp=(*RESULT2)>>6;
            temp=((temp-449)*609)>>10;
            ib=temp;
              
            temp=(*RESULT3)>>6;
            temp=((temp-449)*609)>>10;
            ic=temp;
          
/*****************************得到位置*******************************************/
            
            if(pos>N)pos-=N;  //512-1024的數(shù)處理
    theta=pos-197;
    a_sin=theta;
    a_cos=theta+128;
    b_sin=theta-171;
    b_cos=theta-43;
    c_sin=theta+171;
    c_cos=theta-214;
    if(a_sin<0)a_sin+=N;
    if(a_cos<0)a_cos+=N;
    if(b_cos<0)b_cos+=N;
    if(b_sin<0)b_sin+=N;
    if(c_cos<0)c_cos+=N;
    if(c_sin<0)c_sin+=N;
    
/*****************************得到速度*******************************************/
    
    if(iii<16)            
    {
       w_fk=0;
       w_buf[iii]=pos2;
       iii++;
    }
    else
    {
       w_fk=pos2-w_buf[jjj];
       if(w_fk>700)
         w_fk=w_fk-1024;
       if(w_fk<-700)
         w_fk=w_fk+1024;
       w_fk=w_fk<<1;
       w_buf[jjj]=pos2;
       jjj++;
       if(jjj==16) jjj=0;
    }    
    
/********************************ia,ib,ic變換為id,iq****************************/          
    
    temp=ia;
            temp=temp*MySin[a_cos];
            temp2=ib;
            temp2=temp2*MySin[b_cos];
            temp3=ic;
            temp3=temp3*MySin[c_cos];
temp=(temp+temp2+temp3)*2;
id_fk=temp>>8;    /*得到id_fk,本應移14位,定標減掉6位*/  
    
    temp=ia;
            temp=temp*MySin[a_sin];
            temp2=ib;
            temp2=temp2*MySin[b_sin];
            temp3=ic;
            temp3=temp3*MySin[c_sin];
temp=-(temp+temp2+temp3)*2;
iq_fk=temp>>8;     /*得到iq_fk*,本應移14位,定標減掉6位*/

/********************************id,iq,if濾波**************************************/
          
           if(iiii     {
       id_buf[iiii]=id_fk;
       iq_buf[iiii]=iq_fk;
       if_buf[iiii]=if_fk;
       sumId+=id_fk;
       sumIq+=iq_fk;
       sumIf+=if_fk;
       iiii++;
    }
    else
    {
       sumId=sumId+id_fk-id_buf[jjjj];
       sumIq=sumIq+iq_fk-iq_buf[jjjj];
       sumIf=sumIf+if_fk-if_buf[jjjj];
       id_fk=sumId>>MM;
       iq_fk=sumIq>>MM;
       //if_fk=sumIf>>MM;
       id_buf[jjjj]=id_fk;
       iq_buf[jjjj]=iq_fk;  
       if_buf[jjjj]=if_fk;
       jjjj++;
       if(jjjj==M) jjjj=0;
    }
    
    if(bFirstTime==1)
{
   id_fk=0;
               iq_fk=0;
               if_fk=0;
               bFirstTime=0;
            }                    
            
/************************勵磁電流調(diào)節(jié)器********************************************/
        
   /*if按Q.7/10定標*/
       temp=if_gd-if_fk;
    if(temp>16||temp<-16)
    {  
        if_out=temp*if_kp;  
    }
    else
    {
        if_out=temp*if_kp+if_sum;
        if_sum+=temp*if_ki;    
        if(if_sum>20000)if_sum=20000;
        if(if_sum<0)if_sum=0;
    }
    
    if(if_out>20000)if_out=20000;
    if(if_out<1000)if_out=1000;
   /*Ufm=40V,按Q.9/16定標,=20480*/
    temp=(if_out*51)>>12;      /*(uf/Ufm*256+256)*/
n=temp+256;
*FALPHA=MyAcos[n];       /*得到FALPHA*/       
                 
/*********************************速度調(diào)節(jié)器***************************************/
    
    temp=w_gd-w_fk;
    temp2=w_fk-w_pre;  
    if((temp2<16)&&(temp2>-16))  
      if((temp>16||temp<-16))
      {
                 w_out=temp*w_kp;
              }
      else
      {
         w_out=temp*w_kp+w_sum;
                 w_sum+=temp*w_ki;
         if(w_sum>MaxI)w_sum=MaxI;
         if(w_sum<-MaxI)w_sum=-MaxI;
      }
      
    if(w_out>MaxI)w_out=MaxI;
    if(w_out<-MaxI)w_out=-MaxI;  
    w_pre=w_fk;
          
/*********************************電流調(diào)節(jié)器***************************************/     
    /*iq,id按Q.11/16定標*/  
    iq_gd=w_out;
    temp=iq_gd-iq_fk;
            if(temp>256||temp<-256)
            {
      iq_out=temp*i_kp;
    }
    else
    {
      iq_out=temp*i_kp+iq_sum;
      iq_sum+=temp*i_ki;
      if(iq_sum>MaxU)iq_sum=MaxU;
      if(iq_sum<-MaxU)iq_sum=-MaxU;
    }
    
    temp=id_gd-id_fk;
    if(temp>256||temp<-256)
            {
              id_out=temp*i_kp;
            }
    else
    {
      id_out=temp*i_kp+id_sum;
              id_sum+=temp*i_ki;
      if(id_sum>MaxU)id_sum=MaxU;
      if(id_sum<-MaxU)id_sum=-MaxU;
    }
      
/*********************************電壓前饋**************************************/
          
            temp=-823;
            temp=temp*iq_gd;
            temp=temp>>9;
            temp=temp*w_gd;
            ud_qk=temp;
            
            temp=796;
            temp=temp*w_gd;
            temp=temp*if_gd;
            temp2=285;
            temp2=temp2*iq_gd;
            uq_qk=(temp>>3)+(temp2>>1);
            
/************************ud,uq/ua,ub,uc********************/     
/*Um=135 V,相電壓最大值,按Q.8/16定標,=34572*/

            
   ........................
...........................
...........................
...........................

/****************************檢測通訊*********************************************/

            *R_AALPHA=0xff00;
    *R_BALPHA=0xff00;
    *R_CALPHA=0xff00;
    *R_FALPHA=0xff00;
    if((*R_AALPHA!=*AALPHA)||
       (*R_BALPHA!=*BALPHA)||
       (*R_CALPHA!=*CALPHA)||
       (*R_FALPHA!=*FALPHA))
    {
         /*MAINSTOP=0xffff;
         bStop=1;
         DAC=0;*/
         error=error+1;
    }
       
/****************************觀測4個量*********************************************/     
    
            if(j<5000)
            {
              j++;    
            }
            else
            {
              if(i              {    
                buffer1[i]=theta;
                buffer2[i]=w_gd;
                buffer3[i]=w_fk;
                
                buffer4[i]=iq_gd;
                buffer5[i]=iq_fk;
                buffer6[i]=iq_out>>9;
                
                buffer7[i]=id_gd;
                buffer8[i]=id_fk;
                buffer9[i]=id_out>>9;
                
                buffer10[i]=ud_qk>>9;
                buffer11[i]=uq_qk>>9;
                buffer12[i]=-ud;
                buffer13[i]=-uq;
                
                buffer14[i]=ua;
                buffer15[i]=ub;
                buffer16[i]=uc;
                
                buffer17[i]=ia;
                buffer18[i]=ib;
                buffer19[i]=ic;
                
                buffer20[i]=if_gd;
                buffer21[i]=if_fk;
                
                buffer22[i]=if_out;
                
                i++;            
              }
              if(i==SIZE)
              {
                 /*MAINSTOP=0xffff;
         /*bStop=1;*/
                 i=SIZE;
              }
            }
            
/***************************觀測四路AD********************************************/     
        temp=ia<<2;
    temp=temp+2048;
    DA0=temp;
    
            temp=ia<<2;
    temp=temp+2048;
    DA1=temp;
    
    temp=ia<<2;
    temp=temp+2048;
    DA2=temp;
    
            temp=ia<<2;
    temp=temp+2048;
    DA3=temp;

    DA4=0;
/********************************************************************************/
          
            *ADCTRL2=*ADCTRL2|0X4200; // 復位SEQ1,且清除INT FLAG SEQ1標志寫"1"清0
            asm(" NOP");
            asm(" NOP");
            asm(" NOP");
            asm(" NOP");
            t=15-t;
            DAC=t;
            asm(" clrc INTM");
}  

main()
{
   asm(" setc INTM");
   initial( );
   ValueInit( );
   ADINIT( );
   asm(" clrc INTM");
   ADSOC( );  
   while(!(bStop))
   {}
}
這是別人寫的程序,給大家看看(這只是其中一部分)
0
回復
2006-07-08 16:59
@yaron88
/*regulatorreference*/intw_kp,w_ki;inti_kp,i_ki;inti_kp,i_ki;intif_kp,if_ki;              /*fundmentalvariables*/intua,ub,uc;intud,uq;intia,ib,ic;inttheta=0;  intpos=0,pos2;inti=0,j=0;      //用于設置緩沖時間intii=0,jj=0;          //用于速度給定intiii=0,jjj=0;  //用于求出速度intiiii=0,jjjj=0;  //用于求出id,iqintbFirstTime=1;  intt=5;                    /*標示程序是否在運行*/    intbStop=0;longinttemp,temp2,temp3;  /*臨時的長整數(shù)*/intn;                /*臨時的整數(shù)*/      longintMaxU,MaxI;intUU;interror=0;  intK[2][5]={{100,200,300,500,1000},                {10,15,20,30,50}    /*i_kp*/              };              inta=0,b=0,c=0,d=0;  intbAuto=0;intdirect=1;intkk=8;void  initial(){  asm("clrc      SXM");//符號位擴展  asm("clrc    OVM");//累加器中結果正常溢出  asm("clrc      CNF");//B0被配置為數(shù)據(jù)存儲空間  *SCSR1=0x83CE;//CLKIN=6M,CLKOUT=4*CLKIN=24M  *WDCR=0x0E8;  *IMR=0x0001;//  允許INT1中斷  *IFR=0x0FFFF;//  清除全部中斷標志,"寫1清0"}void  ADINIT(){  *T4CNT=0;//T4計數(shù)器清0  *T4CON=0X150C;//T4為連續(xù)增計數(shù)模式,32分頻,且選用內(nèi)部時鐘源  *T4PER=1830;//設置T4的周期寄存器  *GPTCONB=0X400;//T4周期中斷標志觸發(fā)AD轉換  *EVBIMRB=0X01;//清除EVB中斷標志,寫"1"清0  *EVBIFRB=0X0FFFF;//清除EVB中斷標志,寫"1"清0  *ADCTRL1=0X990;//采樣時間窗口預定標位ACQPS3-ACQPS0為0,//轉換時間預定標位CPS為9,AD為啟動/停止模式,排//序器為級連工作方式,且禁止特殊的兩種工作模式  *ADCTRL2=0X8404;//可以用EVB的一個事件信號觸發(fā)AD轉換,//且用中斷模式1  *MAXCONV=3;//4通道  *CHSELSEQ1=0XA720;}//啟動AD轉換子程序(通過啟動定時器4的方式間接啟動)void  ADSOC(){  *T4CON=*T4CON|0X40;//啟動定時器4}//若是其它中斷則直接返回子程序voidinterruptnothing(){  return;}void  ValueInit(){            if_ki=0;    /*5411.8->42*//*周期為1/1024秒時,21*/    w_ki=0;/*2032.9->798*//*周期為1/1024秒時,399*/    i_ki=0;/*654.7*->82*//*周期為1/1024秒時,41*/        if_kp=500;/*779.853->3119*/      w_kp=200;/*34.56->6949*/    i_kp=100;/*37.62->2408*/        UU=970;    w_gd=192;    id_gd=0;                              if_gd=216;  /*約1.69A*/    MaxU=6925216;    MaxI=10240;/*20480=10A*/        bAuto=0;    kk=4;  }void  interrupt  adint(){                      intflag;            inta_cos,a_sin,b_cos,b_sin,c_cos,c_sin;          staticintid_buf[M],iq_buf[M],if_buf[M],w_buf[16];            staticlongint  sumId=0,sumIq=0,sumIf=0;            *T4CNT=0X0000;            *MAINSTART=0;            *RECSTART=0;                        flag=*EVBIFRB&0X01;          if(flag!=0x01)          {            asm("clrcINTM");//允許總中斷            return;//如果不是定時器1周期中斷,則直接返回          }                    if(bAuto)          {          ii++;            if(ii==1500)            {            a++;            if(a==5)            {              a=0;              b++;            }            if(b==5)              b=5;            //w_kp=K[0][a];            //i_kp=K[1][b];                            ii=0;            if(jj>6;            temp=(temp*325)>>10;            if_fk=temp;                          temp=(*RESULT1)>>6;            temp=((temp-449)*609)>>10;            ia=temp;                                          temp=(*RESULT2)>>6;            temp=((temp-449)*609)>>10;            ib=temp;                          temp=(*RESULT3)>>6;            temp=((temp-449)*609)>>10;            ic=temp;          /*****************************得到位置*******************************************/                        if(pos>N)pos-=N;  //512-1024的數(shù)處理    theta=pos-197;    a_sin=theta;    a_cos=theta+128;    b_sin=theta-171;    b_cos=theta-43;    c_sin=theta+171;    c_cos=theta-214;    if(a_sinMM;      id_buf[jjjj]=id_fk;      iq_buf[jjjj]=iq_fk;        if_buf[jjjj]=if_fk;      jjjj++;      if(jjjj==M)jjjj=0;    }        if(bFirstTime==1){  id_fk=0;              iq_fk=0;              if_fk=0;              bFirstTime=0;            }                                /************************勵磁電流調(diào)節(jié)器********************************************/          /*if按Q.7/10定標*/      temp=if_gd-if_fk;    if(temp>16||temp20000)if_sum=20000;        if(if_sum20000)if_out=20000;    if(if_out>12;      /*(uf/Ufm*256+256)*/n=temp+256;*FALPHA=MyAcos[n];      /*得到FALPHA*/                      /*********************************速度調(diào)節(jié)器***************************************/        temp=w_gd-w_fk;    temp2=w_fk-w_pre;      if((temp2-16))        if((temp>16||tempMaxI)w_sum=MaxI;        if(w_sumMaxI)w_out=MaxI;    if(w_out256||tempMaxU)iq_sum=MaxU;      if(iq_sum256||tempMaxU)id_sum=MaxU;      if(id_sum>9;            temp=temp*w_gd;            ud_qk=temp;                        temp=796;            temp=temp*w_gd;            temp=temp*if_gd;            temp2=285;            temp2=temp2*iq_gd;            uq_qk=(temp>>3)+(temp2>>1);            /************************ud,uq/ua,ub,uc********************/    /*Um=135V,相電壓最大值,按Q.8/16定標,=34572*/              ........................................................................................................./****************************檢測通訊*********************************************/            *R_AALPHA=0xff00;    *R_BALPHA=0xff00;    *R_CALPHA=0xff00;    *R_FALPHA=0xff00;    if((*R_AALPHA!=*AALPHA)||      (*R_BALPHA!=*BALPHA)||      (*R_CALPHA!=*CALPHA)||      (*R_FALPHA!=*FALPHA))    {        /*MAINSTOP=0xffff;        bStop=1;        DAC=0;*/        error=error+1;    }      /****************************觀測4個量*********************************************/                    if(j9;                                buffer7[i]=id_gd;                buffer8[i]=id_fk;                buffer9[i]=id_out>>9;                                buffer10[i]=ud_qk>>9;                buffer11[i]=uq_qk>>9;                buffer12[i]=-ud;                buffer13[i]=-uq;                                buffer14[i]=ua;                buffer15[i]=ub;                buffer16[i]=uc;                                buffer17[i]=ia;                buffer18[i]=ib;                buffer19[i]=ic;                                buffer20[i]=if_gd;                buffer21[i]=if_fk;                                buffer22[i]=if_out;                                i++;                          }              if(i==SIZE)              {                /*MAINSTOP=0xffff;        /*bStop=1;*/                i=SIZE;              }            }            /***************************觀測四路AD********************************************/            temp=ia
還是感覺C語言 親切啊.
0
回復
發(fā)