MatLab与数学建模.pptVIP

MatLab与数学建模 数理与软件工程学院 控制技术研究所 李积宪 模型？ 举例 第五教学楼大厅内的 “交大模型” 军事指挥用的 “沙盘模型” 瘟疫传播模型 模型 突出重点 忽略琐碎 用于分析、研究实际系统的运行机理 数学模型？ 观察、实践 →变量 →关系 数学公式、 逻辑准则、 具体算法 两种情况 系统结构已知 系统结构未知 举例 弹簧-阻尼-质量系统 建立模型 Assumptions: linear systems Newton’s second law of motion Viscous damper 传递函数 Laplace Transform 状态空间描述 仿 真 倒立摆 示意图 数学模型 仿真 Runge-Kutta Method Linear System Principle of superposition Homogeneity Linear Approximations Why does it work? 系统结构未知的情况 系统：Black Box 研究：观察其输入/输出特性 分析方法：时域分析/频域分析 信息可视化 t = [0 .3 .8 1.1 1.6 2.3]’; y = [0.5 0.82 1.14 1.25 1.35 1.40]’; plot(t,y,’o’), grid on Polynomial Regression X = [ones(size(t)) t t.^2]; a = X\y; T = (0:0.1:2.5)’; Y = [ones(size(T)) T T.^2]*a; plot(T,Y,–,t,y,o,), grid on Polynomial curve fitting Matlab script t = [0 .3 .8 1.1 1.6 2.3]; y = [0.5 0.82 1.14 1.25 1.35 1.40]; p2 = polyfit(t, y, 2); p3 = polyfit(t, y, 3); p4 = polyfit(t, y, 4); tt = [0:0.1:3]; y2 = polyval(p2, tt); y3 = polyval(p3, tt); y4 = polyval(p4, tt); plot(t, y, ok, tt, y2, :b, tt, y3, -.r, tt, y4, --g, linewidth, 2); grid on; title(Polynomial curve fitting); legend(raw data,2nd order,3rd order, 4h order); Linear-in-the-Parameters Regression FFT-Based Interpolation Matlab script t = [0:0.25:3]; y = [0.5 0.51 0.82 1.68 3.55 0.20 -3.30 -2.2 -1.18 0.00 0.32 0.49 0.5]; p7 = polyfit(t, y, 7); p12 = polyfit(t, y, 12); fp = interpft(y, 257); tt = [0:3/256:3]; y7 = polyval(p7, tt); y12 = polyval(p12, tt); plot(t, y, ok, tt, y7, -.r, tt, y12, :k, tt, fp, --g, linewidth, 2); grid on; title(Polynomial and fft interpolation); legend(raw data,7th order,12th order, fft); 如果允许做实验… Step Response Impulse Response 有时候“看”不出规律来 可却能够“听出”规律来 复杂系统的模型 * 分析、研究 Difficult!!!! Taylor series expansion about the operating point X = [ones(size(t)) exp(– t) t.*exp(– t)]; a = X\y; T = (0:0.1:2.5); Y = [ones(size(T)) exp(– T) T.exp(– T)]*a; plot(T,Y,–,t,y,o), grid on; Frequency Response *