shamantou blog site

在bbs上看见有人发帖求"求一个最大值为800，周期为20ms的正弦波在前半个周期内跟最大值为1000，周期为0.1ms的三角波的交点"，好像很容易，就当炼手，写的时候才发现原来这个东西要用到数值分析方程求根的递归算法：

// 三角波周期0.1ms,幅值1000
// SIN波周期20ms,幅值800
// 求交点坐标
// 使用二分法迭代求解根
// 2006.05.21

#include "stdafx.h"
#include "conio.h"
#include "math.h"

#define PI 3.14159265358979
#define ep 0.00000000000001 //解的误差设置

double w=2*PI/20;//定义SIN角速度

inline double calc_sin(double t)//计算t时刻sin值
{
return 800*sin(w*t);
}

inline double calc_tri(double t)//计算t时刻三角波值
{
double t1=t-(int)(t*10)/10.0;
if(t1<0.05)
return t1*20000;
else
return (0.1-t1)*20000;
}

inline double diedai(double t1, double t2)//求解在t时刻附近的 Sin(x)-Tri(x)的近似解
{
double t;
if(calc_sin(t1)-calc_tri(t1)==0) return calc_sin(t1)-calc_tri(t1);
if(calc_sin(t2)-calc_tri(t2)==0) return calc_sin(t2)-calc_tri(t2);

do{
t=(t1+t2)/2;
double temp1=calc_sin(t)-calc_tri(t);
double temp2=calc_sin(t1)-calc_tri(t1);
double temp=temp1*temp2;

if(temp1==0)
return t;
if(temp<0)
t2=t;
if(temp>0)
t1=t;
}while (t2-t1>ep);
return (t1+t2)/2;
}

int main(int argc, char* argv[])
{
printf("\tSin(X)-sin波\tTri(X)-三角波 \t交点计算\n\n");
double fx;
for(double t=0;t<=10;t+=0.05)
{
//调用迭代函数
fx=diedai(t,t+.05);
printf("Sin(%.11f)= \tTri(%.11f)= \t%.11f\n",fx,fx,calc_sin(fx));
}

getch();
return 0;
}