Drude模型仿真代码

 imx=[ 4.09643975823885e+002 
 5.11263069660740e+002   
 6.08673401861737e+002   
 7.06083734062734e+002    
 8.03494066263731e+002    
 9.00877928265761e+002    
 9.94053028631932e+002   
 1.10837781797326e+003   
 1.18461199100013e+003    
 1.27352538933251e+003    
 1.35820355583006e+003   
 1.44288172232761e+003    
 1.52755988882516e+003    
 1.61223805532272e+003   
 1.69691622182027e+003   
 1.78159438831782e+003   
 1.86203732298054e+003   
 1.95095072131292e+003];

 imy=[ 1.45586094322141e-001
 1.46615490948665e+000
  1.53380097351510e+000
   1.60144703754355e+000
 1.66909310157200e+000
2.98672078351785e+000
    3.05142571432776e+000
    4.38081792914812e+000
 4.43375832708341e+000
   5.74550374259201e+000
   7.05430802488200e+000
  8.36311230717195e+000
  9.67191658946194e+000
  1.09807208717519e+001
  1.22895251540418e+001
 1.35983294363318e+001
   1.49041925854032e+001
 1.62159380009118e+001];
 plot（imximy‘ko‘‘MarkerFaceColor‘‘k‘‘MarkerSize‘5）; 
% plot（imximy‘r‘）;  
 axis xy;
 axis（[400 2000 -250 20]）;
 
 hold on
 
 rex=[ 
 1.52733195754717e+003 
 1.61072228773585e+003   
 1.69457547169811e+003   
 1.77936910377359e+003   
 1.86275353773585e+003   
 1.94942216981132e+003 
 1.443481721698