poj 1045 Bode Plot

 1#include <stdio.h>
 2#include <math.h>
 3
 4int main()
 5{
 6    float vs, r, c;
 7    float vr,vri,vrr,tao;
 8    float omega;
 9    int n,i;
10
11    /**//*    
12        Vc+Vr=Vs
13        I=C*s*Vc
14        Vr=I*R
15        Vr=s*R*C*(Vs-Vr)
16        Vr=(sRC/1+sRC)*Vs=(jwRC+(wRC)^2)*Vs/(1+(wRC)^2)
17    */    
18    scanf("%f%f%f%d",&vs,&r,&c,&n);
19
20    for(i=0;i<n;i++)
21        {
22        scanf("%f",&omega);
23        tao=omega * r * c;
24               /**//*        
25                vri=tao * vs / (1 + tao * tao);
26        vrr=tao * tao * vs/ (1 + tao * tao);
27        vr=sqrt(vri*vri + vrr*vrr);
28                */
29        vr=tao * vs / sqrt(1+tao *tao);
30
31        printf("%0.3f\n",vr);
32    }
33}
34

用拉氏变换解很简单。

通信没白学啊。
#include <stdio.h>
#include <math.h>

int main()
{
    float vs, r, c;
    float vr,vri,vrr,tao;
    float omega;
    int n,i;

    /*    
        Vc+Vr=Vs
        I=C*s*Vc
        Vr=I*R
        Vr=s*R*C*(Vs-Vr)
        Vr=(sRC/1+sRC)*Vs=(jwRC+(wRC)^2)*Vs/(1+(wRC)^2)
    */    
    scanf("%f%f%f%d",&vs,&r,&c,&n);

    for(i=0;i<n;i++)
    {
        scanf("%f",&omega);
        tao=omega * r * c;
       /*      
        vri=tao * vs / (1 + tao * tao);
        vrr=tao * tao * vs/ (1 + tao * tao);
        vr=sqrt(vri*vri + vrr*vrr);
        */
        vr=tao * vs / sqrt(1+tao *tao);

        printf("%0.3f\n",vr);
    }
}