X = [ -6 -6 -7 0 7 6 6 -3 -3 0 0; ... -7 2 1 8 1 2 -7 -7 -2 -2 -7] dot2dot(X) % A 2-soros X mátrix pontjait köti össze. % Az első pontot a sorozat végén röptében megismétli, hogy % zárt alakzat legyen A1 =[1/2 0; 0 1] dot2dot(A1*X) % x irányban felére zsugorít dot2dot(forgat(X,90)) % 90 fokkal elforgat dot2dot(forgat(X,180)) for i=0:360 dot2dot(forgat(X,i)), pause(0.1), end; clear i; % visszaadjuk i komplex egységvektor szerepét!! x = 0.01:0.01:5; plot(x, fv(x)) fw = @(x) x.^2.*sin(x) + 1; s = -3:0.1:3; plot(s, fw(s)) f = inline('3*sin(2*x^2)'), argnames(f), fplot(f, [0 pi]) x = rand(1,3), y = rand(1,3), ... x(1,4) = x(1,1), y(1,4) = y(1,1), ... plot(x,y,'*-.b') x1 = linspace(-2*pi, 2*pi, 101); ... plot(x1, sin(x1), 'LineWidth', 2, 'Color', [.8 .3 0]); s = -3:0.1:3; plot(s, fw(s), 'Color', [0.7 0.1 0.9]) z = exp(i*linspace(0, 2*pi, 181)); plot(z); axis square; ... % az ábra méretezését négyzetesre alakítjuk ... h = title('Kör', 'Color', 'R'); get(h); % feliratot teszünk az aktív rajzra ... % az objektum mutatót feljegyezzük ... % megjelenítjük az objektum hangolható tulajdonságait ... set(h, 'FontSize', 20) % beállítjuk a betűméretet x = linspace(-pi/2, pi/2, 101); y = sin(x); plot(x, y, 'r'); hold on; plot(y, x, 'b'); axis auto xlabel('x tengely'), ylabel('y tengely') title('Szinusz függvény és inverze') text(-2.5, 1.7, '© Maci Laci') % szöveg, koordinátával legend('sin(x)', 'arcsin(x)', 2);% sarokpozíció megadása, 1-4 grid on, axis equal, hold off; x = [0 1; 0 0]; % 1 hosszú vektor végpontjai az x tengelyen fok = 2; T = [cosd(fok) -sind(fok); sind(fok) cosd(fok)]; % 2 fokos forgatás mátrixa for fi = 0:fok:360 plot(x(1,:),x(2,:),'.-'), axis([-1 1 -1 1]); % rajz pause(0.1); % várakozás x = T*x; % az x vektor új pozíciója (2 fokkal elfordítva) end % ciklus vége x = -10:0.1:10; figure(1), y = sin(x); subplot(1,3,1); plot(x,y); axis square, title('sin(x)'); y = sin(x)./x; subplot(1,3,2); plot(x,y); axis square, title('sin(x)/x'); y = sin(x)./(x.^2+1); subplot(1,3,3); plot(x,y); axis square, title('sin(x)/(x^2+1)'); fplot('ft', [0 pi]), grid on subplot(2,1,1); fplot('ft',[0,4], 'r'); % sima vonalas függvényrajz title('plot: t*exp(-t^2)*sin(4*t)') [x,y] = fplot('ft', [0, pi]); subplot(2,1,2); plot(x,y,'k.-'); title('fplot pontonként: t*exp(-t^2)*sin(4*t)') ezplot('ft') ezplot('sin(x)/x', [-4*pi 4*pi]) ezplot('(x-1)^2/3^2 + (y-2)^2/2^2 = 1', [-3 5, -1 5]), grid on % ellipszis ezplot('0.5*(t-2*sin(t))','0.5*(1-2*cos(t))', [-3 15]), grid on % ciklois ezplot('t*cos(t)','t*sin(t)',[0 4*pi]), grid on % archimédeszi spirális syms x; % x szimbolikus változó lesz fp = x^5-8*x^3+2 % fp függvény (szimbolikusan) fp_diff = diff(fp, x) % fp függvény deriváltja (szimbolikusan) ezplot(fp, [-2.5 2.5]), grid on, hold on % fp rajz h = ezplot(fp_diff, [-2.5 2.5]); % fp_diff rajz, objektumkezelő feljegyzés set(h, 'Color', 'g') % objektum színének változtatása axis([-2.5 2.5 -40 40]) % intervallumhatárok beállítása title('x^5 - 8*x^3 + 2 és deriváltja 5*x^4 - 24*x^2') legend('x^5 - 8*x^3 + 2', '5*x^4 - 24*x^2', 1) % cím, jelmagyarázat (jobb felső sarok) subplot(1,2,1), ezmeshc('y/(1 + x^2 + y^2)', [-5,5,-2*pi,2*pi]) subplot(1,2,2), s='sqrt(1-z^2)*'; h=ezsurf([s 'cos(t)'], [s 'sin(t)'] ,'z', [0 2*pi, -1 1]) % az s karakterlánc sorvektorban tárolódik, így az % [s 'cos(t)'] sorvektor a bővített karakterláncot jelenti subplot(1,2,1), ezpolar('1+cos(t)'), subplot(1,2,2), ezpolar('1+cos(4*t)') ezplot3('t*sin(2*t)', 't*cos(2*t)', '6*pi-t',[0,6*pi], 'animate') fplot('log10', [0.1 2*pi]), hold on, fplot('cos', [0.1 2*pi], 'k'), grid on, title('log10(x) és cos(x) metszéspontjai:') legend('log10(x)','cos(x)',4) % jobb alsóba jelmagyarázat fd = @(x) log10(x)-cos(x) % különbségfüggvény x = [fzero(fd, 2) fzero(@(x) log10(x)-cos(x), 6)] % zérushelyek keresése plot(x, fd(x), '.b') % kék pontok: a különbségfüggvény zérushelyei plot(x, cos(x), 'or') % piros körök: metszéspontok text(x(1)-0.5,-0.1,sprintf('%f',x(1))), text(x(2)-0.5,-0.1,num2str(x(2))) xd = x; % zérushelyek feljegyzése az integráláshoz fplot('f1', [-2.5 2.5]), grid on, hold on [mnx mny kod] = ... fminbnd('f1', 2, 2.5, optimset('TolX', eps)), plot(mnx,mny,'ob') % az optimset beállítás az fzero-nál is használható mxx = fminbnd(@(x) -f1(x), -3, -2), mxy=f1(mxx), plot(mxx, mxy, 'or') hx = text(mxx+.2, mxy, sprintf('(%.2f, %.2f)', mxx, mxy)); set(hx,'fontweight','b','color','r') hn = text(mnx-.4, mny-3, sprintf('(%.2f, %.2f)', mnx, mny)); set(hn, 'color', 'b') % kék koordináták title('f_1(x) = x^5-8*x^3+2 szélsőértékei') % határozott integrálok: quad('f1', mxx, mnx) % saját függvény m-fájllal megadva quad('sin',0,pi,eps) % elemi függvény quad('sin(x)+cos(x/2)',0,pi) % képlettel megadott fv. quad(@(x) log10(x)-cos(x), xd(1), xd(2)) % itt definiált anonymus fv. quad(fd, xd(1), xd(2)) % korábban definiált anonymus fv.