%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%MÉTODO 2%%%%%%%%%%%%%%%%% %%%%%%%%%%%%Cambios en los modos%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %FUNCIÓN CAMBIOS EN LOS MODOS for i=1:1:m divisor=1; %El término "divisor", permite,si se desea, escalar los modos cc1=(1/divisor)*matsin2(:,i); cc2=(1/divisor)*matcon2(:,i); figure(i) hold on for k=1:1:nbarra plot(coordenadas(inic_fin(k,1):inic_fin(k,2),1),zeros(length(coordenadas(inic_fin(k,1):inic_fin(k,2),1))),'k') %hold on plot(coordenadas(inic_fin(k,1):inic_fin(k,2),1),cc1(inic_fin(k,1):inic_fin(k,2)),'b-o') plot(coordenadas(inic_fin(k,1):inic_fin(k,2),1),cc2(inic_fin(k,1):inic_fin(k,2)),'r-o') %legend('Estructura','Modo sin daño','Modo con daño') %axis([0 200 0 20 -0.5 2]) end end disp('La diferencia de los modos en los dos estados da idea de la presencia de daño') disp('La mayor diferencia entre ambas representaciones puede dar idea de la localizacion del daño') %opcion=menu('¿Desea cambial el sentido de algún modo?','Sí',... %'No'); %Cargar archivo con los vectores de frecuencia y amortiguación %if opcion==1 % for i=1:m % k(i)=i; % end %cambiar=menu('¿Qué modo desea cambiar?','Modo',int2str(k)) %else %Salida del programa % break %end