% Comparacao entre Teoria dos Pontos Fixos 
% - Teoria de optimizacao nao linear

clc
clear

global Plothandle m1 k1 c1 ksi1 m2 ome1 ome2 C11 omega iomegmax

m1=10;
k1=500000/4;
ksi1=0.0;
c1=ksi1*2*m1*sqrt(k1/m1);
mu=0.1;
m2=mu*m1;
ome1=2*pi*10;
ome2=2*pi*100;
C11=(ome1+ome2)/2;
omega=0:0.5:2*ome2;

for iomeg=1:length(omega)
    beta(iomeg)=omega(iomeg)/sqrt(k1/m1);
    H1(iomeg)=1./(-m1*(omega(iomeg))^2+k1+2*m1*ksi1*sqrt(k1/m1)*omega(iomeg)*i);
end

subplot(2,1,1)
plot(beta,abs(H1),'k-')
axis([0.5 1.5 0 5e-5])
xlabel('omega/omegan')
ylabel('abs(H1)')
title('Comparacao entre TPF e TPNL com ksi1=0')

pause

hold on

Plothandle = plot(beta,abs(H1),'r-','EraseMode','xor')

%TPNL

lam = [150 100];
lambda=fminsearch('fobj',lam);


% %Teoria dos Pontos Fixos
% 
% alphao=1/(1+mu);
% k2=m2*(alphao*sqrt(k1/m1))^2;
% Xi2o=1/alphao*sqrt(3*mu/(8*(1+mu)^3));
% c2=Xi2o*2*m2*sqrt(k2/m2);
% 
%  for iomeg=1:length(omega)
%      H11(iomeg)=(-omega(iomeg)^2*m2+k2+i*omega(iomeg)*c2)/((-omega(iomeg)^2*m2+k2)*(-omega(iomeg)^2*m1+k1)-omega(iomeg)^2*m2*k2+i*omega(iomeg)*c2*(-omega(iomeg)^2*m1+k1-omega(iomeg)^2*m2));
%  end
% 
% plot(beta,abs(H11),'y-')
% 
% 
% %c2=0
% 
% c2=0;
% 
%  for iomeg=1:length(omega)
%      H11(iomeg)=(-omega(iomeg)^2*m2+k2+i*omega(iomeg)*c2)/((-omega(iomeg)^2*m2+k2)*(-omega(iomeg)^2*m1+k1)-omega(iomeg)^2*m2*k2+i*omega(iomeg)*c2*(-omega(iomeg)^2*m1+k1-omega(iomeg)^2*m2));
%  end
%  
% plot(beta,abs(H11),'b-')
% 
% 
% %c2=infinito
% 
% 
% c2=1000000000000;
% 
%  for iomeg=1:length(omega)
%      H11(iomeg)=(-omega(iomeg)^2*m2+k2+i*omega(iomeg)*c2)/((-omega(iomeg)^2*m2+k2)*(-omega(iomeg)^2*m1+k1)-omega(iomeg)^2*m2*k2+i*omega(iomeg)*c2*(-omega(iomeg)^2*m1+k1-omega(iomeg)^2*m2));
%  end
%  
% plot(beta,abs(H11),'g-')
% legend('sem neutralizador','TPF','TONL','C2=0','C2=infinito')
% 
% 
% 
% 
% 
% %-----------------------------------ksi1 nao igual zero---
% 
% ksi1=0.05;
% c1=ksi1*2*m1*sqrt(k1/m1);
% 
% for iomeg=1:length(omega)
%     beta(iomeg)=omega(iomeg)/sqrt(k1/m1);
%     H1(iomeg)=1./(-m1*(omega(iomeg))^2+k1+2*m1*ksi1*sqrt(k1/m1)*omega(iomeg)*i);
% end
% 
% 
% subplot(2,1,2)
% plot(beta,abs(H1),'k-')
% axis([0.5 1.5 0 5e-5])
% xlabel('omega/omegan')
% ylabel('abs(H1)')
% title('Comparacao entre TPF e TPNL com ksi1 nao igual zero')
% 
% pause
% 
% hold on
% 
% Plothandle = plot(beta,abs(H1),'r-','EraseMode','xor')
% 
% % 
% %TPNL
% 
% lam = [150 100];
% % 
% lambda2=fminsearch('fobj2',lam);
% 
% 
% 
% 
% 
% %Teoria dos Pontos Fixos
% 
% alphao=1/(1+mu);
% k2=m2*(alphao*sqrt(k1/m1))^2;
% Xi2o=1/alphao*sqrt(3*mu/(8*(1+mu)^3));
% c2=Xi2o*2*m2*sqrt(k2/m2);
% 
%  for iomeg=1:length(omega)
%        H11(iomeg)=(-omega(iomeg)^2*m2+k2+i*omega(iomeg)*c2)/((-omega(iomeg)^2*m2+k2)*(-omega(iomeg)^2*m1+k1+i*omega(iomeg)*c1)- ...
%                   omega(iomeg)^2*m2*k2+i*omega(iomeg)*c2*(-omega(iomeg)^2*m1+k1+i*omega(iomeg)*c1-omega(iomeg)^2*m2));
%  end
% 
% plot(beta,abs(H11),'y-')
% 
% legend('sem neutralizador','TPF','TONL')
% 
% 
% %--------------------------------------------------------------------------
% 
% % 
% % omegan2=lambda(1);
% % ksi2=lambda(2)/C11;
% % c20=2*ksi2*m2*omegan2
% % k20=omegan2^2*m2
% % om=omega(iomegmax)
% % % c2 =76.2917
% % % k2 =3.9386e+004
% % % om=236;
% % 
% % C2 =c20-50:c20+50;
% % K2 =k20-5e4:0.1e4:k20+5e4;                
% % [c2,k2] = meshgrid(C2,K2); 
% % 
% % 
% % for iomeg=1:length(omega)
% %     H11(iomeg)=(-omega(iomeg)^2*m2+k2+i*omega(iomeg)*c2)/((-omega(iomeg)^2*m2+k2)*(-omega(iomeg)^2*m1+k1)- ...
% %                   omega(iomeg)^2*m2*k2+i*omega(iomeg)*c2*(-omega(iomeg)^2*m1+k1-omega(iomeg)^2*m2));
% % end
% % H11=max(abs(H11));
% % 
% % 
% % 
% % subplot(2,2,3)
% % plot3(c2,k2,H11)
% % hold on
% % H110=(-om^2*m2+k20+i*om*c20)./((-om^2*m2+k20).*(-om^2*m1+k1)-om^2*m2.*k20+i*om.*c20.*(-om^2*m1+k1-om^2*m2));
% % plot3(c20,k20,abs(H110),'go')
% % 
% % xlabel('c2')
% % ylabel('k2')
% % zlabel('abs(H11)')
% % rotate3d
% % 




  





 















 
