3943.水温控制系统程序附件

1、附：系统程序主函数文件#include lpc23xx.h/* lpc23xx/24xx peripheral registers*/#include type.h#include irq.h#include target.h#include timer.h#include rtc.h#include define.h#include coefficient.h#include keyled_define.h#include lcd_menu_resource.h#include math.h #include string.h#include stdio.h#define samplecyc

2、lenum 5#define pi 3.14159265#define sqrt_2 1.41421357#define pointsnumber 64extern _armabi float roundf(float /*x*/);double caculate_frequency(int piont,double frequency_0);void delay_for_ad(void);void adgpio_initial(void);enum system_state systemstate; /系统状态，分为测频阶段和参量计算两个状态int ad0_value512,ad1_valu

3、e512; /采样值word totalpointsnumber;dword adc_timerinterval;byte timescount_for_frequency; / 测频状态时记录采集周期数byte newvaluecreated_flag; /新菜了一个周期波形的标志byte frequnecystateenter_flag; /进入频率测试状态的标志byte elsestateenter_flag; /点入电参量计算的标志dword file_data_buf_byte_count=0;byte usbconnect_flag=0; /u盘连接标志，当u盘连接时，标志置1in

4、t upointsnumber*3+1; int ipointsnumber*3+1;float a_for_u11;float b_for_u11;/傅氏分解的实部、虚部float a_for_i11;float b_for_i11;/傅氏分解的实部、虚部struct voltage_coe u_coe_temp10,u_coe_value, i_coe_temp10,i_coe_value;/各谐波的有效值,记录十次取平均值float u_value=0,i_value=0,pf_value=0,power_value=0,q_value=0,s_value=0; /分别为电压、电流、有功

5、、无功和视在功率float umax_value=0,imax_value=0,powermax_value=0,qmax_value=0,smax_value=0; /分别为电压、电流、有功、无功和视在功率最大值float umin_value=0,imin_value=0,powermin_value=0,qmin_value=0,smin_value=0; /分别为电压、电流、有功、无功和视在功率最小值double frequency_value; /频率float u_temp10,i_temp10,power_temp10,q_temp10,s_temp10; /分别为电压、电流、有

6、功、无功和视在功率byte valuetempcount=0;int main(void) word i=0;byte j;byte keynumber; word temp; scs|=1; initkeyledrelay_pin(); lcd_init();/初始化clearscreen(); /清屏menu_pageid=pageid_mainpage;key_databuf=0x3f;enter_flag=1;esc_flag=0;keylcd_deal(nokey);backpara(pageid_vrefadj);/* rtc时钟初始化 */ systemtime.rtc_sec=

7、0; systemtime.rtc_min=0; systemtime.rtc_hour=0; systemtime.rtc_mday=1; systemtime.rtc_mon=6; systemtime.rtc_year=2008;backpara(pageid_timeset); rtcinit(); /* initialize timer0,for sample pf */ init_timer ( 0, fpclk); /起动定时器0，作为测功率因数用 enable_timer(0); adgpio_initial(); usb_init(); frequency_value=50;

8、 timescount_for_frequency=0; systemstate=undetect_state; /初始设置为非检测状态 / frequnecystateenter_flag=1; elsestateenter_flag=1; while(1) if(elsestateenter_flag=1) elsestateenter_flag=0; timer1_counter_for_points=0; adc_timerinterval= ( fpclk/(totalpointsnumber_for_else*frequency_value)-1); init_timer ( 1,

9、 adc_timerinterval ); /将在tmer.c文件的timer0handler (void) _irq 函数中处理adc 的采样 enable_timer(1); totalpointsnumber=totalpointsnumber_for_else; keynumber=keyscan(); /扫描按键 if(keynumber!=nokey) keylcd_deal(keynumber); if(newvaluecreated_flag=1) newvaluecreated_flag=0; timescount_for_frequency+; if(timescount_

10、for_frequency=3) timescount_for_frequency=0; timer1_counter_for_points=0; caculate_value(); keylcd_deal(nokey); /存储记录数据file_data_buffile_data_buf_byte_count+=systemtime.rtc_mon;file_data_buffile_data_buf_byte_count+=systemtime.rtc_mday;file_data_buffile_data_buf_byte_count+=systemtime.rtc_hour;file_

11、data_buffile_data_buf_byte_count+=systemtime.rtc_min;file_data_buffile_data_buf_byte_count+=(int)(u_value*10)8;file_data_buffile_data_buf_byte_count+=(int)(u_value*10)&0xff);file_data_buffile_data_buf_byte_count+=(int)(u_value*10)8;file_data_buffile_data_buf_byte_count+=(int)(i_value*100)&0xff);file

12、_data_buffile_data_buf_byte_count+=(int)(i_value*100)8;file_data_buffile_data_buf_byte_count+=(int)(power_value)&0xff);file_data_buffile_data_buf_byte_count+=(int)(power_value)8;file_data_buffile_data_buf_byte_count+=(int)(q_value)&0xff);file_data_buffile_data_buf_byte_count+=(int)(q_value)8;file_da

13、ta_buffile_data_buf_byte_count+=(int)(s_value)&0xff);file_data_buffile_data_buf_byte_count+=(int)(s_value)8; if(file_data_buf_byte_count=file_data_buf_len)file_data_buf_byte_count=0; if(usbconnect_flag =1 ) /检测到u盘插入,对u盘进行处理 usbconnect_flag=0; if(parachanged_flag=1) /如果正在修改参数，先将其参数撤销，在更改到usb通讯界面 keyl

14、cd_deal(key_esc); clearscreen(); disp_string(0,2,16,menu_usbconnect); usb_deal(); keylcd_deal(nokey); adc_timerinterval= ( fpclk/(totalpointsnumber_for_else*frequency_value)-1); init_timer ( 1, adc_timerinterval ); /将在tmer.c文件的timer0handler (void) _irq 函数中处理adc 的采样 enable_timer(1); / elsestateenter_

15、flag=1; else reset_timer(1); enable_timer(1); void caculate_value(void)word i,j;float temp;double ratio; u_tempvaluetempcount=0;i_tempvaluetempcount=0; for(j=0;j10;j+) a_for_uj=0; b_for_uj=0; a_for_ij=0; b_for_ij=0; power_tempvaluetempcount=0;q_tempvaluetempcount=0;s_tempvaluetempcount=0; for(i=0;i(

16、pointsnumber*3+1);i+) u_tempvaluetempcount=u_tempvaluetempcount+ad0_valuei*ad0_valuei*coefficienti;i_tempvaluetempcount=i_tempvaluetempcount+ad1_valuei*ad1_valuei*coefficienti; for(j=0;j10;j+) a_for_uj=a_for_uj+sqrt_2*ad0_valuei*cos(j*i*2*pi/pointsnumber)*coefficienti;b_for_uj=b_for_uj+sqrt_2*ad0_va

17、luei*sin(j*i*2*pi/pointsnumber)*coefficienti; a_for_ij=a_for_ij+sqrt_2*ad1_valuei*cos(j*i*2*pi/pointsnumber)*coefficienti;b_for_ij=b_for_ij+sqrt_2*ad1_valuei*sin(j*i*2*pi/pointsnumber)*coefficienti; power_tempvaluetempcount=power_tempvaluetempcount+ad0_valuei*ad1_valuei*coefficienti; if(i(samplecycl

18、enum-2) /采集完十次，进行求平均值 valuetempcount=0; u_value=0;i_value=0; /各数据在求平均值之前清0 power_value=0;q_value=0; s_value=0; u_coe_value.degree1=0;u_coe_value.degree2=0;u_coe_value.degree3=0;u_coe_value.degree4=0;u_coe_value.degree5=0;u_coe_value.degree6=0;u_coe_value.degree7=0;u_coe_value.degree8=0;u_coe_value.d

19、egree9=0;u_coe_value.degree10=0; i_coe_value.degree1=0;i_coe_value.degree2=0;i_coe_value.degree3=0;i_coe_value.degree4=0;i_coe_value.degree5=0;i_coe_value.degree6=0;i_coe_value.degree7=0;i_coe_value.degree8=0;i_coe_value.degree9=0;i_coe_value.degree10=0; for(i=0;isamplecyclenum;i+) u_value=u_value+u

20、_tempi; i_value=i_value+i_tempi; power_value=power_value+power_tempi; q_value=q_value+q_tempi; s_value=s_value+s_tempi; u_coe_value.degree1=u_coe_value.degree1+u_coe_tempi.degree1; u_coe_value.degree2=u_coe_value.degree2+u_coe_tempi.degree2; u_coe_value.degree3=u_coe_value.degree3+u_coe_tempi.degree

21、3; u_coe_value.degree4=u_coe_value.degree4+u_coe_tempi.degree4; u_coe_value.degree5=u_coe_value.degree5+u_coe_tempi.degree5; u_coe_value.degree6=u_coe_value.degree6+u_coe_tempi.degree6; u_coe_value.degree7=u_coe_value.degree7+u_coe_tempi.degree7; u_coe_value.degree8=u_coe_value.degree8+u_coe_tempi.d

22、egree8; u_coe_value.degree9=u_coe_value.degree9+u_coe_tempi.degree9; u_coe_value.degree10=u_coe_value.degree10+u_coe_tempi.degree10; i_coe_value.degree1=i_coe_value.degree1+i_coe_tempi.degree1; i_coe_value.degree2=i_coe_value.degree2+i_coe_tempi.degree2; i_coe_value.degree3=i_coe_value.degree3+i_coe

23、_tempi.degree3; i_coe_value.degree4=i_coe_value.degree4+i_coe_tempi.degree4; i_coe_value.degree5=i_coe_value.degree5+i_coe_tempi.degree5; i_coe_value.degree6=i_coe_value.degree6+i_coe_tempi.degree6; i_coe_value.degree7=i_coe_value.degree7+i_coe_tempi.degree7; i_coe_value.degree8=i_coe_value.degree8+

24、i_coe_tempi.degree8; i_coe_value.degree9=i_coe_value.degree9+i_coe_tempi.degree9; i_coe_value.degree10=i_coe_value.degree10+i_coe_tempi.degree10; u_value=(u_value/samplecyclenum)*100; /得到平局值并乘以pt ct变比 i_value=(i_value/samplecyclenum)*10; power_value=(power_value/samplecyclenum)*100*10; q_value=(q_va

25、lue/samplecyclenum)*100*10; s_value=u_value*i_value; pf_value=power_value/s_value; u_coe_value.degree1=(u_coe_value.degree1/samplecyclenum)*100; u_coe_value.degree2=(u_coe_value.degree2/samplecyclenum)*100; u_coe_value.degree3=(u_coe_value.degree3/samplecyclenum)*100; u_coe_value.degree4=(u_coe_valu

26、e.degree4/samplecyclenum)*100; u_coe_value.degree5=(u_coe_value.degree5/samplecyclenum)*100; u_coe_value.degree6=(u_coe_value.degree6/samplecyclenum)*100; u_coe_value.degree7=(u_coe_value.degree7/samplecyclenum)*100; u_coe_value.degree8=(u_coe_value.degree8/samplecyclenum)*100; u_coe_value.degree9=(

27、u_coe_value.degree9/samplecyclenum)*100; u_coe_value.degree10=(u_coe_value.degree10/samplecyclenum)*100; i_coe_value.degree1=(i_coe_value.degree1/samplecyclenum)*10; i_coe_value.degree2=(i_coe_value.degree2/samplecyclenum)*10; i_coe_value.degree3=(i_coe_value.degree3/samplecyclenum)*10; i_coe_value.

28、degree4=(i_coe_value.degree4/samplecyclenum)*10; i_coe_value.degree5=(i_coe_value.degree5/samplecyclenum)*10; i_coe_value.degree6=(i_coe_value.degree6/samplecyclenum)*10; i_coe_value.degree7=(i_coe_value.degree7/samplecyclenum)*10; i_coe_value.degree8=(i_coe_value.degree8/samplecyclenum)*10; i_coe_v

29、alue.degree9=(i_coe_value.degree9/samplecyclenum)*10; i_coe_value.degree10=(i_coe_value.degree10/samplecyclenum)*10; if(systemstate=detect_state) /检测状态下记录更新最大最小值 if(umax_value u_value) umin_value=u_value; if(imax_value i_value) imin_value=i_value; if(powermax_value power_value) powermin_value=power_

30、value; if(qmax_value q_value) qmin_value=q_value; if(smax_value s_value) smin_value=s_value; else valuetempcount+;double caculate_frequency(int piont,double frequency_0)word i; double angle1,angle2,angle2_difference,frequency,frequency_difference;double r;/r为修正系数double a1,b1,a2,b2; for(i=0;itotalpoi

31、ntsnumber_for_caculate;i+) a1=a1+2*pionti*cos(i*2*pi/totalpointsnumber_for_caculate); b1=b1+2*pionti*sin(i*2*pi/totalpointsnumber_for_caculate); a2=a2+2*pionti+totalpointsnumber_for_caculate*cos(i*2*pi/totalpointsnumber_for_caculate); b2=b2+2*pionti+totalpointsnumber_for_caculate*sin(i*2*pi/totalpointsnumber_for_caculate); a1=a1/totalpointsnumber_for_caculate; b