% Sinuskurve und ihre Schnittpunkte %

% --- Parameter --- %

B=9;                         % Gerade für Schnittpunkt
x=[0:1:100];                 % Definition des cycles von 0 bis 100%
P=B+0.1*B;                 % Maximum der Sinuskurve [°]
N=-3;                      % Minimum der Sinuskurve [°]
off=(P+N)/2             % Offset
a=(P-N)/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)


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

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


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

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