Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function Marcin_Piasecki_cz3
- parametry = [1 2 3 4 5];
- while 1
- nr=input('Podaj numer wymuszenia [1, 2, 3, 4, 5]');
- As = ismember(nr,parametry);
- if As
- break;
- end
- end
- if nr==1
- wykres(1e3,1e3,0.032e-6,0.032e-6,0,10);
- elseif nr==2
- wykres(1,1.14963e4,0.032e-6,0.032e-10,0,10000);
- elseif nr==3
- wykres(1,1e3,0.032e-6,0.032e-6,0,10000);
- elseif nr==4
- wykres(1,1e4,0.032e-6,0.032e-10,0,10000);
- elseif nr==5
- wykres(1,1e5,0.032e-6,0.032e-6,0,1000);
- end
- end
- function wykres(R1,R2,C1,C2,fmin,fmax)
- n = 1000;
- G = [];
- f = [];
- h = (fmax-fmin)/n;
- f = fmin:h:fmax;
- fr = [];
- db =1;
- for i=1:n+1
- G(i) = g1(R1,R2,C1,C2,f(i));
- if G(i)>-3.047 & G(i)<-2.8
- db = i;
- end
- end
- for i=1:n+1
- fr(i) = f(db);
- end
- figure
- semilogx(f,G,fr,G);
- title('Wykres 5')
- legend('F(x)');
- xlabel("f(Hz)");
- ylabel("y(x)");
- end
- function [result] = g1(R1,R2,C1,C2,x)
- result = 20*log10(abs(1/(R1*R2*C1*C2)/((j*2*pi*x)^2+(j*2*pi*x)*(1/(R1*C1)+1/(R2*C1*R2*C2))+1/(R1*R2*C1*C2))));
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement