matlab에서 전기 에너지 플로팅

Question

나는 matlab를 사용하여 회로의 에너지 (W)를 표시하려고합니다. 방정식은 다음과 같습니다

는, 전류 전압을 플로팅 전원이 어려운 일이 아니었다 (나는 내가 쓴 생각 올바른 코드를 알려주세요. 전자 당신은 잘못 생각하는 경우) :!

%Variables 
t1=20e-3;t2=60e-3;N=1000;f=50; 
t=linspace(t1,t2,N); 
Vmax = 10; %Peak Voltage 
Imax = 2; %Peak Current 
f = 50; %Frequency in Hertz 
omega = 2 * pi * f; 
phi_default = -pi/4; %Constant for phase difference 

%Equations 
Vsrc = Vmax * sin(omega * t); %Equation for Voltage 
Isrc = Imax * sin((omega * t) + phi_default); %Equation for Current 
Psrc = Vsrc .* Isrc; %Equation for Power 

%Plotting the result 
figure(1); %Window 1 
plot(t,Vsrc, 'b', t, Isrc, 'r', t, Psrc, 'g'); %Plotting Voltage, Current and Power 

하는 것은 이제 문제는 내가 MATLAB에서 에너지를 그릴 수있는 방법을 알아낼 수 있다는 것입니다. 나는 Δt (Delta t)가 어쩌면 같은해야 뭔가를 diff(t) 기능과 전체 에너지를 사용하여 정의해야한다고 생각 : T 기간 인 경우

그리고 평균 전력 Pavr이 에 의해 계산된다.

질문 : 그럼, 의 총 순간 에너지를 범위에서 위상차의 함수로 어떻게 그릴 수 있습니까? for loop은 각 위상 값을 계산해야합니다.

또한 평균 전력을 같은 기간에 표시하고 싶습니다.

감사합니다.

Answer 1

에너지는 일정 기간 동안 전력을 통합하여 계산됩니다. 여기에서 불연속 값을 다루므로 합계로 근사값을 구해야합니다. 귀하의 경우 시간에 따른 순시 에너지

Winst = (t2 - t1)/N * cumsum(Psrc);

cumsum 당신에게 Psrc(t2 - t1)/N 및 시간 간격이고 벡터의 모든 요소의 누적 합을 제공한다.

Φ의 함수로 에너지를 계산하려면 고려할 기간을 고쳐야합니다. 당신은 주어진 파이 먼저 에너지를 계산하고 다시 일정 시간 간격 (T1)를 통해 통합 ... T2 :

phi = -pi/2:phi_step:0 
Wphi = zeros(1, length(phi)); 
for k = 1:length(phi) 
    Vsrc = Vmax * sin(omega * t); 
    Isrc = Imax * sin((omega * t) + phi(k)); 
    Psrc = Vsrc .* Isrc; 
    Wphi(k) = (t2 - t1)/N * sum(Psrc); 
end 

plot(phi, Wphi); 

평균 전력이 Psrc의 평균 값 단순히 :

Pavg = sum(Psrc)/length(Psrc)