wouter
/
Matlab-Prologix-scripts
1
0
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
42 Commits
2 Branches
60MB
Tree: 32323781b2
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '32323781b2'
${ noResults }
Matlab-Prologix-scripts/MOD2-files/EXP5/Transfer_Functions_Matlab/Transfer_Function_RLC.m

334 lines
9.5KB
Raw Blame History 

  1. clear;

  2. clc;

  3. 

  4. fprintf('Bode diagram of the RLC-circuit\n\n');

  5. fprintf('Figure 1: Bode diagram theory+asymptote + average measurement\n')

  6. fprintf('Figure 2: Bode diagram average measurement+SD + theory\n');

  7. fprintf('Figure 3: Bode diagram average measurement+SD\n\n');

  8. fprintf('N.B. SD=0 is not displayed in the figures!\n\n');

  9. 

  10. fprintf('\nWatch out! This script searches the file in "D:\"\n\n');

  11. 

  12. f_name=input('Filename: ','s');

  13. f_name1=sprintf('%s%s',f_name,'_theorie');

  14. f_name2=sprintf('%s%s',f_name,'_meting');

  15. 

  16. wd=cd;

  17. cd 'D:\'

  18. 

  19. [th]=dlmread(f_name1,' ',0,0);

  20. [met]=dlmread(f_name2,' ',0,0);

  21. 

  22. cd(wd);

  23. 

  24. f_th=th(1,:);

  25. amp_th=th(2,:);

  26. fase_th=th(3,:);

  27. asym=th(4,:);

  28. if mean(met(1,:))==0;

  29.     nr_met=0;

  30. end;

  31. 

  32. if mean(met(1,:))~=0;

  33.     size_met=size(met);

  34.     nr_met=(size_met(1)-1)/2;

  35.     nr_punt=size_met(2);

  36.     f_m=met(1,:);

  37.     

  38.     for i=1:nr_met;

  39.         for j=1:nr_punt;

  40.             amp_met(i,j)=met(i+1,j);

  41.             fase_met(i,j)=met(i+1+nr_met,j);

  42.         end;

  43.     end;

  44.     if nr_met~=1;

  45.         for i=1:nr_punt;

  46.             av_amp(i)=mean(amp_met(:,i));

  47.             av_log_amp(i)=20*log10(av_amp(i));

  48.             av_fase(i)=mean(fase_met(:,i));

  49.         end;

  50.         sd_amp=std(amp_met,0,1);

  51.         for i=1:nr_punt;

  52.             if sd_amp(i)==0;

  53.                 sd_log_amp(i)=0;

  54.             elseif av_amp(i)==1;

  55.                 av_amp(i)=1.0001;

  56.             else

  57.                 sd_log_amp(i)=2/log(av_amp(i))*sd_amp(i);

  58.             end;

  59.         end;

  60.         sd_fase=std(fase_met,0,1);

  61.     else 

  62.         av_log_amp=20*log10(amp_met);

  63.         av_fase=fase_met;

  64.         sd_log_amp=zeros(1,nr_punt);

  65.         sd_fase=zeros(1,nr_punt);

  66.     end;

  67. end;

  68. 

  69. figure(1);

  70. set(clf,'PaperType','A4');

  71. set(clf,'Color',[1,1,1]);

  72. %set(clf,'PaperPosition',[1.13386 3.97146 5.5 4.125]);

  73. subplot(2,1,1);

  74. if mean(met(1,:))~=0;

  75.     graf=semilogx(f_th,amp_th,'r',f_th,asym,'g',f_m,av_log_amp,'.-b','LineWidth',1);

  76.     legend([graf],'theory','asymp.','measurement',3);

  77. else 

  78.     graf=semilogx(f_th,amp_th,'r',f_th,asym,'g','LineWidth',1);

  79.     legend([graf],'theory','asymp.',3);

  80. end;

  81. ylabel('Amplitude [dB]','Fontsize',10)

  82. 

  83. % grid bepaling

  84. startf=floor(log10(f_th(1)));

  85. sizef=size(f_th);

  86. endf=ceil(log10(f_th(sizef(2))));

  87. xf(1)=10^startf;

  88. i=startf+1;

  89. count=2;

  90. while i<=endf;

  91.     xf(count)=10^(i);

  92.     count=count+1;

  93.     i=i+1;

  94. end;

  95. set(gca, 'XTick',xf,'Fontsize',9);

  96. 

  97. if mean(met(1,:))~=0;

  98.     min_amp_th=floor(min(amp_th));

  99.     min_amp_met=floor(min(av_log_amp-abs(sd_log_amp)));

  100.     min_amp=min(min_amp_th,min_amp_met);

  101.     max_amp_th=ceil(max(amp_th));

  102.     max_amp_met=ceil(max(av_log_amp+abs(sd_log_amp)));

  103.     max_amp=max(max_amp_th,max_amp_met)+2;

  104. else

  105.     min_amp=floor(min(amp_th));

  106.     max_amp=ceil(max(amp_th))+2;

  107. end;

  108.     

  109. axis([10^startf,10^endf,min_amp,max_amp]);

  110. grid on;

  111. title('Bode diagram RLC-circuit','Fontsize',10,'Fontweight','bold');

  112. 

  113. subplot(2,1,2);

  114. if mean(met(1,:))~=0;

  115.     graf=semilogx(f_th,fase_th,'r',f_m,av_fase,'.-b','LineWidth',1);

  116.     legend([graf],'theory','measurement',3);

  117. else

  118.     graf=semilogx(f_th,fase_th,'LineWidth',1);

  119.     legend([graf],'theory',3);

  120. end;

  121. xlabel('Frequency [Hz]','Fontsize',10);

  122. ylabel('Phase [rad]','Fontsize',10);

  123. set(gca, 'XTick',xf,'Fontsize',9);

  124. 

  125. if mean(met(1,:))~=0;

  126.     min_fase_th=floor(min(fase_th));

  127.     min_fase_met=floor(min(av_fase-abs(sd_fase)));

  128.     min_fase=min(min_fase_th,min_fase_met);

  129.     max_fase_th=ceil(max(fase_th));

  130.     max_fase_met=ceil(max(av_fase+abs(sd_fase)));

  131.     max_fase=max(max_fase_th,max_fase_met);

  132. else

  133.     min_fase=floor(min(fase_th));

  134.     max_fase=ceil(max(fase_th));

  135. end;  

  136.     

  137. axis([10^startf,10^endf,min_fase,max_fase]);

  138. grid on;

  139. 

  140. if nr_met~=0;

  141.     figure(2);

  142.     set(clf,'PaperType','A4');

  143.     set(clf,'Color',[1,1,1]);

  144. %    set(clf,'PaperPosition',[1.13386 3.97146 5.5 4.125]);

  145.     subplot(2,1,1);

  146.     graf=semilogx(f_m,av_log_amp,'.-b',f_th,amp_th,'-r','LineWidth',1);

  147.     legend([graf],'measurement','theory',3);

  148.     hold on;

  149.     width_sd=(log10(f_m(2))-log10(f_m(1)))/3;

  150.     for i=1:nr_punt;

  151.         if sd_log_amp(i)~=0;

  152.             x=[f_m(i),f_m(i)];

  153.             y=[av_log_amp(i)+sd_log_amp(i),av_log_amp(i)-sd_log_amp(i)];

  154.             xl1=[10^(log(f_m(i))/log(10)-width_sd),10^(log(f_m(i))/log(10)+width_sd)];

  155.             yl1=[av_log_amp(i)+sd_log_amp(i),av_log_amp(i)+sd_log_amp(i)];

  156.             yl2=[av_log_amp(i)-sd_log_amp(i),av_log_amp(i)-sd_log_amp(i)];

  157.             plot(x,y,'-b');

  158.             plot(xl1,yl1,'-b');

  159.             plot(xl1,yl2,'-b');        

  160.         end;

  161.     end;

  162.     hold off;

  163.     ylabel('Amplitude [dB]','Fontsize',10);

  164. 

  165.     % grid bepaling

  166.     startfm=floor(log10(f_m(1)));

  167.     sizefm=size(f_m);

  168.     endfm=ceil(log10(f_m(sizefm(2))));

  169. %    xfm(1)=10^startfm;

  170.     xfm(1)=f_m(1);

  171.     i=startfm+1;

  172.     countm=2;

  173.     flag=0;

  174.     while flag~=1;

  175.         if xfm(countm-1)<f_m(nr_punt);

  176.             xfm(countm)=10^(i);

  177.             countm=countm+1;

  178.             i=i+1;

  179.         else flag=1;

  180.         end;

  181.     end;

  182.     xfm(countm-1)=f_m(nr_punt);

  183.     %    set(gca,'XTick',xfm,'Fontsize',9);

  184.     

  185.     fmin_th=f_th(1);

  186.     count_min=1;

  187.     while fmin_th<f_m(1);

  188.         count_min=count_min+1;

  189.         fmin_th=f_th(count_min);

  190.     end;

  191. 

  192.     size_fth=size(f_th);

  193.     size_fm=size(f_m);

  194.     fmax_th=f_th(size_fth(2));

  195.     count_max=size_fth(2);

  196.     while fmax_th>f_m(size_fm(2));

  197.         count_max=count_max-1;

  198.         fmax_th=f_th(count_max);

  199.     end;

  200. 

  201.     amp_th2=amp_th(count_min:count_max);

  202.     min_amp_met=floor(min(av_log_amp-abs(sd_log_amp)));

  203.     min_amp_th2=floor(min(amp_th2));

  204.     max_amp_met=ceil(max(av_log_amp+abs(sd_log_amp)));

  205.     max_amp_th2=ceil(max(amp_th2));

  206.     min_amp_met2=floor(min(min_amp_met,min_amp_th2))-2;

  207.     max_amp_met2=ceil(max(max_amp_met,max_amp_th2))+2;

  208. 

  209.     set(gca,'Fontsize',9);

  210.     axis([xfm(1),xfm(countm-1),min_amp_met2,max_amp_met2]);

  211.     grid on;

  212.     title('Bode diagram RLC-circuit','Fontsize',10,'Fontweight','bold');

  213. 

  214.     subplot(2,1,2);

  215.     

  216.     graf=semilogx(f_m,av_fase,'.-b',f_th,fase_th,'-r','LineWidth',1);

  217.     legend([graf],'measurement','theory',3);

  218.     hold on;

  219.     width_sd=(log10(f_m(2))-log10(f_m(1)))/3;

  220.     for i=1:nr_punt;

  221.         if sd_fase(i)~=0;

  222.             x=[f_m(i),f_m(i)];

  223.             y=[av_fase(i)+sd_fase(i),av_fase(i)-sd_fase(i)];

  224.             xl1=[10^(log10(f_m(i))-width_sd),10^(log10(f_m(i))+width_sd)];

  225.             yl1=[av_fase(i)+sd_fase(i),av_fase(i)+sd_fase(i)];

  226.             yl2=[av_fase(i)-sd_fase(i),av_fase(i)-sd_fase(i)];

  227.             plot(x,y,'-b');

  228.             plot(xl1,yl1,'-b');

  229.             plot(xl1,yl2,'-b');        

  230.         end;

  231.     end;

  232.     hold off;

  233.     

  234.     xlabel('Frequency [Hz]','Fontsize',10);

  235.     ylabel('Phase [rad]','Fontsize',10);

  236. %    set(gca, 'XTick',xfm,'Fontsize',9);

  237.  

  238.     min_fase_met=floor(min(av_fase-abs(sd_fase)));

  239.     max_fase_met=ceil(max(av_fase+abs(sd_fase)));

  240.     fase_th2=fase_th(count_min:count_max);

  241.     min_fase_th2=floor(min(fase_th2));

  242.     max_fase_th2=ceil(max(fase_th2));

  243.     min_fase_met2=floor(min(min_fase_met,min_fase_th2))-1;

  244.     max_fase_met2=ceil(max(max_fase_met,max_fase_th2));

  245.     

  246.     set(gca,'Fontsize',9);

  247.     axis([f_m(1),f_m(nr_punt),min_fase_met2,max_fase_met2]);

  248.     grid on;

  249. end;

  250. 

  251. if nr_met~=0;

  252.     figure(3);

  253.     set(clf,'PaperType','A4');

  254.     set(clf,'Color',[1,1,1]);

  255. %    set(clf,'PaperPosition',[1.13386 3.97146 5.5 4.125]);

  256.     subplot(2,1,1);

  257.     graf=semilogx(f_m,av_log_amp,'.-b','LineWidth',1);

  258. %    legend([graf],'meting',3);

  259.     hold on;

  260.     width_sd=(log10(f_m(2))-log10(f_m(1)))/3;

  261.     for i=1:nr_punt;

  262.         if sd_log_amp(i)~=0;

  263.             x=[f_m(i),f_m(i)];

  264.             y=[av_log_amp(i)+sd_log_amp(i),av_log_amp(i)-sd_log_amp(i)];

  265.             xl1=[10^(log(f_m(i))/log(10)-width_sd),10^(log(f_m(i))/log(10)+width_sd)];

  266.             yl1=[av_log_amp(i)+sd_log_amp(i),av_log_amp(i)+sd_log_amp(i)];

  267.             yl2=[av_log_amp(i)-sd_log_amp(i),av_log_amp(i)-sd_log_amp(i)];

  268.             plot(x,y,'-b');

  269.             plot(xl1,yl1,'-b');

  270.             plot(xl1,yl2,'-b');        

  271.         end;

  272.     end;

  273.     hold off;

  274.     ylabel('Amplitude [dB]','Fontsize',10);

  275. 

  276.     % grid bepaling

  277.     startfm=floor(log10(f_m(1)));

  278.     sizefm=size(f_m);

  279.     endfm=ceil(log10(f_m(sizefm(2))));

  280. %    xfm(1)=10^startfm;

  281.     xfm(1)=f_m(1);

  282.     i=startfm+1;

  283.     countm=2;

  284.     flag=0;

  285.     while flag~=1;

  286.         if xfm(countm-1)<f_m(nr_punt);

  287.             xfm(countm)=10^(i);

  288.             countm=countm+1;

  289.             i=i+1;

  290.         else flag=1;

  291.         end;

  292.     end;

  293.     xfm(countm-1)=f_m(nr_punt);

  294. %    set(gca,'XTick',xfm,'Fontsize',9);

  295.     

  296.     min_amp_met=floor(min(av_log_amp-abs(sd_log_amp)))-2;

  297.     max_amp_met=ceil(max(av_log_amp+abs(sd_log_amp)))+2;

  298. 

  299.     set(gca,'Fontsize',9);

  300.     axis([xfm(1),xfm(countm-1),min_amp_met,max_amp_met]);

  301.     grid on;

  302.     title('Bode diagram RLC-circuit (measurement)','Fontsize',10,'Fontweight','bold');

  303. 

  304.     subplot(2,1,2);

  305.     

  306.     graf=semilogx(f_m,av_fase,'.-b','LineWidth',1);

  307. %    legend([graf],'meting',3);

  308.     hold on;

  309.     width_sd=(log10(f_m(2))-log10(f_m(1)))/3;

  310.     for i=1:nr_punt;

  311.         if sd_fase(i)~=0;

  312.             x=[f_m(i),f_m(i)];

  313.             y=[av_fase(i)+sd_fase(i),av_fase(i)-sd_fase(i)];

  314.             xl1=[10^(log10(f_m(i))-width_sd),10^(log10(f_m(i))+width_sd)];

  315.             yl1=[av_fase(i)+sd_fase(i),av_fase(i)+sd_fase(i)];

  316.             yl2=[av_fase(i)-sd_fase(i),av_fase(i)-sd_fase(i)];

  317.             plot(x,y,'-b');

  318.             plot(xl1,yl1,'-b');

  319.             plot(xl1,yl2,'-b');        

  320.         end;

  321.     end;

  322.     hold off;

  323.     

  324.     xlabel('Frequency [Hz]','Fontsize',10);

  325.     ylabel('Phase [rad]','Fontsize',10);

  326. %    set(gca, 'XTick',xfm,'Fontsize',9);

  327.  

  328.     min_fase_met=floor(min(av_fase-abs(sd_fase)));

  329.     max_fase_met=ceil(max(av_fase+abs(sd_fase)))+2;

  330.  

  331.     set(gca,'Fontsize',9);

  332.     axis([f_m(1),f_m(nr_punt),min_fase,max_fase]);

  333.     grid on;

  334. end;