Đang tải... (xem toàn văn)
Mô tả chương trình mô phỏng sự tạo ảnh X quang thông qua mô phỏng sự lan truyền của photon trong mô hình vật chiếu có dạng hình học đơn giản như khối cầu, khối lập phương hoặc các kết hợp của chúng bằ.
BM Vật lý Ứng dụng, Khoa KHUD, ĐHBK TP.HCM Bài tập lớn Matlab – Vật lý A1 Bài tập 16: Vẽ quỹ đạo chuyển động ném xiên trọng trường bỏ qua lực cản xác định vài thông số liên quan” Yêu cầu Sử dụng Matlab để giải toán sau: “Hai vật ném lên từ mặt đất với vận tốc ban đầu v Coi trọng trường trái đất bỏ qua sức cản khơng khí Với góc ném α1 (góc đại lượng nhập vào lần chạy matlab) cho trước, tìm góc ném α2 để vật có tầm ném xa?” Điều kiện 1) Sinh viên cần có kiến thức lập trình MATLAB 2) Tìm hiểu lệnh Matlab liên quan symbolic đồ họa Nhiệm vụ Xây dựng chương trình Matlab: 1) Nhập giá trị ban dầu (những đại lượng đề cho) 2) Thiết lập phương trình tương ứng Sử dụng lệnh symbolic để giải hệ phương trình Xuất kết hình 3) Vẽ hình quỹ đạo vật tương ứng với góc ném Chú ý: Sinh viên dùng cách tiếp cận khác Tài liệu tham khảo: A L Garcia and C Penland, MATLAB Projects for Scientists and Engineers, Prentice Hall, Upper Saddle River, NJ, 1996 http://www.algarcia.org/fishbane/fishbane.html CÁC EM K22 KHƠNG NÊN BÊ Y NGUN VÀO NỘP NHÉ, CĨ GÌ LÀ TẤT CẢ CHÚNG TA ĐỀU BỊ BẾ ĐẤY!!! THAM KHẢO THUI NHEN HÃY CHỪA ĐƯỜNG SỐNG CHO THẾ HỆ SAU VÌ TƯƠNG LAI RA TRƯỜNG ĐÚNG HẠN 1/1 BM Vật lý Ứng dụng, Khoa KHUD, ĐHBK TP.HCM Bài tập lớn Matlab – Vật lý A1 CODE: close all; clc; clear; % Khởi tạo giá trị ban đầu v0 = 15; g = 10; a1 = 60; % Viết phương trình cần dùng syms a t vx = v0*cos(a/180*pi); vy = v0*sin(a/180*pi); x = vx*t; y = vy*t - 1/2*g*t^2; L = v0^2*sin(2*a/180*pi)/g; % Tìm góc alpha hai a2 = solve([L == subs(L, a, a1), a ~= a1, a > 0, a < 90], a); % Kết phương trình theo góc alpha thứ x1 = subs(x, a, a1); y1 = subs(y, a, a1); % Tính thời gian chuyển động t1 = single(solve(x1 == subs(L, a, a1), t)); % Kết phương trình theo góc alpha thứ hai x2 = subs(x, a, a2); y2 = subs(y, a, a2); % Tính thời gian chuyển động t2 = single(solve(x2 == subs(L, a, a1), t)); % Kiểm tra thời gian có âm hay khơng if (t1 < 0) || (t2 < 0) disp("Error") return end 1/2 ... Tính thời gian chuyển động t1 = single(solve(x1 == subs(L, a, a1), t)); % Kết phương trình theo góc alpha thứ hai x2 = subs(x, a, a2); y2 = subs(y, a, a2); % Tính thời gian chuyển động t2 = single(solve(x2