用C++ 生成的正弦波数据文件

#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; }