zbibo.c

来自「傅里叶变换」· C语言 代码 · 共 64 行

#include <math.h>
#include <stdio.h>
#include "lcd_86.h"
#include "mmdrv.h"
#include "mcbsp.h"
#include "cpu_reg.h"    
#include "adda.h"
int x[10] = {3,1,-3,-4,12,3,12,-4,-3,1};

main()
{
	int b[100],b1[99],r;
	int n1,n,i,j,nn,k,*buf;
	char bchar[6];
	
//	init_lcd();
	buf = x;		
	for(k=0;k<2;k++)
	{
		n = *buf++;
		n1 = n+1;
  		LCD_clear(0);
		cur_row=0;        
		cur_col=6;        
    		LCD_pr_chars("polynomial of b[i]",18);
		for(i=0;i<n1;i++)
		{
			b[i] = *buf++;
    			Int_to_String(b[i],bchar);
    			cur_row=1+i;        
			cur_col=1;        
    			LCD_pr_chars(bchar,5);
		}
		j=1;
		do
		{
			r = b[0]/b[n1-1];
			if (abs(r)>=1)
			{
				cur_row=1+n1;        
				cur_col=1;        
    			LCD_pr_chars("filter is not stable",20);
    			break;
    		}
    		nn = n1 - 1;
    		if (nn==1)
    		{
    			cur_row=1+n1;        
				cur_col=1;        
    			LCD_pr_chars("filter is stable",16);
    			break;
    		}
    		for(i=0;i<n1-1;i++)
    			b1[i] = b[n1-1]*b[i+1]-b[0]*b[n1-1-i];
    		b[n1-1]=0;
    		for(i=0;i<n1-1;i++)
    			b[i] = b1[i];
    		n1=n1-1;
    		j = j+1;
    	}while(1);
    }
}