Subversion Repositories svnkaklik

Rev

Rev 969 | Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | Download

set terminal png

set key off
set grid
set angles degrees
set pointsize 3

f(x)=I0 * (sin(x))**4 + c
g(x)=I0 * (sin(x) * cos(x))**2 + a


set xlabel "fi[°]"
set ylabel "I[-]"

fit f(x) "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 via I0, c
set output "polarizace.png"
plot "./data/Malusuv_zakon/polarizace_vertikalne_stupne.txt" using 1:2 with points, f(x)

fit g(x) "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 via I0, a
set output "polarizace_horizontalne.png"
plot "./data/Malusuv_zakon/polarizace_horizontalne_stupne.txt" using 1:2 with points, g(x)

set output "podelny_rez.png"
set xlabel "x[-]"
set ylabel "I[-]"
plot "podelny_rez.txt" using 1:2 with points

set output "pricny_rez.png"
set xlabel "x[-]"
set ylabel "I[-]"
plot "./data/rozlozeni_pole/pricny_rez.txt" using 1:2 with points

set samples 400
set angles radians
a=4
d=100
c=6

I0=3

h(x)=I0 * sin(c*x+d) + a

fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 via I0,a,d,c
set output "stojata_vlna.png"
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:2 with points, h(x)


I0=2

fit h(x) "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 via I0,a,d,c
set output "stojata_vlna_deska.png"
plot "./data/stojata_vlna/stojata_vlna_vzdalenost.txt" using 1:3 with points, h(x)



set output "hrana.png"
plot "./data/difrakce/hrana.txt" using 1:(abs($2)) with points

set output "pasek.png"
plot "./data/difrakce/pasek.txt" using 1:(abs($2)) with points

set key under
set xlabel "fi[°]"
set ylabel "I[-]"
set output "sterbina.png"


g(x)= I1*(sin(d1*pi/l1*sin(x))/(d1*pi/l1*sin(x)))**2

sinc(x) = sin(pi*x)/(pi*x)
i(x)=I2*(sinc(d2/l2*sin(x)))**2

set angles degrees

I1=10
l1=0.03
d1=0.04

I2=10
l2=0.03
d2=0.04


fit [5:60] g(x) "sterbina.txt" using 1:2 via I1,d1,l1
fit [5:60] i(x) "sterbina.txt" using 1:3 via I2,d2,l2
plot "sterbina.txt" using 1:2 with points title "40mm", "" using 1:(abs($3)) with points title "60mm", i(x), g(x)