wouter
/
Matlab-Prologix-scripts
1
0
Atdalīts 0
Kods Problēmas 0 Izmaiņu pieprasījumi 0 Laidieni 0 Vikivietne Aktivitāte
Nevar pievienot vairāk kā 25 tēmas Tēmai ir jāsākas ar burtu vai ciparu, tā var saturēt domu zīmes ('-') un var būt līdz 35 simboliem gara.
66 Revīzijas
2 Atzari
60MB
Koks: 94f9b26577
Atzari Tagi
${ item.name }
Izveidot atzaru ${ searchTerm }
no '94f9b26577'
${ noResults }
Matlab-Prologix-scripts/MOD2-files/EXP5/Transfer_Functions_Matlab/Transfer_Function_RC.m

336 rindas
9.5KB
Neapstrādāts Vainot Vēsture 

  1. % clear;

  2. % clc;

  3. 

  4. fprintf('Bode diagram of the RC-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. 

  25. 

  26. f_th=th(1,:);

  27. amp_th=th(2,:);

  28. fase_th=th(3,:);

  29. asym=th(4,:);

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

  31.     nr_met=0;

  32. end;

  33. 

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

  35.     size_met=size(met);

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

  37.     nr_punt=size_met(2);

  38.     f_m=met(1,:);

  39.     

  40.     for i=1:nr_met;

  41.         for j=1:nr_punt;

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

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

  44.         end;

  45.     end;

  46.     if nr_met~=1;

  47.         for i=1:nr_punt;

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

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

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

  51.         end;

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

  53.         for i=1:nr_punt;

  54.             if sd_amp(i)==0;

  55.                 sd_log_amp(i)=0;

  56.             elseif av_amp(i)==1;

  57.                 av_amp(i)=1.0001;

  58.             else

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

  60.             end;

  61.         end;

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

  63.     else 

  64.         av_log_amp=20*log10(amp_met);

  65.         av_fase=fase_met;

  66.         sd_log_amp=zeros(1,nr_punt);

  67.         sd_fase=zeros(1,nr_punt);

  68.     end;

  69. end;

  70. 

  71. figure(1);

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

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

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

  75. subplot(2,1,1);

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

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

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

  79. else 

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

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

  82. end;

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

  84. 

  85. % grid bepaling

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

  87. sizef=size(f_th);

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

  89. xf(1)=10^startf;

  90. i=startf+1;

  91. count=2;

  92. while i<=endf;

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

  94.     count=count+1;

  95.     i=i+1;

  96. end;

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

  98. 

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

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

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

  102.     min_amp=min(min_amp_th,min_amp_met);

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

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

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

  106. else

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

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

  109. end;

  110.     

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

  112. grid on;

  113. title('Bode diagram RC-circuit','Fontsize',10,'Fontweight','bold');

  114. 

  115. subplot(2,1,2);

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

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

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

  119. else

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

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

  122. end;

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

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

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

  126. 

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

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

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

  130.     min_fase=min(min_fase_th,min_fase_met);

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

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

  133.     max_fase=max(max_fase_th,max_fase_met);

  134. else

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

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

  137. end;  

  138.     

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

  140. grid on;

  141. 

  142. if nr_met~=0;

  143.     figure(2);

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

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

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

  147.     subplot(2,1,1);

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

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

  150.     hold on;

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

  152.     for i=1:nr_punt;

  153.         if sd_log_amp(i)~=0;

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

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

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

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

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

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

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

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

  162.         end;

  163.     end;

  164.     hold off;

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

  166. 

  167.     % grid bepaling

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

  169.     sizefm=size(f_m);

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

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

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

  173.     i=startfm+1;

  174.     countm=2;

  175.     flag=0;

  176.     while flag~=1;

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

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

  179.             countm=countm+1;

  180.             i=i+1;

  181.         else flag=1;

  182.         end;

  183.     end;

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

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

  186.     

  187.     fmin_th=f_th(1);

  188.     count_min=1;

  189.     while fmin_th<f_m(1);

  190.         count_min=count_min+1;

  191.         fmin_th=f_th(count_min);

  192.     end;

  193. 

  194.     size_fth=size(f_th);

  195.     size_fm=size(f_m);

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

  197.     count_max=size_fth(2);

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

  199.         count_max=count_max-1;

  200.         fmax_th=f_th(count_max);

  201.     end;

  202. 

  203.     amp_th2=amp_th(count_min:count_max);

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

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

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

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

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

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

  210. 

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

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

  213.     grid on;

  214.     title('Bode diagram RC-circuit','Fontsize',10,'Fontweight','bold');

  215. 

  216.     subplot(2,1,2);

  217.     

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

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

  220.     hold on;

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

  222.     for i=1:nr_punt;

  223.         if sd_fase(i)~=0;

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

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

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

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

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

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

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

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

  232.         end;

  233.     end;

  234.     hold off;

  235.     

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

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

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

  239.  

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

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

  242.     fase_th2=fase_th(count_min:count_max);

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

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

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

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

  247.     

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

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

  250.     grid on;

  251. end;

  252. 

  253. if nr_met~=0;

  254.     figure(3);

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

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

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

  258.     subplot(2,1,1);

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

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

  261.     hold on;

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

  263.     for i=1:nr_punt;

  264.         if sd_log_amp(i)~=0;

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

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

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

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

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

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

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

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

  273.         end;

  274.     end;

  275.     hold off;

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

  277. 

  278.     % grid bepaling

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

  280.     sizefm=size(f_m);

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

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

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

  284.     i=startfm+1;

  285.     countm=2;

  286.     flag=0;

  287.     while flag~=1;

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

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

  290.             countm=countm+1;

  291.             i=i+1;

  292.         else flag=1;

  293.         end;

  294.     end;

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

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

  297.     

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

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

  300. 

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

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

  303.     grid on;

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

  305. 

  306.     subplot(2,1,2);

  307.     

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

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

  310.     hold on;

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

  312.     for i=1:nr_punt;

  313.         if sd_fase(i)~=0;

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

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

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

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

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

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

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

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

  322.         end;

  323.     end;

  324.     hold off;

  325.     

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

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

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

  329.  

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

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

  332.  

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

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

  335.     grid on;

  336. end;