#include #include #include int main( int c, char *y[] ){ double amplituda,frekvence,delka,deltaT,t,vzorky,um,f,g,dsp,zp1,zp2,x,theta; int a; amplituda=atof(y[1]); frekvence=atof(y[2]); delka=atof(y[3]); deltaT=0.000125; if(frekvence<0)return 1; if(delka<0)return 1; //printf("%f %f %f\n",amplituda,frekvence,delka); //printf("%f\n",deltaT); vzorky=delka/(deltaT*1000); //printf("%f\n",vzorky); theta=2*M_PI*(frekvence/8000); f=sin(theta); g=2*cos(theta); zp1=0; zp2=0; for(a=0;a<=vzorky;a++){ x=zp1; t=a*deltaT; dsp=um*f+zp1*f+zp2*g; um=amplituda*sin(2.0*M_PI*frekvence*t); dsp=um+zp1*f+zp2*g; printf("%lf %lf %lf\n " ,t,um,dsp); zp1=um; if(a<=1)continue; zp2=x; } }