Bode Plot - 1045

Bode Plot

Time Limit: 1000MS		Memory Limit: 10000K
Total Submissions: 13495		Accepted: 8523

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 cos wt 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 9
0.01
0.031623
0.1
0.31623
1.0
3.1623
10.0
31.623
100.0

Sample Output

0.010
0.032
0.100
0.302
0.707
0.953
0.995
1.000
1.000

分析：

此题的关键在于进行微积分的公式推导，其中涉及复合函数求导和三角函数的计算等问题，具体推导过程如下：

最后其实就是:VR = CRW*VS*sqrt(1/pow(CRW) + 1)

代码实现如下：

 
 

  
  
/***********************************************************************
    Copyright (c) 2015,wangzp
    All rights no reserved.
  
    Name: 《Bode Plot》In PEKING UNIVERSITY ACM
    ID:   PROBLEM 1045
    问题简述:  根据已知条件求电压值，关键在于公式的推导。
 
    Date: Aug 07, 2015 
 
***********************************************************************/
#include 
    
    
     
     
#include 
     
     
      
      

int main(void)
{
	double Vr,Vs,R,C,w;
	double tmp;
	int n,m;
	scanf("%lf %lf %lf %d",&Vs,&R,&C,&n);
	for (m = 0;m < n;m++)
	{
		scanf("%lf",&w);
		tmp = C * R * w;
		Vr = ( C * R * w * Vs) / sqrt(1 +  pow(tmp,2));
		printf("%.3lf\n",Vr);
	}
	return 0;
}