clc;clear all;
format long;

%%%% Vstupni udaje ziskane navrhem v EOD %%%%
r_0=38; %vnitrni polomer civky
dd=1; %prumer dratu civky
d=126; %delka civky
s=20; %sirka vinuti
uhel=360/180*pi; %uhel vinuti
az=1714.3825; %buzeni civky v amperzavitech při pracovní vzdálenosti 20mm
%rezerva=10;
ro=17.5e-6; %merny el. odpor materialu vynuti (med)

%%%% Vypocet proudu v dratech civky %%%%
vod=(d/dd); %pocet zavitu ve vodorovnem smeru
hor=(s/dd); %pocet zavitu v horizontalnim smeru
N = vod*hor % maximalni pocet zavitu vinuti
%M=N-rezerva % pocet zavitu brany s rezervou

I=(az/N) %proud v cívce
%I = (amm*(pi*(dd/2)^2))

%%%% Vypocet delky dratu %%%%
for i=1:hor
    x(i)=(uhel*(r_0+dd/2+(i-1)*dd))*(d/dd);
end
X=sum(x) %delka dratu (mm)

%%%% Vypocet odporu dratu bez uvazovani zahrivani %%%%
R=ro*(X/(pi*(dd/2)^2)) %odpor v ohmech

%%%% Potrebne napeti %%%%
U=R*I %napeti ve Voltech

%%%% Ztratovy vykon %%%%
ZU=R*I^2 %ve Wattech
    