POJ 1045 Bode Plot

Bode Plot

Time Limit: 1000MS Memory Limit: 10000KTotal Submissions: 14537 Accepted: 9069

Description

Consider the AC circuit below. We will assume that the circuit is in steady-state. Thus, the voltage at nodes 1 and 2 are given by v₁ = V_S coswt and v₂ = V_Rcos (wt + q ) where V_S is the voltage of the source, w is the frequency (in radians per second), and t is time. V_R is the magnitude of the voltage drop across the resistor, and q is its phase. 

You are to write a program to determine V_R for different values of w. You will need two laws of electricity to solve this problem. The first is Ohm's Law, which states v₂ = iR where i is the current in the circuit, oriented clockwise. The second is i = C d/dt (v₁-v₂) which relates the current to the voltage on either side of the capacitor. "d/dt"indicates the derivative with respect to t. 

Input

The input will consist of one or more lines. The first line contains three real numbers and a non-negative integer. The real numbers are V_S, R, and C, in that order. The integer, n, is the number of test cases. The following n lines of the input will have one real number per line. Each of these numbers is the angular frequency, w. 

Output

For each angular frequency in the input you are to output its corresponding V_R on a single line. Each V_R value output should be rounded to three digits after the decimal point.

Sample Input

1.0 1.0 1.0 90.010.0316230.10.316231.03.162310.031.623100.0

Sample Output

0.0100.0320.1000.3020.7070.9530.9951.0001.000

//这道题纯属公式题，推导出计算VR的公式即可。利用iR=CR d/dt (VScos(wt)-VRcos (wt + q ))=VRcos (wt + q ) 得到CRW(VRsin(wt+q)-VSsin(wt))=VRcos (wt + q )  .之后分别令t=0和wt+q=0产生两个式子，合并之后得VR=CRW*Vs/（sqrt(CRW*CRW)+1）

#include <iostream>#include <cstdio>#include <cmath>using namespace std;int main(int argc, char const *argv[]){int n;double a,b,c;scanf("%lf %lf %lf %d",&a,&b,&c,&n);while(n--){double w,value;scanf("%lf",&w);value=a*b*c*w/sqrt(1+pow(b*c*w,2));printf("%.3f\n",value);}return 0;}