blob: c81b4d6c2f10aca2f00b724a9bbbc9860b74f90a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
close();
clear;
clc;
//resistance 'R', rate of electrical energy 'rate', time of operation 't', time period 'T'
T = 0.01; //s
R = 20; //ohm
rate = 0.06; //$ per kWh
time = 24; //hours
//from graph, voltage as function of time 'v'
function V = v(t)
V = 10^4 * t; //V
endfunction
//current as a function of time 'i'
function I = i(t)
I = v(t)/R; //A
endfunction
//power as afunction of time 'p'
function power = p(t)
power = v(t)*i(t); //W
endfunction
//power 'P'
P = 1/T * intg(0,T,p)/1000; //kW
energy = P*time; //kWh
//cost of electrical energy 'cost'
cost = rate*energy; //$
mprintf("Electrical energy costs $ %0.2f for %d hours",cost,time);
|
