原创力文档- 1、本文档共33页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验6
#includexxgc.h
double func1(double t)
{return (t=0)?sin(t):0;}
double func2(double t)
{return (t=0)?exp(-1*t/2):0;}
double reversal(double t,double *tao,double (*func)(),void *arg)
{return (*func)(*tao-t,arg);}
double f1CONVf2(double t)
{double sum=0,a=-5,dt=0.01;int N,k=0;
t=t-a*2;
N=t/dt;
for(;kN;k++)
sum+=func1(a+k*dt)*func2(a+(N-k)*dt);
return sum*dt;}
main()
{double tao=0,right;
struct WIN2 *p;
int yo;
initgd();
window2(dynamic convolution demo,
-3,2,15,-4,
t,f,RED,WHITE);
plotxy2(BLUE,func1);
p=getwin2();
right=p-x2;
yo=p-yo;
for(instKeyCtr();keyCtr(tao,-10,10,50);)
{clearviewport();
plotfb2(GREEN,func1);
plotxy2(BLUE,reversal,tao,func2);
p-x2=tao; p-yo+=80;
plotfb2(BLUE,f1CONVf2);
p-x2=right;
p-yo=yo;}
closegraph();}
实验7
1)
#includexxgc.h
double winfunc(double t,double tao[1])
{double u;
u=fmod(t,4);
if(u=0utao[0]) return u/tao[0];
else if(u=tao[0]u4) return (u-4)/(tao[0]-4);
else if(u0utao[0]-4) return u/(tao[0]-4);
else if(u=tao[0]-4u-4) return (4+u)/tao[0];
}
main()
{double a,b,tao,T=4;
dcomplex c[20];
int i,N=20;
initgd();
window3(amplitude charistic,
-1,-0.5,-0.2,21,0.5,0.2,
N,r,i,BLUE,WHITE);
for(tao=0,instKeyCtr();keyCtr(tao,0,4,20);)
{for(i=0;iN;i++)
{a=fcosCoef(i,T,winfunc,tao);
b=fsinCoef(i,T,winfunc,tao);
c[i].r=a;c[i].i=-b;}
c[0].r=c[0].r/2;
frame3(N,r,i,BLUE);xyz3(BLUE);
plotgri3(BLUE,RED,c,N);}
closegraph();
}
2)
#includexxgc.h
double winfunc(double t,double h[1])
{return (sin(t)+h[0]-1=0)?(sin(t)+h[0]-1):0;
}
main()
{double a,b,h,T=6.28;
dcomplex c[20];
int i,N=20;
initgd();
window3(amplitude charistic,
-1,-0.2,-0.2,21,0.2,0.2,
N,r,i,BLUE,WHITE);
for(h=0,instKeyCtr();keyCtr(h,0,1,20);)
{for(i=0;iN;i++)
{a=fcosCoef(i,T,winfunc,h);
b=fsinCoef(i,T,winfunc,h);
c[i].r=a;c[i].i=-b;}
c[0].r=c[0].r/2;
frame3(N,r,i,BLUE);xyz3(BLUE);
plotgri3(BLUE,RED,c,N);}
closegraph();
}
实验8
1)
#includexxgc.h
double etut(double t,double *tao)
{if(t=-1**tao/2t=*tao/2)
文档评论（0）