% Orbitdarstellung % Voraussetzung: Belegung von 8 Kanälen (alles Beschl.Sensoren), CH1 - % CH8. % Ohne Analogsignale für Geschwindigkeit und % Wickeldurchmesser % Grafische Darstellung für VariFlex/ VariStep fs= VCAP_SAMPLERATE; % Erzeugter Datensatz in SigLab N= 1*2048; % Datenpunkte; Einstellung in SigLab: vna-Prog. Ts= 1/fs; % Def,: fs= Samplerate t= Ts*N; % Zeit für Datenaufnahme eines Frames Frames= length(VCAP_DATA)/N; % Anzahl der Frames im kompletten Datensatz %#### Eingabedaten ######### % Kalibrierfaktoren (damit werden die Signale multipliziert): x_V= 300; % Geschwindigkeit, 300 (m/min)/V x_D= 150; % Durchmesser, 150 (mm)/V x_x= 19607; % Beschleunigung, PCB-Aufnehmer (Beschleunigungssignale) x_B_K= 10000; % B&K-Aufnehmer (Beschleunigungssignale); Einstellung am Ladungsverstärker beachten x_Dornier= 0.25; % Induktiver Aufnehmer (Wegsignale), 0.25 (mm)/V x_Laser= 0.5; % Laser (Wegsignale), 0.5 (mm)/V % Hiermit wird ein Frame aus dem kompletten Datensatz ausgewählt k= 5; % Frame Nr aus dem Datensatz Dact= 812; % [mm], Wahl des Wickeldurchmessers für die grafische Darstellung Vact= 1000; % [m/min], Geschwindigkeit (Wird nur als Info für die Darstellung benutzt) fact= (1000/60)*Vact/(Dact*pi); % Mit Dact und Vact wird die Wickeldrehfrequenz berechnet % Maschinengeometrie: TW1= 620; % Durchmesser [mm] von Tragwalze 1 TW2= 620; % Durchmesser [mm] von Tragwalz2 2 DW= 812; % Durchmesser der Druckwalze S= 30; % Lichte Spalt [mm] zwischen den Tragwalzen TWH= 0.5*(TW1-TW2); % Höhenversatz zwischen beiden Tragwalzen [mm] %##### Ende der Eingabedaten ######### %#### Beachten #### % Bei gemischten Signalen an den Kanälen: Beschleunigungs- und Wegsignalen % müssen die Wegsignale durch zweifaches Differenzieren in Beschleunigungs- % signale umgerechnet werden. % #### Zuordnung der Kanäle ohne Korrektur % der Sensorrichtung: CH1= VCAP_DATA(:,1)*x_x; % Kanal 1, Beschleunigungssensor CH2= VCAP_DATA(:,2)*x_x; % Kanal 2, Beschleunigungssensor CH3= VCAP_DATA(:,3)*x_x; % Kanal 3, Beschleunigungssensor CH4= VCAP_DATA(:,4)*x_x; % Kanal 4, Beschleunigungssensor CH5= VCAP_DATA(:,5)*x_x; % Kanal 5, Beschleunigungssensor CH6= VCAP_DATA(:,6)*x_x; % Kanal 6, Beschleunigungssensor CH7= VCAP_DATA(:,7)*x_x; % Kanal 7, Beschleunigungssensor CH8= VCAP_DATA(:,8)*x_x; % Kanal 8, Beschleunigungssensor CH9= VCAP_DATA(:,9)*x_x; % Kanal 9, Beschleunigungssensor CH10= VCAP_DATA(:,10)*x_x; % Kanal 10, Beschleunigungssensor CH11= VCAP_DATA(:,11)*x_x; % Kanal 11, Beschleunigungssensor CH12= VCAP_DATA(:,12)*x_x; % Kanal 12, Beschleunigungssensor CH13= VCAP_DATA(:,13)*x_x; % Kanal 13, Beschleunigungssensor CH14= VCAP_DATA(:,14)*x_x; % Kanal 14, Beschleunigungssensor CH15= VCAP_DATA(:,15)*10; % Kanal 15, Beschleunigungssensor %CH16= VCAP_DATA(:,16)*x_x; % Kanal 16, Beschleunigungssensor %###### Ohne Analogsignale für Durchmesser und Geschwindigkeit %CH9= VCAP_DATA(:,9)*x_V; % Analogsignal für Geschwindigkeit %CH10= VCAP_DATA(:,10)*x_D; % Analogsignal für Durchmesser % Definition der Zeitachse. t_scale= 0:Ts:t; % Definition der Frequenzachse f_scale= fs/N*(0:N/2-1); % Zeitsignale: CH1_t = CH1(k*N+1:(k+1)*N); CH2_t = CH2(k*N+1:(k+1)*N); CH3_t = CH3(k*N+1:(k+1)*N); CH4_t = CH4(k*N+1:(k+1)*N); CH5_t = CH5(k*N+1: (k+1)*N); CH6_t = CH6(k*N+1: (k+1)*N); CH7_t = CH7(k*N+1: (k+1)*N); CH8_t = CH8(k*N+1: (k+1)*N); CH9_t = CH9(k*N+1: (k+1)*N); CH10_t = CH10(k*N+1: (k+1)*N); CH11_t = CH11(k*N+1: (k+1)*N); CH12_t = CH12(k*N+1: (k+1)*N); CH13_t = CH13(k*N+1: (k+1)*N); CH14_t = CH14(k*N+1: (k+1)*N); CH15_t = CH15(k*N+1: (k+1)*N); %CH16_t = CH16(k*N+1: (k+1)*N); %### Ohne Analogsignale für Geschwindigkeit und Wickeldurchmesser %CH9_t = CH9(k*N+1: (k+1)*N); %CH10_t = CH10(k*N+1: (k+1)*N); % Aktuelle Geschwindigkeit, Wickeldurchmesser und % Wickeldrehfrequenz aus dem gewählten Frame: %Vact= 1/length(CH9_t)*sum(CH9_t); %Dact= 1/length(CH10_t)*sum(CH10_t); %fact= (1000/60)*Vact/(Dact*pi); % Frequenzspektren % Berücksichtigung der Phasenlage (Montage der Aufnehmer) beachten: % Bei Drehung einzelner Sensoren um 180° wird der fft-Wert mit % (-1) oder mit exp(i*pi) multipliziert: % z.Besp.: x-Richtung: CH1_f = fft(CH1_t,N)/(N/2) ist um 180° verdreht % gegenüber den anderen Sensoren in x-Richtung montiert, dann % CH1_f = fft(CH1_t,N)/(N/2)* (-1); % ist identisch mit % CH1_f = fft(CH1_t,N)/(N/2)*exp(i*pi); % Dann zeigen alle Sensoren in die gleiche % Richtung. CH1_f = fft(CH1_t,N)/(N/2); % Kanal 1 CH2_f = fft(CH2_t,N)/(N/2); % Kanal 2 CH3_f = fft(CH3_t,N)/(N/2); % Kanal 3 CH4_f = fft(CH4_t,N)/(N/2); % Kanal 4 CH5_f = fft(CH5_t,N)/(N/2); % Kanal 5 CH6_f = fft(CH6_t,N)/(N/2)*(-1); % Kanal 6 Richtungskorrektur CH7_f = fft(CH7_t,N)/(N/2)*(-1); % Kanal 7 Richtungskorrektur CH8_f = fft(CH8_t,N)/(N/2); % Kanal 8 CH9_f = fft(CH9_t,N)/(N/2); % Kanal 9 CH10_f = fft(CH10_t,N)/(N/2)*(-1); % Kanal 10 Richtungskorrektur CH11_f = fft(CH11_t,N)/(N/2); % Kanal 11 CH12_f = fft(CH12_t,N)/(N/2); % Kanal 12 CH13_f = fft(CH13_t,N)/(N/2); % Kanal 13 CH14_f = fft(CH14_t,N)/(N/2); % Kanal 14 CH15_f = fft(CH15_t,N)/(N/2); % Kanal 15 %CH16_f = fft(CH16_t,N)/(N/2); % Kanal 16 % ##### Liegen gemischte Signal vor ?? % z.Bsp. Beschleunigungs- und Wegsignale % Die Spektren von Wegsignalen werden durch doppeltes Differenzieren % in Beschleunigungsspektren umgerechnet: % A) Differenzieren: x(t) --> x'(t) % CHx_f= CHx_f.*(i*2*pi*f_scale'); von f_scale den transponierten Vektor % B) Doppeltes Differenzieren: x(t) --> x''(t) % CHx_f= CHx_f.*(i*2*pi*f_scale').^2; % #### Beachte: % Wenn differenziert wird, dann muß die Abzissensalierung in % der nachfolgenden figure (5) der veränderten Amplitudengröße % angepaßt werden. % Doppeltes differenzieren: Umrechnung von % Wegsignale in Beschleunigungssignale: %CH1_f = CH1_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 1 %CH2_f = CH2_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 2 %CH3_f = CH3_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 3 %CH4_f = CH4_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 4 %CH5_f = CH5_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 5 mit Richtungskorrektur %CH6_f = CH6_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 6 %CH7_f = CH7_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 7 %CH8_f = CH8_f(1:N/2).*(i*2*pi*f_scale').^2; % Kanal 8 % ##### Darstellung der Zeit- und Frequenzsignale: % A) Zeitsignale aller Kanäle (CH1 - CH8) figure (5) set(5,'Color','White') %subplot(3,1,1) % Darstellung der Zeitsignale vom gewählten Frame %plot(t_scale(1:end-1),CH1_t,'c') %h= title(['Gewählter Frame: ', num2str(k),' Durchmesser: ', num2str(Dact),' [mm]'... % ' Geschwindigkeit: ', num2str(Vact),' [m/min]']); %set(h,'Fontsize',14) %hold on %plot(t_scale(1:end-1), CH2_t,'b') %plot(t_scale(1:end-1), CH3_t,'r') %plot(t_scale(1:end-1), CH4_t,'k') %plot(t_scale(1:end-1), CH5_t,'b') %plot(t_scale(1:end-1), CH6_t,'r') %plot(t_scale(1:end-1), CH7_t,'k') %plot(t_scale(1:end-1), CH8_t,'g') %legend('CH1','CH2','CH3','CH4',... % 'CH5','CH6','CH7','CH8') %hold off subplot(3,1,1) % Frequnzspektrum von CH1 bis CH4 plot(f_scale, abs(CH15_f(1:N/2)),'r') h= title(['Trigersignal S-Walze ']); set(h,'FontSize',14) axis([0 250 0 10]); % Skalierung des Frequenzspektrums hold on legend('Triggersignal') h= ylabel('[V]'); set(h,'FontSize',10) h= xlabel('Frequenz [Hz]'); set(h,'FontSize',10) %delta = 10; %set(gca,'xtick',([min(xlim):delta:max(xlim)])) grid on hold off % B) Frequenzsignale der Kanäle (CH1 - CH4) subplot(3,1,2) % Frequnzspektrum von CH1 bis CH4 plot(f_scale, abs(CH1_f(1:N/2)),'r') h= title(['Frequenzspektrum Soft Kalander Varel PM4 Führerseite (FS)']); set(h,'FontSize',14) axis([0 250 0 10000]); % Skalierung des Frequenzspektrums hold on plot(f_scale, abs(CH2_f(1:N/2)),'g') plot(f_scale, abs(CH3_f(1:N/2)),'b') plot(f_scale, abs(CH4_f(1:N/2)),'c') plot(f_scale, abs(CH5_f(1:N/2)),'m') plot(f_scale, abs(CH6_f(1:N/2)),'k') legend('S-Walze Außenlager horizontal TS', 'S-Walze Außenlager vertikal TS', 'S-Walze Achse horizontal TS', 'S-Walze Achse vertikal TS',... 'Heizwalze Außenlager horizontal TS', 'Heizwalze Außenlager vertikal TS' ) h= ylabel('[mm/s²]'); set(h,'FontSize',10) h= xlabel('Frequenz [Hz]'); set(h,'FontSize',10) %delta = 10; %set(gca,'xtick',([min(xlim):delta:max(xlim)])) grid on hold off % B) Frequenzsignale der Kanäle (CH5 - CH6) subplot(3,1,3) % Frequenzspektrum von CH5 bis CH8 plot(f_scale, abs(CH7_f(1:N/2)),'r') h= title(['Frequenzspektrum Soft Kalander Varel PM4 Antriebsseite (TS)']); set(h,'FontSize',14) axis([0 250 0 10000]); % Skalierung des Frequenzspektrums hold on plot(f_scale, abs(CH8_f(1:N/2)),'g') plot(f_scale, abs(CH9_f(1:N/2)),'b') plot(f_scale, abs(CH10_f(1:N/2)),'c') plot(f_scale, abs(CH11_f(1:N/2)),'m') plot(f_scale, abs(CH12_f(1:N/2)),'k') legend('S-Walze Außenlager horizontal FS', 'S-Walze Außenlager vertikal FS', 'S-Walze Achse horizontal FS', 'S-Walze Achse vertikal FS',... 'Heizwalze Außenlager horizontal FS', 'Heizwalze Außenlager vertikal FS' ) h=ylabel('[mm/s²]'); set(h,'FontSize',10) h= xlabel('Frequenz [Hz]'); set(h,'FontSize',10) %delta = 10; %set(gca,'xtick',([min(xlim):delta:max(xlim)])) grid on hold off %############################### % A) Die Darstellung der Orbits erolgt mit Bescheunigungssignalen, % B) Mit der Zoom-Funktion wird aus der obigen Grafik die Frequenz % ausgewählt. % Die Kommastelle der Frequenz muß stimmen !!!!! % C) Nach der Eingabeaufforderung wird die ausgewöählte Frequenz % in das Command Window eingetippt. %### An dieser Stelle: Programm STOPPT, % Programm erwartet die Eingabe einer Frequenz pause; % Prog. wird angehalten f_select= input('Eingabe einer selektierten Frequenz aus dem dargestellen Spektrum: '); L= find(f_scale == f_select); % Suche nach dem Index zu der Frequenz % Für die Visualisierzung können Maßstabsfaktoren gesetzt werden, % sie verändern nur die Orbitgröße nicht die Phasenlage: quo=1; % für die Tragwalzen quw=1; % für den Wickel qud=1; % für die Druckwalze % Amplituden von der selektierten Frequenz: A1= 1/quo *abs(CH1_f(L)); % Kanal 1 A2= 1/quo *abs(CH2_f(L)); % Kanal 2 A3= 1/quo *abs(CH3_f(L)); % Kanal 3 A4= 1/quo *abs(CH4_f(L)); % Kanal 4 A5= 1/quw *abs(CH5_f(L)); % Kanal 5 A6= 1/quw *abs(CH6_f(L)); % Kanal 6 A7= 1/qud *abs(CH7_f(L)); % Kanal 7 A8= 1/qud *abs(CH8_f(L)); % Kanal 8 A9= 1/qud *abs(CH9_f(L)); % Kanal 9 A10= 1/qud *abs(CH10_f(L)); % Kanal 10 A11= 1/qud *abs(CH11_f(L)); % Kanal 11 A12= 1/qud *abs(CH12_f(L)); % Kanal 12 Amplituden = [A1;A2;A3;A4;A5;A6;A7;A8;A9;A10;A11;A12]; % Amplituden von der selektierten Frequenz % Phasenlage von der selektierten Frequenz: phi1= angle(CH1_f(L)); % Kanal 1 phi2= angle(CH2_f(L)); % Kanal 2 phi3= angle(CH3_f(L)); % Kanal 3 phi4= angle(CH4_f(L)); % Kanal 4 phi5= angle(CH5_f(L)); % Kanal 5 phi6= angle(CH6_f(L)); % Kanal 6 phi7= angle(CH7_f(L)); % Kanal 7 phi8= angle(CH8_f(L)); % Kanal 8 phi9= angle(CH9_f(L)); % Kanal 9 phi10= angle(CH10_f(L)); % Kanal 10 phi11= angle(CH11_f(L)); % Kanal 11 phi12= angle(CH12_f(L)); % Kanal 12 Phase = [phi1;phi2;phi3;phi4;phi5;phi6;phi7;phi8;phi9;phi10;phi11;phi12]; % Phasenlagen in Bogenmaß %#### Darstellung der Maschiengeometrie % Die obigen Eingaben von TW1, TW2, DW und Dact werden verechnet conf= pi/180; % Umrechung von Winkel und Bogenmaß D= 0.5*(TW1 + TW2) + S; DD= D^2; C= sqrt(D^2-TWH^2); X= C/D; % Walzenabstand, Höhenwickel; ALPHA= acos(X)/conf; F=0.5*(TW1 + Dact); FF=F^2; E=0.5*(TW2 + Dact); EE=E^2; WW2=(DD + EE - FF)/(2*D*E); WW3= (FF + EE - DD)/(2*F*E); WW2= acos(WW2)/conf; WW3= acos(WW3)/conf; % Winkel im Wickelzentrum zu den NIP's Wwickel=WW3; % Berechnung der Winkel; ZETAA = WW3 - 90 + WW2 - ALPHA; PHII = 90 + ALPHA - WW2; TAU = 180 - PHII - ZETAA; % Zurückgelegter Winkel auf den Tragwalzen % PHI: NIP-Winkel auf TW2 % ZETA: NIP-Winkel auf TW1 ZETA = 90 - ZETAA; PHI = 90 - PHII; % WINKEL=[ZETA(x);PHI(x);Wwickel(x)] %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Orbitdarstellung der Schwingbewegung tt= 0:0.0001:1/f_select; % Orbit-Berechnung von der selektierten Frequenz CH1_H = A1*cos(2*pi*f_select*tt + phi1); % Kanal 1 horizontal CH2_V = A2*sin(2*pi*f_select*tt + phi2); % Kanal 2 vertikal CH3_H = A3*cos(2*pi*f_select*tt + phi3); % Kanal 3 horizontal CH4_V = A4*sin(2*pi*f_select*tt + phi4); % Kanal 4 vertikal CH5_H = A5*cos(2*pi*f_select*tt + phi5); % Kanal 5 horizontal CH6_V = A6*sin(2*pi*f_select*tt + phi6); % Kanal 6 vertikal CH7_H = A7*cos(2*pi*f_select*tt + phi7); % Kanal 7 horizontal CH8_V = A8*sin(2*pi*f_select*tt + phi8); % Kanal 8 vertikal CH9_H = A9*cos(2*pi*f_select*tt + phi9); % Kanal 5 horizontal CH10_V = A10*sin(2*pi*f_select*tt + phi10); % Kanal 6 vertikal CH11_H = A11*cos(2*pi*f_select*tt + phi11); % Kanal 7 horizontal CH12_V = A12*sin(2*pi*f_select*tt + phi12); % Kanal 8 vertikal % Geometrie der Tragwalzenkonturen AW = Dact/2* cos(2*pi*f_select*tt) + 0.5*(TW2+Dact)*cos(PHI*conf); % Wickel: x-Koordinate BW = Dact/2* sin(2*pi*f_select*tt) + 0.5*(TW2+Dact)*sin(PHI*conf); % Wickel: y-Koordinate AT1 = TW1/2*cos(2*pi*f_select*tt) + C; % TW1: x-Koordinate BT1 = TW1/2*sin(2*pi*f_select*tt) - TWH; % TW1: y-Koordinate AT2 = TW2/2* cos(2*pi*f_select*tt); % TW2: x-Koordinate BT2 = TW2/2* sin(2*pi*f_select*tt); % TW2: y-Koordinate DW1 = 0.5*(TW2 + Dact)*cos(PHI*conf) + DW/2*cos(2*pi*f_select*tt); % DW-x-Koordinate DW2 = 0.5*(Dact + DW) + 0.5*(TW2 + Dact)*sin(PHI*conf) + DW/2*sin(2*pi*f_select*tt); % DW-y-Koordinate %### Darstellung der Orbitformen figure (7) set(7,'Color','White','name','Orbitdarstellung') ax= -TW2/2; ay= -TW1; % Koodinatenverschiebung h= plot(AT2-1000,BT2-500,'r'); set(h,'LineWidth',2) hold on achsen= 2000; % Standardskalierung %achsen= 4000; axis([-achsen achsen -achsen achsen]) axis square % Damit bleibt die Grafik unverzerrt axis off h= plot(AT1+300,BT1+10,'r'); set(h,'LineWidth',2) h= plot(AW-1320,BW-420,'k'); % Darstellung des Wickels set(h,'LineWidth',2) h= plot(DW1+620,DW2-720,'k'); % Darstellung der Druckwalze set(h,'LineWidth',2) % Kommentar zu der Darstellung: h=title(['Softkalander PM4 Varel, Orbitdarstellung, ', 'Frame: ', num2str(k)]); set(h,'FontSize',20) %xlabel(['Durchmesser : ',num2str(Dact), ' Gesch. : ',num2str(Vact),' select Freq. : ', num2str(f_select)]) h=text(-2000, -1850,['Selektierte Frequenz: ', num2str(f_select),' [Hz] ']); set(h,'FontSize',16) h=text(-2000, -2000,['Durchmesser Heizwalze: ', num2str(Dact),' [mm]']); set(h,'Fontsize',16) h=text(-2000, -2150,['Durchmesser S-Walze: ', num2str(TW1),' [mm]']); set(h,'Fontsize',16) % Orbitdarstellungen a= plot(CH1_H - 1000, CH2_V - 500,'b'); % S-Walze FS set(a,'LineWidth',2); a1= plot(CH3_H - 1000, CH4_V - 500,'g'); % S-Walze FS-Achse set(a,'LineWidth',2); a2= plot(CH7_H + C+ 300, CH8_V - TWH + 10,'b'); % S-Walze TS set(a,'LineWidth',2); a3= plot(CH9_H + C+ 300, CH10_V - TWH + 10,'g'); % S-Walze TS-Achse set(a,'LineWidth',2); a4= plot(CH5_H -1320 +0.5*(TW2 + Dact)*cos(PHI*conf), CH6_V - 420 + 0.5*(TW2 + Dact)*sin(PHI*conf),'m'); % Heizwalze FS set(a,'LineWidth',2); a5= plot(CH11_H + 620 + 0.5*(TW2 + Dact)*cos(PHI*conf), CH12_V + 0.5*(Dact + DW) - 720 + 0.5*(TW2 + Dact)*sin(PHI*conf),'m'); % Heizwalze TS set(a,'LineWidth',2); %legend ('S-Walze Außenlager', 'S-Walze Achse', 'Heizwalze Außenlager' ) legend([a a1 a4],{'S-Walze Außenlager', 'S-Walze Achse', 'Heizwalze Außenlager'}) %nur selektierte Plots in der Legende darstellen grid on; h= text(-1100, -1000,'Führerseite (FS)'); set(h,'Fontsize',16) h= text(800, -500,'Antriebsseite (TS)'); set(h,'Fontsize',16) % Markierungen, Umlaufsinn der Orbits: A= 1; % Grün: Markierung für den Start Z= 20; % Rot: Markierung für die Umlaufrichtung % Tragwalze 2 (kleine Walze): [CH1_H, CH2_V] % Tragwalze 1 (große Walze): [CH3_H, CH4_V] % Wickkung: [CH5_H, CH6_V] % Druckwalze: [CH7_H, CH8_V] plot(CH1_H(A) - 1000, CH2_V(A) - 500,'g.','Markersize',20) %S-Walze FS plot(CH1_H(Z) - 1000, CH2_V(Z) - 500,'r.','Markersize',20) %S-Walze FS plot(CH3_H(A) - 1000, CH4_V(A) - 500,'g.','Markersize',20) %S-Walze FS Achse plot(CH3_H(Z) - 1000, CH4_V(Z) - 500,'r.','Markersize',20) %S-Walze FS Achse plot(CH7_H(A) + C+ 300, CH8_V(A) - TWH + 10,'g.','Markersize',20) % S-Walze TS plot(CH7_H(Z) + C+ 300, CH8_V(Z) - TWH + 10,'r.','Markersize',20) % S-Walze TS plot(CH9_H(A) + C+ 300, CH10_V(A) - TWH + 10,'g.','Markersize',20) % S-Walze TS plot(CH9_H(Z) + C+ 300, CH10_V(Z) - TWH + 10,'r.','Markersize',20) % S-Walze TS plot(CH5_H(A) - 1320 + 0.5*(TW2 + Dact)*cos(PHI*conf), CH6_V(A) - 420 + 0.5*(TW2 + Dact)*sin(PHI*conf),'g.','Markersize',20) plot(CH5_H(Z) - 1320 + 0.5*(TW2 + Dact)*cos(PHI*conf), CH6_V(Z) - 420 + 0.5*(TW2 + Dact)*sin(PHI*conf),'r.','Markersize',20) plot(CH11_H(A) + 620 + 0.5*(TW2 + Dact)*cos(PHI*conf), CH12_V(A) + 0.5*(DW + Dact) - 720 + 0.5*(TW2 + Dact)*sin(PHI*conf),'g.','Markersize',20) plot(CH11_H(Z) + 620 + 0.5*(TW2 + Dact)*cos(PHI*conf), CH12_V(Z) + 0.5*(DW + Dact) - 720 + 0.5*(TW2 + Dact)*sin(PHI*conf),'r.','Markersize',20) hold off