TH Tin học kỹ thuật matlab 1

Tin học kỹ thuật Cơ sở MatLab Thiết kế giao diện Simulink Symbolic Communication System / DSP System / Simscape / Control System Bài 1 BÀI 1: CƠ SỞ VỀ MATLAB Khai báo biến, phép toán toán tử Các lệnh thường dùng Các phát biểu Ma trận Script function Các hàm vẽ Bài Khai báo biến, phép toán toán tử Khai báo biến - Phân biệt chữ hoa chữ thường - Không cần phải khai báo kiểu biến - Tên biến phải bắt đầu ký tự khơng có khoảng trắng - Không đặt tên trùng với tên đặc biệt MatLab - Để khai báo biến toàn cục (sử dụng tất chương trình con), phải dùng thêm từ khố global phía trước Bài Khai báo biến, phép toán toán tử Các phép toán toán tử Các hằng: + , - , * , / , \ (chia trái) , ^ (mũ) , ' (chuyển vị hay số phức liên hiệp) pi 3.14159265 i số ảo Các toán tử quan hệ : j tương tự i eps sai số 2-52 realmin số thực nhỏ = 2-1022 realmax số thực lớn = 21023 inf vô lớn NaN Not a number < , , >= , == , ~= Các toán tử logic : & , | (or) , ~ (not) Bài Khai báo biến, phép toán toán tử >>a= 1; >>b=-2; >>c=1; >>delta= b^2- 4*a*c; % Tinh delta >>x1=(-b+ sqrt(delta) )/(2*a) >>x2=(-b- sqrt(delta) )/(2*a) Bài Các hàm toán học cos(x) tính cos sin(x) tính sin tan(x) tính tang acos(x) tính cos-1 exp(x) hàm mũ số e asin(x) tính sin-1 sqrt(x) bậc hai x atan(x) tính tang-1 𝑒 𝑥 +𝑒 −𝑥 log(x) logarit số e cosh(x) tính log10(x) logarit số 10 coth(x) tính cosh(x)/sinh(x) abs(x) module số phức x (giá trị tuyệt đối số thực) angle(x) argument số phức a sinh(x) tính 𝑒 𝑥 −𝑒 −𝑥 tanh(x) tính sinh(x)/cosh(x) conj(x) số phức liên hợp x acosh(x) tính cosh-1 imag(x) phần ảo x acoth(x) tính coth-1 real(x) phần thực x asinh(x) tính sinh-1 sign(x) dấu x Bài atanh(x) tính tanh-1 Các lệnh thường dùng >>help tên_hàm % tham khảo help hàm >>lookfor ‘chuỗi’ %Tìm kiếm >>clc % Xố hình >>clear tên_biến % Xố biến >>clear all %Xoá tất biến >> clf %Xoá figure >>save % >>load % >>who % liệt kê biến nhớ >>whos % liệt kê chi tiết biến nhớ >>which >>what Bài Các phát biểu if if end elseif if else else end end Bài Các phát biểu Ví dụ: Xét chương trình đốn tuổi sau: disp(‘Xin chao! Han hanh duoc lam quen’); %Xuất chuỗi hình x = fix(30*rand); disp(‘Tuoi toi khoang ‐ 30’); g = input(‘Xin nhap tuoi cua ban: ‘); if g < x disp(‘Ban tre hon toi’); elseif g > x disp(‘Ban lon hon toi’); else disp(‘Ban bang tuoi toi’); end Bài Các phát biểu switch case n1 case n2 case nn Otherwise end Bài 10 Định dạng LaTex Phân số: \frac{tử số}{mẫu số} Bài 40 Định dạng LaTex Dấu căn: \sqrt{} hay \sqrt[]{} \sqrt{a^2+b^2} \qquad \sqrt[n]{a^n+b^n} Tích phân: \int Tổng: \sum \int^a_bx^2dx \sum_{i=-\infty}^{\infty}\int^a_bx^2dx Bài 41 Các dấu ngoặc: \left \right, sau dấu ngoặc (, ), |, [, ] Ma trận: \matrix (phân biệt phần tử hàng dấu & xuống hàng \cr) \left| \matrix{1 & a &c \cr & \cr c & & d} \right| \left( \matrix{\frac{a}{b} & a &c \cr & \cr c & & d} \right) Bài 42 Các ký tự khác: Bài 43 Bài 44 Bài 45 Bài 46 Bài 47 ∞ න 𝑒 −𝑥 + 𝑒 𝑥 2𝑥 + −∞ sin 𝜔 𝑥 − 𝜏 𝑑𝑥 ['$$\int_{-\infty}^{+\infty}\frac{e^{x}+e^x}{\sqrt{2x+1}}sin\left(\omega(x - \tau)\right)dx$$'] 𝑥 𝑦ሷ = 𝑥 𝑥 𝑥3 𝑥2 𝑥3 𝑥 𝑥3 𝑥 𝑥2 ['$$\ddot y = x\left[{\matrix{x & x^2 & x^3 \cr x^2 & x^3 & x \cr x^3 & x & x^2}}\right]$$'] ሶ ∆𝑍𝑒 + ሶ 𝑘 𝑍𝑒 ሶ =− 𝐽Ԧ𝑛𝑔 𝑗𝜔𝜀 ['$$\Delta\dot {\vec{Z_e}} + k^2\dot {\vec{Z_e}} = \frac{1}{j\omega\epsilon}\dot{\vec{J}}_{ng}$$'] Bài 48 𝐸= 𝑄𝑧 4𝜋𝜀 𝑎2 + 𝑧2 𝑖𝑧 N −1 X (k ) =  x(n)e − j 2kn / N n =0 Bài 49 t = 0:0.005:4*pi; fm =100; ft = 10;V0 = 5; vam = V0*(cos(ft*t)).*cos(fm*t); subplot(2,1,1); plot(t,vam) axis([0 -10 10]); subplot(2,1,2); pwelch(vam); 𝑉𝐴𝑀 𝑡 = 𝑉0 𝑐𝑜𝑠Ω𝑡𝑐𝑜𝑠𝜔0 𝑡 10 -5 -10 0.5 1.5 2.5 3.5 Welch Power Spectral Density Estimate Power/frequency (dB/rad/sample) 40 20 -20 -40 -60 0.1 0.2 0.3 0.4 0.5 0.6 Normalized Frequency ( rad/sample) 0.7 0.8 0.9 501 𝑉𝐴𝑀 𝑡 = 𝑉0 + 𝑚𝑐𝑜𝑠Ω𝑡 𝑐𝑜𝑠𝜔0 𝑡 t = 0:0.001:2*pi; fm = 200;m = 0.5; ft = 40;V0 = 5; vam = V0*(1+m*cos(ft*t)).*cos(fm*t); subplot(2,1,1); plot(t,vam) axis([0 -10 10]); subplot(2,1,2); pwelch(vam); 10 -5 -10 0.2 0.4 0.6 0.8 1.2 1.4 1.6 1.8 0.7 0.8 0.9 Welch Power Spectral Density Estimate Power/frequency (dB/rad/sample) 40 20 -20 -40 -60 0.1 0.2 0.3 0.4 0.5 0.6 Normalized Frequency ( rad/sample) 51 t = 0:.001:2*pi; y = square(2*pi*t);y=(1+y)/2; xc = sin(2*pi*5*t); xask = xc.*y; subplot(2,1,1); plot(t,xask) axis([0 -2 2]); subplot(2,1,2); pwelch(xask); Tín hiệu nhị phân (đơn cực) -1 -2 0.5 1.5 2.5 3.5 Welch Power Spectral Density Estimate Power/frequency (dB/rad/sample) 20 -20 -40 -60 -80 0.1 0.2 0.3 0.4 0.5 0.6 Normalized Frequency ( rad/sample) 0.7 0.8 0.9 52 y2 = (1-y)/2; m(t) + Tín hiệu nhị phân xc2 = sin(2*pi*5*t); BFSK (đơn cực) x2 = xc2.*y2; Acos 2𝜋𝑓1 𝑡 xfsk = x1+x2; + ഥ Đảo 𝑚(𝑡) subplot(2,1,1); bit plot(t,xfsk) Acos 2𝜋𝑓2 𝑡 axis([0 -2 2]); t = 0:.001:2*pi; subplot(2,1,2); y = square(2*pi*t); pwelch(xfsk); y1=(1+y)/2; xc1 = sin(2*pi*10*t); x1 = xc1.*y1; -1 -2 0.5 1.5 2.5 3.5 Welch Power Spectral Density Estimate Power/frequency (dB/rad/sample) 20 -20 -40 -60 -80 0.1 0.2 0.3 0.4 0.5 0.6 Normalized Frequency ( rad/sample) 0.7 0.8 0.9 53 Tín hiệu nhị phân (lưỡng cực) t = 0:.001:2*pi; y = square(2*pi*t); xc = sin(2*pi*10*t); xpsk = xc.*y; subplot(2,1,1); plot(t,xpsk) axis([0 -2 2]); subplot(2,1,2); pwelch(xpsk); BPSK Acos 2𝜋𝑓𝑐 𝑡 -1 -2 0.5 1.5 2.5 3.5 Welch Power Spectral Density Estimate Power/frequency (dB/rad/sample) 20 -20 -40 -60 -80 0.1 0.2 0.3 0.4 0.5 0.6 Normalized Frequency ( rad/sample) 0.7 0.8 0.9 54