/**************************************************************** 外部晶振8M PA0~3：四位數碼管的位選 PB0~7：數碼管的8位段選 外部中斷0用于計數 定時器0溢出中斷的定時為1ms 說明 ：檢測到水流較小時，繼電器延時1秒關閉 ******************************************************************/ #include<iom16v.h> #include<macros.h> #define uchar unsigned char #define uint unsigned int char led_7[10]={0x3F,0x06,0x5B,0x4F,0x66,0x6D,0x7D,0x07,0x7F,0x6F}; //數碼管段選 char position[4]={0xfe,0xfd,0xfb,0xf7};//數碼管位選 uint sumnum=0; //用于記錄1000ms內進入中斷的次數 uint time=0; //記錄進入比較定時器0的次數 uint num=0; //記錄1ms內進入中斷的次數 uint count=0; //進入外部中斷0的次數 uchar flag; uint sumnum1; //記錄100ms內的數目 /***************************函數聲明***************************/ void delay(); void display(uint m ); void init(); void init_0(); void init_2(); void _delay_us(uint l) { unsigned int i; for(i=0;i<l;i++) { asm("nop"); } } /**************************主函數***********************************/ void main() { init(); init_0(); init_2(); while(sumnum<5) { PORTD=0XBF; segdisplay(sumnum1); } while(1) { segdisplay(sumnum1); } } /*************************掃描數碼管時的延時函數*********************/ void delay() { uchar i,j; for(i=6;i>0;i--) for(j=225;j>0;j--); } /************************數碼管顯示函數*****************************/ void segdisplay( int temp) { int seg[]={0x3f,0x06,0x5b,0x4f,0x66,0x6d,0x7d,0x07,0x7f,0x6f}; int temp1,temp2,temp3,temp4; temp1=temp/1000; temp2=(temp/100)%10; temp3=(temp/10)%10; temp4=temp%10; DDRB=0xff; DDRA|=0x0f; PORTA=~BIT(3); PORTB=seg[temp1]; _delay_us(100); PORTA=~BIT(2); PORTB=seg[temp2]; _delay_us(100); PORTA=~BIT(1); PORTB=seg[temp3]; _delay_us(100); PORTA=~BIT(0); PORTB=seg[temp4]; _delay_us(100); } /***********************管腳初始化函數*********************/ void init() { DDRD|=0X40; //PD4 設置為輸出 PORTD=0XBF; DDRA=0XFF; DDRB=0XFF; PORTA=0XFF; PORTB=0XFF; } /***********************外部中斷0初始化*********************/ void init_0() { MCUCR=0X02; //INT0為下降沿觸發 GICR=0X40; //使能INT0中斷 SREG=0X80; //使能總中斷 } /**********************定時器2初始化***********************/ void init_2() { TCCR0=0x03; // 內部時鐘，64 分頻（8M/64=125KHz） TCNT0=0x83; //裝初值 TIMSK=0x01; // 允許 T/C0溢出中斷中斷 } /***********************外部中斷0子函數********************/ #pragma interrupt_handler int0_isr:2 void int0_isr(void) { count++; } /*********************定時計數器0溢出中斷子函數*****************/ #pragma interrupt_handler int0_over:10 void int0_over(void) { TCNT0=0x83; //重裝初值 if((time%100) == 0) sumnum1 = num; if(time == 1000) { sumnum=num; if(sumnum<10) { if((flag==1)&&(sumnum<10)) { PORTD=0XFF; flag=0; } flag++; } else PORTD=0XBF; num=0; time=0; } num+=count; count=0; ++time; }
標簽： C語言
上傳時間： 2016-03-09
上傳用戶：彥 yan
最簡單的阿爾法 貝塔 伽馬 濾波程序 設置了初值 和真實值 并對真實值和估計值進行了比較
標簽： 卡爾曼濾波
上傳時間： 2016-10-21
上傳用戶：260970449
針對傳統的系統級標定方法狀態變量維數高、標定參數可觀測性差等特點，提出一種十位置系統級標 定方法。該方法以分立式標定結果為初值，以速度誤差和姿態誤差作為觀測量，合并加速度計標定誤差和光 纖陀螺標定誤差，降低Kalman濾波器維數。設計10個位置對SINS標定誤差進行估計，然后將估計值進行解 耦，計算SINS標定參數。仿真和轉臺實驗結果表明：十位置系統級標定方法可以一次性標定出標度因數、安 裝誤差和零位等24個標定參數。
上傳時間： 2016-10-23
上傳用戶：260970449
/*import java.util.Scanner; //主類 public class student122 { //主方法 public static void main(String[] args){ //定義7個元素的字符數組 String[] st = new String[7]; inputSt(st); //調用輸入方法 calculateSt(st); //調用計算方法 outputSt(st); //調用輸出方法 } //其他方法 //輸入方法 private static void inputSt(String st[]){ System.out.println("輸入學生的信息："); System.out.println("學號 姓名 成績1，2，3"); //創建鍵盤輸入類 Scanner ss = new Scanner(System.in); for(int i=0; i<5; i++){ st[i] = ss.next(); //鍵盤輸入1個字符串 } } //計算方法 private static void calculateSt(String[] st){ int sum = 0; //總分賦初值 int ave = 0; //平均分賦初值 for(int i=2;i<5;i++) { /計總分，字符變換成整數后進行計算 sum += Integer.parseInt(st[i]); } ave = sum/3; //計算平均分 //整數變換成字符后保存到數組里 st[5] = String.valueOf(sum); st[6] = String.valueOf(ave); } //輸出方法 private static void outputSt(String[] st){ System.out.print("學號 姓名 "); //不換行 System.out.print("成績1 成績2 成績3 "); System.out.println("總分 平均分");//換行 //輸出學生信息 for(int i=0; i<7; i++){ //按格式輸出，小于6個字符，補充空格 System.out.printf("%6s", st[i]); } System.out.println(); //輸出換行 } }*/ import java.util.Scanner; public class student122 { public static void main(String[] args) { // TODO 自動生成的方法存根 String[][] st = new String[3][8]; inputSt(st); calculateSt(st); outputSt(st); } //輸入方法 private static void inputSt(String st[][]) { System.out.println("輸入學生信息："); System.out.println("班級 學號 姓名 成績：數學 物理 化學"); //創建鍵盤輸入類 Scanner ss = new Scanner(System.in); for(int j = 0; j < 3; j++) { for(int i = 0; i < 6; i++) { st[j][i] = ss.next(); } } } //輸出方法 private static void outputSt(String st[][]) { System.out.println("序號 班級 學號 姓名 成績：數學 物理 化學 總分 平均分"); //輸出學生信息 for(int j = 0; j < 3; j++) { System.out.print(j+1 + ":"); for(int i = 0; i < 8; i++) { System.out.printf("%6s", st[j][i]); } System.out.println(); } } //計算方法 private static void calculateSt(String[][] st) { int sum1 = 0; int sum2 = 0; int sum3 = 0; int ave1 = 0; int ave2 = 0; int ave3 = 0; for(int i = 3; i < 6; i++) { sum1 += Integer.parseInt(st[0][i]); } ave1 = sum1/3; for(int i = 3; i < 6; i++) { sum2 += Integer.parseInt(st[1][i]); } ave2 = sum2/3; for(int i = 3; i < 6; i++) { sum3 += Integer.parseInt(st[2][i]); } ave3 = sum3/3; st[0][6] = String.valueOf(sum1); st[1][6] = String.valueOf(sum2); st[2][6] = String.valueOf(sum3); st[0][7] = String.valueOf(ave1); st[1][7] = String.valueOf(ave2); st[2][7] = String.valueOf(ave3); } }
上傳時間： 2017-03-17
上傳用戶：simple
用迭代法求方程 x3-x2-1=0 在[1.3, 1.6]內的一個實根,選初值x0 =1.3,迭代一步。
標簽： gggggg
上傳時間： 2017-05-09
上傳用戶：mic0000
對于非線性方程組F(X)=0，用fsolve函數求其數值解。fsolve函數的調用格式為： X=fsolve('fun',X0,option) 其中X為返回的解，fun是用于定義需求解的非線性方程組的函數文件名，X0是求根過程的初值，option為最優化工具箱的選項設定。最優化工具箱提供了20多個選項，用戶可以使用optimset命令將它們顯示出來。如果想改變其中某個選項，則可以調用optimset()函數來完成。例如，Display選項決定函數調用時中間結果的顯示方式，其中‘off’為不顯示，‘iter’表示每步都顯示，‘final’只顯示最終結果。optimset(‘Display’,‘off’)將設定Display選項為‘off’。
標簽： MATLAB fsolve 非線性 方程 代碼 源程序
上傳時間： 2017-07-31
上傳用戶：ChrisQQQ
四階龍格-庫塔法源程序代碼，用于求解一階或高階常微分方程(組)的初值問題。
標簽： 源程序
上傳時間： 2017-12-29
上傳用戶：liudansn
采用四階龍格庫塔法求解初值微分方程，采用Fortran語言進行編譯
標簽： Fortran
上傳時間： 2018-06-07
上傳用戶：zhuchenmp
function [alpha,N,U]=youxianchafen2(r1,r2,up,under,num,deta) %[alpha,N,U]=youxianchafen2(a,r1,r2,up,under,num,deta) %該函數用有限差分法求解有兩種介質的正方形區域的二維拉普拉斯方程的數值解 %函數返回迭代因子、迭代次數以及迭代完成后所求區域內網格節點處的值 %a為正方形求解區域的邊長 %r1,r2分別表示兩種介質的電導率 %up,under分別為上下邊界值 %num表示將區域每邊的網格剖分個數 %deta為迭代過程中所允許的相對誤差限 n=num+1; %每邊節點數 U(n,n)=0; %節點處數值矩陣 N=0; %迭代次數初值 alpha=2/(1+sin(pi/num));%超松弛迭代因子 k=r1/r2; %兩介質電導率之比 U(1,1:n)=up; %求解區域上邊界第一類邊界條件 U(n,1:n)=under; %求解區域下邊界第一類邊界條件 U(2:num,1)=0;U(2:num,n)=0; for i=2:num U(i,2:num)=up-(up-under)/num*(i-1);%采用線性賦值對上下邊界之間的節點賦迭代初值 end G=1; while G>0 %迭代條件：不滿足相對誤差限要求的節點數目G不為零 Un=U; %完成第n次迭代后所有節點處的值 G=0; %每完成一次迭代將不滿足相對誤差限要求的節點數目歸零 for j=1:n for i=2:num U1=U(i,j); %第n次迭代時網格節點處的值 if j==1 %第n+1次迭代左邊界第二類邊界條件 U(i,j)=1/4*(2*U(i,j+1)+U(i-1,j)+U(i+1,j)); end if (j>1)&&(j U2=1/4*(U(i,j+1)+ U(i-1,j)+ U(i,j-1)+ U(i+1,j)); U(i,j)=U1+alpha*(U2-U1); %引入超松弛迭代因子后的網格節點處的值 end if i==n+1-j %第n+1次迭代兩介質分界面（與網格對角線重合）第二類邊界條件 U(i,j)=1/4*(2/(1+k)*(U(i,j+1)+U(i+1,j))+2*k/(1+k)*(U(i-1,j)+U(i,j-1))); end if j==n %第n+1次迭代右邊界第二類邊界條件 U(i,n)=1/4*(2*U(i,j-1)+U(i-1,j)+U(i+1,j)); end end end N=N+1 %顯示迭代次數 Un1=U; %完成第n+1次迭代后所有節點處的值 err=abs((Un1-Un)./Un1);%第n+1次迭代與第n次迭代所有節點值的相對誤差 err(1,1:n)=0; %上邊界節點相對誤差置零 err(n,1:n)=0; %下邊界節點相對誤差置零 G=sum(sum(err>deta))%顯示每次迭代后不滿足相對誤差限要求的節點數目G end
標簽： 有限差分
上傳時間： 2018-07-13
上傳用戶：Kemin
51里，從中斷開始，很多初學者就會感到很迷茫，不知道自己在學什么，特上傳此資料。單片機在執行程序的過程中,暫時中斷執行當前的程序,轉而去執行其他的應急處理程序,稱為中斷。例如。你正在餐廳吃飯，有朋友喊你去接電話，于是你就收到了來自朋友的一個中斷(可以叫做外部中斷)，當你準備去接電話時，突然肚子疼，需要上廁所（內部中斷），這又是一個中斷，我們把引起中斷的事件叫中斷源（例如接電話、上廁所等，外部引起的叫外部中斷，內部引起的叫內部中斷），產生中斷就要去處理它，這稱為中斷的響應。在接電話、上廁所這些中斷源中，顯然上廁所更需要立刻處理，這就是中斷的優先級。 51單片機的中斷系統十分重要，分為外部中斷和定時器中斷。中斷發生CPU在處理某一事件A時，發生了另一事件B請求CPU迅速去處理中斷響應和中斷服務CPU暫時中斷當前的工作，轉去處理事件B中斷返回待CPU將事件B處理完畢后，再回到原來事件A被中斷的地方繼續處理事件A這一過程稱為中斷單片機有10個寄存器主要與中斷程序的書寫控制有關中斷允許控制寄存器IE定時器控制寄存器TCON串口控制寄存器SCON中斷優先控制寄存器IP定時器工作方式控制寄存器TMOD定時器初值賦予寄存器（TH0/TH1,TL0/TL1）
上傳時間： 2022-02-16
上傳用戶：qdxqdxqdxqdx
