最近无聊,研究上C语言的一些程序,感觉波形生成还是挺有用的
这里主要生成了正弦,方波,三角波 ,锯齿波,指数波形等
首先看看效果,
这里把生成的数据用python画了出来,程序在最后
正弦,调用参数
wave_gen(1000,10,100,sin);
波形长度1000个(1s),频率10hz,幅值100,类型为sin
方波,调用参数
wave_gen(1000,5,0.1,squ);
波形长度1000,频率5hz,占空比10%,类型为方波
锯齿波,调用参数
wave_gen(1000,5,0.1,saw);
波形长度1000(ms),频率5hz,斜率0.1(相对于时间ms而言),波形类型为锯齿波
三角波,调用参数
wave_gen(1000,5,10,tri);
波形长度1000ms,频率5hz,斜率10(相对于时间ms而言),波形类型为三角波,其实也就是特殊的锯齿波
x/2的积分,即x^2指数函数
调用参数wave_gen(1000,5,2,inte);
波形长度1000,频率5hz,2次幂,积分(指数)波形
剩下的一些波形还有点问题,就不演示了
贴下代码
#include "stdio.h"
#include "stdlib.h"
#include "string.h"
#include "math.h"
#define PI 3.14159265
float *ydata,*xdata;
#define sin 0 //正弦
#define squ 1 //方波
#define saw 2 //锯齿波
#define tri 3 //三角波
#define inte 4 //积分 x^2
#define diff 5
#define gauss 6 // 高斯白噪声
void wave_gen(int len,float fre,float rate,char type){//生成波形长度ms,波形频率 ,波形参数,波形类型
int i,wave_cnt=0,wave_sw_value=1000/fre;
ydata = (float *)malloc(sizeof(float)*len);
xdata = (float *)malloc(sizeof(float)*len);
for(i=0;i<len;i++){
if(wave_cnt>=wave_sw_value) wave_cnt=0;
switch(type){
case sin:
ydata[i] = rate*sinf(2*wave_cnt*PI/wave_sw_value);
break;
case squ:
ydata[i] = wave_cnt < rate*wave_sw_value ? 1:0;
break;
case saw:
ydata[i] = rate*wave_cnt;
break;
case tri:
ydata[i] = rate*(wave_cnt>wave_sw_value/2 ? wave_sw_value-wave_cnt: wave_cnt);
break;
case inte:
ydata[i] = (pow(wave_cnt,rate));
break;
case diff:
ydata[i] = 2*wave_cnt;
break;
case gauss:
ydata[i] = (20.0 / (10 * sqrt(2.0 * 3.14159))) * exp(-((wave_cnt/1.0 - 10) * (wave_cnt/1.0 - 10)) / (2.0 * 10 * 10));
break;
}
xdata[i] = i;
wave_cnt++;
}
}
int main(){
int i=0;
wave_gen(1000,5,2,inte);
for(i=0;i<1000;i++){
printf("%.2f,",ydata[i]);
}
return 0;
}