sin_flt.cpp

一些在CCS上可以运行的初级DSP汇编程序

#include <math.h>
#include<stdlib.h>
#include<iostream.h>
#include<stdio.h>

void sine(double f,double ph,double fs,double x[],int n)
{
	double PI;
	int k;
	PI=4*atan(1.0);
	f=2*PI*f/fs;
	ph=ph*PI/180;
	FILE *fp;
	fp=fopen("iirsin.dat","w");
	if(fp==NULL)
	{
		printf("Cannot open the file!\n");
		exit(0);
	}
	for(k=0;k<n;k++)
	{
		x[k]=sin(k*f+ph);
		fprintf(fp,"%lf\n",x[k]);
		if(-1<x[k]&x[k]<0)  
			x[k]=2+x[k];
	}
	fclose(fp);
}

main()
{
	int n,result,i;
	double f,fs,ph;
	//double (*x) ;
	//x= new double[n];
	//n=10;
	FILE *fp1;
	ph=0;
	printf("Please Input N:\n");
	result=scanf("%d",&n);
	double (*x) ;
	x= new double[n];
	int (*X);
	X= new int[n];

	printf("Please Input f:\n");
	result=scanf("%lf",&f);

	printf("Please Input fs:\n");
	result=scanf("%lf",&fs);

	sine(f,ph,fs,x,n);

	for(i=0;i<n;i++)
	{
		X[i]=(int)(32768*x[i]);
		if(x[i]==1)	
			X[i]=32767;
		if(x[i]==-1)
			X[i]=32768;
	}

	fp1=fopen("IN.dat","w");

	if(fp1==NULL)
	{
		printf("Cannot open file !");
		exit(0);
	}

	for(i=0;i<n;i++)
	{
		fprintf(fp1,"%x\n",X[i]);
	}
	fclose(fp1);
	/*fp2=fopen("in.dat","w");

	if(fp2==NULL)
	{
		printf("Cannot open file !");
		exit(0);
	}
	double a;
	a=0;
	for(i=0;i<n;i++)
	{
		fprintf(fp2,"%x\n%x\n",X[i],a);
	}
	fclose(fp2);*/
	delete[] x;
	delete[] X;
	return 1;
}