可積系統與數值演算法_可積系統是什麼

A. PID演算法溫控C語言

#include<reg51.h>
#include<intrins.h>
#include<math.h>
#include<string.h>
struct PID {
unsigned int SetPoint; // 設定目標 Desired Value
unsigned int Proportion; // 比例常數 Proportional Const
unsigned int Integral; // 積分常數 Integral Const
unsigned int Derivative; // 微分常數 Derivative Const
unsigned int LastError; // Error[-1]
unsigned int PrevError; // Error[-2]
unsigned int SumError; // Sums of Errors
};
struct PID spid; // PID Control Structure
unsigned int rout; // PID Response (Output)
unsigned int rin; // PID Feedback (Input)
sbit data1=P1^0;
sbit clk=P1^1;
sbit plus=P2^0;
sbit subs=P2^1;
sbit stop=P2^2;
sbit output=P3^4;
sbit DQ=P3^3;
unsigned char flag,flag_1=0;
unsigned char high_time,low_time,count=0;//占空比調節參數
unsigned char set_temper=35;
unsigned char temper;
unsigned char i;
unsigned char j=0;
unsigned int s;
/***********************************************************
延時子程序,延時時間以12M晶振為准,延時時間為30us×time
***********************************************************/
void delay(unsigned char time)
{
unsigned char m,n;
for(n=0;n<time;n++)
for(m=0;m<2;m++){}
}
/***********************************************************
寫一位數據子程序
***********************************************************/
void write_bit(unsigned char bitval)
{
EA=0;
DQ=0; /*拉低DQ以開始一個寫時序*/
if(bitval==1)
{
_nop_();
DQ=1; /*如要寫1，則將匯流排置高*/
}
delay(5); /*延時90us供DA18B20采樣*/
DQ=1; /*釋放DQ匯流排*/
_nop_();
_nop_();
EA=1;
}
/***********************************************************
寫一位元組數據子程序
***********************************************************/
void write_byte(unsigned char val)
{
unsigned char i;
unsigned char temp;
EA=0;
TR0=0;
for(i=0;i<8;i++) /*寫一位元組數據，一次寫一位*/
{
temp=val>>i; /*移位操作，將本次要寫的位移到最低位*/
temp=temp&1;
write_bit(temp); /*向匯流排寫該位*/
}
delay(7); /*延時120us後*/
// TR0=1;
EA=1;
}
/***********************************************************
讀一位數據子程序
***********************************************************/
unsigned char read_bit()
{
unsigned char i,value_bit;
EA=0;
DQ=0; /*拉低DQ，開始讀時序*/
_nop_();
_nop_();
DQ=1; /*釋放匯流排*/
for(i=0;i<2;i++){}
value_bit=DQ;
EA=1;
return(value_bit);
}
/***********************************************************
讀一位元組數據子程序
***********************************************************/
unsigned char read_byte()
{
unsigned char i,value=0;
EA=0;
for(i=0;i<8;i++)
{
if(read_bit()) /*讀一位元組數據，一個時序中讀一次，並作移位處理*/
value|=0x01<<i;
delay(4); /*延時80us以完成此次都時序，之後再讀下一數據*/
}
EA=1;
return(value);
}
/***********************************************************
復位子程序
***********************************************************/
unsigned char reset()
{
unsigned char presence;
EA=0;
DQ=0; /*拉低DQ匯流排開始復位*/
delay(30); /*保持低電平480us*/
DQ=1; /*釋放匯流排*/
delay(3);
presence=DQ; /*獲取應答信號*/
delay(28); /*延時以完成整個時序*/
EA=1;
return(presence); /*返回應答信號，有晶元應答返回0,無晶元則返回1*/
}
/***********************************************************
獲取溫度子程序
***********************************************************/
void get_temper()
{
unsigned char i,j;
do
{
i=reset(); /*復位*/
}while(i!=0); /*1為無反饋信號*/
i=0xcc; /*發送設備定位命令*/
write_byte(i);
i=0x44; /*發送開始轉換命令*/
write_byte(i);
delay(180); /*延時*/
do
{
i=reset(); /*復位*/
}while(i!=0);
i=0xcc; /*設備定位*/
write_byte(i);
i=0xbe; /*讀出緩沖區內容*/
write_byte(i);
j=read_byte();
i=read_byte();
i=(i<<4)&0x7f;
s=(unsigned int)(j&0x0f);
s=(s*100)/16;
j=j>>4;
temper=i|j; /*獲取的溫度放在temper中*/
}
/*====================================================================================================
Initialize PID Structure
=====================================================================================================*/
void PIDInit (struct PID *pp)
{
memset ( pp,0,sizeof(struct PID));
}
/*====================================================================================================
PID計算部分
=====================================================================================================*/
unsigned int PIDCalc( struct PID *pp, unsigned int NextPoint )
{
unsigned int dError,Error;
Error = pp->SetPoint - NextPoint; // 偏差
pp->SumError += Error; // 積分
dError = pp->LastError - pp->PrevError; // 當前微分
pp->PrevError = pp->LastError;
pp->LastError = Error;
return (pp->Proportion * Error // 比例項
+ pp->Integral * pp->SumEror // 積分項
+ pp->Derivative * dError); // 微分項
}
/***********************************************************
溫度比較處理子程序
***********************************************************/
compare_temper()
{
unsigned char i;
if(set_temper>temper)
{
if(set_temper-temper>1)
{
high_time=100;
low_time=0;
}
else
{
for(i=0;i<10;i++)
{ get_temper();
rin = s; // Read Input
rout = PIDCalc ( &spid,rin ); // Perform PID Interation
}
if (high_time<=100)
high_time=(unsigned char)(rout/800);
else
high_time=100;
low_time= (100-high_time);
}
}
else if(set_temper<=temper)
{
if(temper-set_temper>0)
{
high_time=0;
low_time=100;
}
else
{
for(i=0;i<10;i++)
{ get_temper();
rin = s; // Read Input
rout = PIDCalc ( &spid,rin ); // Perform PID Interation
}
if (high_time<100)
high_time=(unsigned char)(rout/10000);
else
high_time=0;
low_time= (100-high_time);
}
}
// else
// {}
}
/*****************************************************
T0中斷服務子程序，用於控制電平的翻轉 ,40us*100=4ms周期
******************************************************/
void serve_T0() interrupt 1 using 1
{
if(++count<=(high_time))
output=1;
else if(count<=100)
{
output=0;
}
else
count=0;
TH0=0x2f;
TL0=0xe0;
}
/*****************************************************
串列口中斷服務程序，用於上位機通訊
******************************************************/
void serve_sio() interrupt 4 using 2
{
/* EA=0;
RI=0;
i=SBUF;
if(i==2)
{
while(RI==0){}
RI=0;
set_temper=SBUF;
SBUF=0x02;
while(TI==0){}
TI=0;
}
else if(i==3)
{
TI=0;
SBUF=temper;
while(TI==0){}
TI=0;
}
EA=1; */
}
void disp_1(unsigned char disp_num1[6])
{
unsigned char n,a,m;
for(n=0;n<6;n++)
{
// k=disp_num1[n];
for(a=0;a<8;a++)
{
clk=0;
m=(disp_num1[n]&1);
disp_num1[n]=disp_num1[n]>>1;
if(m==1)
data1=1;
else
data1=0;
_nop_();
clk=1;
_nop_();
}
}
}
/*****************************************************
顯示子程序
功能：將占空比溫度轉化為單個字元，顯示占空比和測得到的溫度
******************************************************/
void display()
{
unsigned char code number[]={0xfc,0x60,0xda,0xf2,0x66,0xb6,0xbe,0xe0,0xfe,0xf6};
unsigned char disp_num[6];
unsigned int k,k1;
k=high_time;
k=k%1000;
k1=k/100;
if(k1==0)
disp_num[0]=0;
else
disp_num[0]=0x60;
k=k%100;
disp_num[1]=number[k/10];
disp_num[2]=number[k%10];
k=temper;
k=k%100;
disp_num[3]=number[k/10];
disp_num[4]=number[k%10]+1;
disp_num[5]=number[s/10];
disp_1(disp_num);
}
/***********************************************************
主程序
***********************************************************/
main()
{
unsigned char z;
unsigned char a,b,flag_2=1,count1=0;
unsigned char phil[]={2,0xce,0x6e,0x60,0x1c,2};;
TMOD=0x21;
TH0=0x2f;
TL0=0x40;
SCON=0x50;
PCON=0x00;
TH1=0xfd;
TL1=0xfd;
PS=1;
EA=1;
EX1=0;
ET0=1;
ES=1;
TR0=1;
TR1=1;
high_time=50;
low_time=50;
PIDInit ( &spid ); // Initialize Structure
spid.Proportion = 10; // Set PID Coefficients
spid.Integral = 8;
spid.Derivative =6;
spid.SetPoint = 100; // Set PID Setpoint
while(1)
{
if(plus==0)
{
EA=0;
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(plus==0)
{
set_temper++;
flag=0;
}
}
else if(subs==0)
{
for(a=0;a<5;a++)
for(b=0;a<102;b++){}
if(subs==0)
{
set_temper--;
flag=0;
}
}
else if(stop==0)
{
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(stop==0)
{
flag=0;
break;
}
EA=1;
}
get_temper();
b=temper;
if(flag_2==1)
a=b;
if((abs(a-b))>5)
temper=a;
else
temper=b;
a=temper;
flag_2=0;
if(++count1>30)
{
display();
count1=0;
}
compare_temper();
}
TR0=0;
z=1;
while(1)
{
EA=0;
if(stop==0)
{
for(a=0;a<5;a++)
for(b=0;b<102;b++){}
if(stop==0)
disp_1(phil);
// break;
}
EA=1;
}
}

B. 可積系統是什麼

在經典力學有限維Hamiltion系統中，Liouvile定理指出:一個定義在某一區域的N維Hamiltion系統，只要存在N個在該區域上彼此對合的運動積分就是完全可積的，即解可用積分表示。

C. 數字信號處理給定y(n)=x(n)+x(n+1),試判定該系統是否是因果穩定系統,並說明理由!

該系統是否輸入有界啊?你沒說,所以不能判斷是否穩定（穩定系統,需要在域上絕對可和或可積）.但是該系統輸出信號會出現在輸入信號前一時刻（也就是說當還沒有輸入信號激勵時,也會出現輸出響應信號）,所以為非因果系統是肯定.所以可以判定該系統不是因果穩定系統

D. 孤立子與可積系統與數學金融哪個好呢

是的。通過計算機計算和關於淺水波的實驗觀測，表明孤立波碰撞後仍保持各自原來的形狀和速度，猶如粒子，因而稱為孤立子，隨著研究的深入，發現除KdV方程外，還有一系列在應用中十分重要的非線性演化方程，孤立子解反映了自然界的一種相當普遍的非線性現象；並發展了一套求解這類非線性微分方程的強有力的解法，因而受到廣泛的重視。孤立子被應用於粒子物理、固體物理以及各種非線性物理問題中

E. 如何用matlab求積分

一、符號積分
符號積分由函數int來實現。該函數的一般調用格式為：
int(s)：沒有指定積分變數和積分階數時，系統按findsym函數指示的默認變數對被積函數或符號表達式s求不定積分；
int(s,v)：以v為自變數，對被積函數或符號表達式s求不定積分；
int(s,v,a,b)：求定積分運算。a,b分別表示定積分的下限和上限。該函數求被積函數在區間[a,b]上的定積分。a和b可以是兩個具體的數，也可以是一個符號表達式，還可以是無窮(inf)。當函數f關於變數x在閉區間[a,b]上可積時，函數返回一個定積分結果。當a,b中有一個是inf時，函數返回一個廣義積分。當a,b中有一個符號表達式時，函數返回一個符號函數。

2.數值積分的實現方法
基於變步長辛普生法，MATLAB給出了quad函數來求定積分。該函數的調用格式為：
[I,n]=quad('fname',a,b,tol,trace)
基於變步長、牛頓－柯特斯(Newton-Cotes)法，MATLAB給出了quadl函數來求定積分。該函數的調用格式為：
[I,n]=quadl('fname',a,b,tol,trace)
其中fname是被積函數名。a和b分別是定積分的下限和上限。tol用來控制積分精度，預設時取tol=0.001。trace控制是否展現積分過程，若取非0則展現積分過程，取0則不展現，預設時取trace=0。返回參數I即定積分值，n為被積函數的調用次數。

F. 數學與應用數學專業考研想考中科大，請問考哪些科目

620 數學分析
842 線性代數與解析幾何
資料可以到科大科院考研網看一下。

今年的招生簡章已經公布了，目前要考的科目已經確定了，可以放真復習了。。

G. 侏羅紀公園6.2積分系統演算法

如果你想要知道自己的積分密碼，那麼大家可以注意到一下規律1：每9位數就會變更一次頭尾的代碼，但期間頭尾的代碼是不變的如56~64的都是B??U2：中間的字母是按照36進制來遞進的，如37: -load BGBE 38: -load BHCE 從這里可以看出。但是有個別字母由於作者編碼原因刻意漏失，每一層漏失數不一，26個字母加10個數字，合36個，其中有5個是漏失的。3：按照積分獎勵原則，在一定積分內首字母是相同的(數值為31)，如：32~63 第一個字母都是B 64~95 第一個字母都是C20: -load ZVN
30: -load 84D
32: -load BA6D
37: -load BGBE
38: -load BHCE
41: -load BLF3
43: -load BNH3
45: -load BPJ3
55: -load BYST
56: -load B1VU
58: -load B3XU
65: -load CB5K
67: -load CD7K
69: -load CF9K
71: -load CJB9
73: -load CLEA
75: -load CNGA
79: -load CRKA
81: -load CUMZ
82: -load CVN0
84: -load CXP0
86: -load CZR0
88: -load C1T0
90: -load C4WP
92: -load C6YQ
93: -load C7ZQ
95: -load C91Q
97: -load DB3Q
99: -load DD5Q
101: -load DSIV
103: -load DUKV
108: -load DZPV
110: -load D2SL
117: -load D9ZL
119: -load EB1M
121: -load EE4B
126: -load EJ9B
128: -load EMBC
130: -load EPD1
135: -load EUI1
137: -load EWK2
143: -load E3RR
150: -load FBZH

導航:首頁 > 源碼編譯 > 可積系統與數值演算法

可積系統與數值演算法

與可積系統與數值演算法相關的資料