Tổng hợp code các bài tập thực hành Matlab

Đây là các code một số bài tập mâtb cho dẫn vật lý, rất hay và dễ hiểu mà mình sưa tầm được...
Nội dung trích xuất từ tài liệu:
Tổng hợp code các bài tập thực hành Matlab Cao V¨n ChÝnh_B1K52 Tổng hợp code các bài thực hành MatlabBai 1clcx=-2:0.01:2y=exp(-0.3.*x.^2- 0.5.*x +1).*cos(20.*x)subplot(2,1,1)% lenh subplot nay chia man%hinh do thi theo 2 hang 1 cot vi tri 1plot(x,y,m)%lenh nay ve do thi vao phan %hinh minh chia ra vào subplot 1%hinhcon 2subplot(2,1,2)% ve dang ham nay thi khai bao truc tiep ham trong dau fplot(exp(-0.3.*x.^2- 0.5.*x+1).*cos(20.*x),[-2 2])title(do thi ham so )% lenh nay dien%tieu de cho do thixlabel(truc x)% ghi tieu de cho truc xxlabel(truc y)% ghi tieu de cho truc yBai 2%hinh con 1x=0:0.05:4;y=3.*sin(pi.*x);y1=exp(-0.2.*x)subplot(2,1,1);plot(x,y,r,x,y1,b)gtext(giao diem 1)% gan chu thich vao%cac giao diem giua 2 do thigtext(giao diem 2)gtext(giao diem 3)gtext(giao diem 4)% hinh con 2& 3x1=0:0.05:10;y2=exp(-0.2.*x1).*cos(x1)subplot(2,2,3);bar(x1,y2)% ve do thi dang cotsubplot(2,2,4);stairs(x1,y2)%ve do thi dang bac thang--------------------------------------------------------------------------------------------Bai 3% Chu y ham function o day chi chay dc tren Window thui,voi may tinh tren truong thi khongchay duoc dau!!!%cau afunction tong=bai3a(n,x)s=0;for k=1:n s=s+((-1).^k./k).*((-(x-1).^k).*(x>0 & x=1));endtong=s% bai3bfunction tong=bai3b(x)ak=1;ss=1e-6; s1=0;k=0;while any(abs(ak)>ss) k=k+1 ak=((-1).^k./k).*((-(x-1).^k).*(x>0 & x=1)); s1=s1+ak;end tong=s1%bai 3cx=0.1:0.1:10;plot(x,bai3a(1,x),r);hold on;plot(x,bai3a(5,x),c);hold on;plot(x,bai3a(10,x),b);hold on;plot(x,bai3b(x),m);hold on;legend(s1,s5,s10,s)gtext(s1) ;gtext(s5)gtext(s10) ;gtext(s)Bai 4 Cao V¨n ChÝnh_B1K52[x,y]=meshgrid(-2:.1:2);z=2.*x.*y./(x.^2 +y.^2 +1);subplot(2,2,1);mesh(x,y,z);% do thi meshtitle(mesh)% gan chu thich meshsubplot(2,2,2);surf(x,y,z)% do thi surftitle(surf)subplot(2,2,3);contour(x,y,z)title(contour)%do thi duong dong muc contoursubplot(2,2,4);contour3(x,y,z)title(contour3)%do thi duong dong muc 3chieuBai 5% Cau asubplot(2,2,1);[teta,phi]=meshgrid(0:pi/30:pi, 0:pi/30:2*pi);x0=0;y0=0;z01=-1;z02=1;r=1;x1=x0+r*cos(phi).*sin(teta);y1=x0+r*sin(phi).*sin(teta);z11=z01+r*cos(teta);z12=z02+r*cos(teta);mesh(x1,y1,z11);hold onmesh(x1,y1,z12);% Cau bsubplot(2,2,2);R=2;r=0.25;[teta,phi]=meshgrid(0:pi/30:2*pi);x=(R+r*cos(teta)).*cos(phi)y=(R+r*cos(teta)).*sin(phi)z=r*sin(teta);mesh(x,y,z);% Cau csubplot(2,1,2);a1=[0,1/sqrt(3),-1/2*sqrt(3)];b1=[0,0,1/2];c1=[sqrt(2/3),0,0];fill3(a1,b1,c1,r);hold on;a2=[1/sqrt(3),-1/2*sqrt(3),-1/2*sqrt(3)];b2=[0,1/2,-1/2];c2=[0,0,0];fill3(a2,b2,c2,y);a3=[-1/2*sqrt(3),-1/2*sqrt(3),0];b3=[1/2,-1/2,0];c3=[0,0,sqrt(2/3)];fill3(a3,b3,c3,g);a4=[-1/2*sqrt(3),0,1/sqrt(3)];b4=[-1/2,0,0];c4=[0,sqrt(2/3),0];fill3(a4,b4,c4,r);Bai 6:function t=bai6(n,x)if (n Cao V¨n ChÝnh_B1K52plot(x,bai6(4,x));ylim([-1.5 1.5])%gioi han do thi[x0,y0]=ginput;% lay toa do = cach click chuotf=inline(bai6(4,x),x);% de y vao ham nay khi tim khong diemfor k=1:length(x0);%do rong cua khong diem [xn(k) fxn(k)]=fzero(f,x0(k));%ham tim khong diem thuc fzero text(xn(k),fxn(k),< 0 diem)%end% Tim cac max[x1,y1]=ginput;f1=inline((-1)*bai6(4,x),x);for i=1:length(x1); [xn(i) fxn(i)]=fminsearch(f1,x1(i)); text(xn(i),(-1)*fxn(i),< max)%chu y co -1end% Tim cac min[x1,y1]=ginput;f1=inline(bai6(4,x),x);for i=1:length(x1); [xn(i) fxn(i)]=fminsearch(f1,x1(i));% ham tim cuc tieu text(xn(i),fxn(i),< min)endBai 7function [x0,fx0,n,xn,fxn]=bai7(a,b,fx,dfx,tol);%tol la nhap so lieu dang sai so%thay vi lap dao ham cua fx ta tinh luon dfxfx=inline(fx);dfx=inline(dfx);%inline la%1 ham moi ma minh tu dat no dang nhu 1%function hoac @ hamxi=linspace(a,b,100);yi=subs(fx,x,xi);%ham dao bienplot(xi,yi,r);grid on[x0,fx0]=ginput(5);for i=1:length(x0) xn(i)=x0(i);%dau ra xn n(i)=0; k(i)=fx(x0(i))./dfx(x0(i));%chu y chuyen ve while abs(k(i)) Cao V¨n ChÝnh_B1K52figure(2)[x,y]=meshgrid(-2:.2:2); z=x+i*y;%z la bieu thuc phucf=z.^2+cos(log(abs(z)+1)+3)-2.^z-i;contour(x,y,real(f),[0 0],r);hold on% ve theo phan thuccontour(x,y,imag(f),[0 0],b);hold off% ve theo phan ao[x0,y0]=ginput; z0=x0+i*y0;fz=inline(z^2+cos(log(abs(z)+1)+3)-2^z-i,z);for k=1:length(z0) [zn,fzn]=fsolve(fz,z0(k))endfprintf( Cac cap nghiem cua phuong trinh 2 la: );znfznBai 9% phan aclear allR=1r=R.*sin(pi/10)./cos(pi/5)k=0:9;x=R.*sin(k.*pi./5).*((-1).^k==1)+r.*sin(k.*pi./5).*((-1).^k==-1);y=R.*cos(k.*pi./5).*((-1).^k==1)+r.*cos(k.*pi./5).*((-1).^k==-1);disp(Toa do d ...