% Sinuskurve und ihre Schnittpunkte %

% --- Parameter --- %

B=9;                         % Gerade für Schnittpunkt
x=[0:1:100];                 % Definition des cycles von 0 bis 100%
max=B+0.1*B;                 % Maximum der Sinuskurve [°]
min=-3;                      % Minimum der Sinuskurve [°]
off=(max+min)/2;             % Offset
a=(max-min)/2;               % Amplitude
phase=asin(off/a)*(50/pi);   % Phase
f=1;                         % Frequenz

% ---- Sinusfunktion ---- %

y_1=a*sin((x*f-phase)*pi/50)+off;     


% ---- Plot ---- %

figure(1)
plot(x,y_1)

title('Sinus 1', 'FontSize',14)
xlabel('cycle[%]','FontSize',12)
ylabel ('[°]','FontSize',12)


% --- Ausgabe der enthaltenen Variablen als double --- %
a=double (a)
phase=double (phase)
off=double (off)
B=double(B)



% ----Berechnung Schnittpunkte mit der x-Achse ---- %

[x,y]=solve('y=a*sin((x*1-phase)*pi/50)+ off','y=0'); %Berechnung der SP
x_SP=x(1)  %Ausgabe x-Wert von SP1 im Command Window
x_SP= double(x_SP) % Ausgabe als double-Wert


% ----Berechnung Schnittpunkte mit y=B (9) ---- %

[x,y]=solve('y=a*sin((x*1-phase)*pi/50)+ off','y=B'); %Berechnung der x-Werte von y=B
x_B_1=x(1) 
x_B_1=double(x_B_1)
x_B_2=x(2) 
x_B_2=double(x_B_2)