MATLAB 07111014第5次上机作业 曹红兴
matlab上机实验答案-整理版
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y)subplot(2,1,1)plot(x,y)subplot(2,1,2)plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi;rho=a*cos(b+n*theta);polar(theta,rho)3. 列出求下列空间曲面交线的程序clcclearx=[-5:0.5:5];[X,Y]=meshgrid(x);z1=X.^2-2*Y.^2;z2=X.*2-Y.*3;xlabel('x')ylabel('y')zlabel('z')surf(X,Y,z1)hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5);x1=X(k)y1=Y(k)z3=x1.^2-2*y1.^2hold onplot3(x1,y1,z3,'*')⎪⎭⎫ ⎝⎛+⋅=-334sin 234πt e y t ()θρn b a +=cos 2212y x z -=y x z 322-=4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clcclearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2));plot(x,y,'b*-');title('绘图');xlabel('x 坐标');ylabel('y 坐标');legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y x clccleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8];c=b/a;x=c(1,1)y=c(1,2)z=c(1,3)w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
MATLAB上机实习指导书
MATLAB上机实习指导书目录实习内容一MATLAB基本操作-----------------------------------2 实习内容二数值数组及其运算---------------------------------8 实习内容三MATLAB图形绘制基础-------------------------------15 实习内容四SIMULINK仿真基础---------------------------------22 附录1 流程控制------------------------------------------------24 附录2 SIMULINK的库模块-------------------------------------- 25附录3 转义符号------------------------------------------------27实习内容一MATLAB基本操作一、实习目的掌握MATLAB的启动和退出；熟悉MATLAB的命令窗口；熟悉其它常用窗口和工具栏。
二、实习内容与步骤1、进入MATLAB的开发环境。
方法一：点击桌面上的快捷方式或matlab\文件夹下的快捷方式图标。
方法二：开始→所有程序→MA TLAB6.5→MATLAB6.5.EXE。
方法三：点击matlab\bin\win32文件夹下的matlab.exe。
这三种方法的当前目录不同。
优先选用方法一。
可见到如下交互界面（见图1-1）：图1-1 Matlab默认桌名平台2、退出MATLAB的开发环境。
退出Matlab的方式很多：（1）在Matlab命令窗口的“File”菜单下选择“Exit Matlab”；（2）快捷键“Ctrl+q”；（3）在命令窗口输入“quit”命令；（4）在命令窗口输入“exit”命令；（5）用鼠标单击命令窗口右上角的“X”按钮；（6）用鼠标双击命令窗口左上角的图标“”。
MATLAB上机内容及作业
MATLAB上机内容及作业无约束优化求解函数fminsearch和fminunc求解无约束非线性优化问题的函数有fminsearch 函数、fminunc 函数。
函数fminsearch和fminunc功能相同，但fminunc函数可以得到目标函数在最优解处的梯度和Hessian矩阵值。
无约束优化数学模型为：min f(X) X∈R n求解无约束非线性优化问题的步骤为：第一步：先编写目标函数的M文件；第二步：再在命令窗口中调用相应的优化函数。
1、fminsearch函数调用格式为[x, fval]=fminsearch(@myfun, x0)输出参数的含义：x：返回最优解的设计变量的值；fval：在最优设计变量值时，目标函数的最小值；exitflag：返回算法终止的标志，有以下几种情况，＞0 表示算法收敛于最优解处；＝0 表示算法已经达到迭代的最大次数；＜0 表示算法不收敛。
output：返回优化算法信息的一个数据结构，有以下信息：output.iteration 表示迭代次数output.algorithm 表示所采用的算法output.funcCount 表示函数评价次数输入参数的含义：@myfun：目标函数的M文件，在其前要加“@”，或表示为'myfun' ，myfun自己可以任意命名；x0：在调用该优化函数时，需要先对设计变量赋一个初始值；2、fminunc函数的调用格式[x, fval]=fminunc (@myfun, x0)grad：返回目标函数在最优解处的梯度信息；hessian：返回目标函数在最优解处的hessian矩阵信息。
其余含义同上。
3、实例已知某一优化问题的数学模型为：min f(X)=3x12+2x1x2+x22X∈R n用MA TLAB程序编写的代码为：第一步：首先编写目标函数的.m文件，并保存为examplefsearch.m文件（先单击file菜单，后点击New 命令中的M—file，即可打开M文件编辑窗口进行代码的编辑，在英文状态下输入程序代码），代码为：function f=examplefsearch（x）f=3*x(1)^ 2+2*x(1)*x(2)+x(2)^2;第二步：在Command窗口中调用fminsearch函数，代码为：x0=[1;1]; %赋初值[x,fval]=fminsearch(@examplefsearch,x0) %回车即可调用fminsearch函数，得到结果输出最优解结果为：x=1.0e-0.08* -0.7914 0.2260 %分别为x1和x2的最优点的值（近似为0）fval=1.5722e-016 %对应最优点的最优目标函数值（近似为0）4、作业已知几个优化问题的数学模型分别为：（1）min f(X)=0.1935x1 x22 x32(4+6x4) X∈R4（2）min f(X)= (x13+cos x2+log x3）/ e x1 X∈R3（3）min f(X)=2x13+4x1x23 -10x1x2+x33X∈R3试用MATLAB编程分别求解上述优化问题的最优解。
Matlab上机作业部分参考答案
-2
-1.7
-1.4
-1.1
-0.8
-0.5
-0.2
0.1
0.4
பைடு நூலகம்0.7
1
1.3
.10289 .11741 .13158 .14483 .15656 .16622 .17332 .1775
.17853 .17635 .17109 .16302
【求解】用下面的语句可以立即得出给定样本点数据的三次插值与样条 插值，得出的结果如，可见，用两种插值方法对此例得出的结果几乎一 致，效果均很理想。 >> x=[-2,-1.7,-1.4,-1.1,-0.8,-0.5,-0.2,0.1,0.4,0.7,1,1.3,... 1.6,1.9,2.2,2.5,2.8,3.1,3.4,3.7,4,4.3,4.6,4.9]; y=[0.10289,0.11741,0.13158,0.14483,0.15656,0.16622,0.17332,... 0.1775,0.17853,0.17635,0.17109,0.16302,0.15255,0.1402,... 0.12655,0.11219,0.09768,0.08353,0.07019,0.05786,0.04687,... 0.03729,0.02914,0.02236]; x0=-2:0.02:4.9; y1=interp1(x,y,x0,'cubic'); y2=interp1(x,y,x0,'spline'); plot(x0,y1,':',x0,y2,x,y,'o')
S 2i 1 2 4 8 262 263 用数值方法可以求出，
i 0
63
试不采用循环的形式求出和式的数值解。由于数值方法采用 double 形式进行计算的，难以保证有效位数字，所以结果不 一定精确。试采用符号运算的方法求该和式的精确值。 【求解】用符号运算的方式可以采用下面语句 >> sum(2.^[1:63]) ans = 1.844674407370955e+019 >> sum(sym(2).^[1:63]) ans = 18446744073709551614
MATLAB课后上机实验报告
《MATLAB程序设计与应用》实验指导书实验一 matlab 集成环境使用与运算基础1，先求下列表达式的值，然后显示matlab 工作空间的使用情况并保存全部变量。
（1）0122sin851z e =+程序：.>> z1=2*sin(85*pi/180)/(1+exp(2)) 结果： z1 =0.2375（2）22121(0.4552i z In x x +⎡⎤=+=⎢⎥-⎣⎦其中 程序：>> x=[2,1+2*i;-0.45,5];>> z2=0.5*log(x+sqrt(1+x*x)) 结果： z2 =0.7114 - 0.0253i 0.8968 + 0.3658i 0.2139 + 0.9343i 1.1541 - 0.0044i（3）0.3,9.2,8.2,...,8.2,9.2,0.3,23.0)3.0sin(23.03.03---=+++-=-a aIn a e e z a a 提示：利用冒号表达式生成a 向量，求各点函数值时用点乘运算。
程序：>> a=-3.0:0.1:30;>> z3=(exp(0.3*a)-exp(-0.3*a))/2.*sin((a+0.3)*pi/180)+log((0.3+a)/2) 结果： z3 =1.0e+003 *Columns 1 through 40.0003 + 0.0031i 0.0003 + 0.0031i 0.0003 + 0.0031i 0.0002 + 0.0031iColumns 5 through 80.0002 + 0.0031i 0.0001 + 0.0031i 0.0001 + 0.0031i 0.0000 + 0.0031i Columns 9 through 12-0.0000 + 0.0031i -0.0001 + 0.0031i -0.0001 + 0.0031i -0.0002 + 0.0031i Columns 13 through 16-0.0003 + 0.0031i -0.0003 + 0.0031i -0.0004 + 0.0031i -0.0005 + 0.0031i Columns 17 through 20-0.0006 + 0.0031i -0.0007 + 0.0031i -0.0008 + 0.0031i -0.0009 + 0.0031i Columns 21 through 24-0.0010 + 0.0031i -0.0012 + 0.0031i -0.0014 + 0.0031i -0.0016 + 0.0031i Columns 25 through 28-0.0019 + 0.0031i -0.0023 + 0.0031i -0.0030 + 0.0031i -0.0370 Columns 29 through 32-0.0030 -0.0023 -0.0019 -0.0016 Columns 33 through 36-0.0014 -0.0012 -0.0010 -0.0009 Columns 37 through 40-0.0008 -0.0007 -0.0006 -0.0005 Columns 41 through 44-0.0004 -0.0003 -0.0003 -0.0002 Columns 45 through 48-0.0001 -0.0001 -0.0000 0.0000 Columns 49 through 520.0001 0.0001 0.0002 0.0002 Columns 53 through 560.0003 0.0003 0.0003 0.0004 Columns 57 through 600.0004 0.0005 0.0005 0.0005 Columns 61 through 640.0006 0.0006 0.0006 0.0007 Columns 65 through 680.0007 0.0007 0.0008 0.0008 Columns 69 through 720.0008 0.0008 0.0009 0.0009 Columns 73 through 760.0009 0.0010 0.0010 0.0010 Columns 77 through 800.0011 0.0011 0.0011 0.0011 Columns 81 through 840.0012 0.0012 0.0012 0.0013 Columns 85 through 880.0013 0.0013 0.0013 0.0014 Columns 89 through 920.0014 0.0014 0.0015 0.0015 Columns 93 through 960.0015 0.0016 0.0016 0.0016 Columns 97 through 1000.0017 0.0017 0.0017 0.0018 Columns 101 through 1040.0018 0.0018 0.0019 0.0019 Columns 105 through 1080.0020 0.0020 0.0020 0.0021 Columns 109 through 1120.0021 0.0022 0.0022 0.0023 Columns 113 through 1160.0023 0.0024 0.0024 0.0025 Columns 117 through 1200.0025 0.0026 0.0026 0.0027 Columns 121 through 1240.0027 0.0028 0.0029 0.0029 Columns 125 through 1280.0030 0.0031 0.0031 0.0032 Columns 129 through 1320.0033 0.0034 0.0034 0.0035 Columns 133 through 1360.0036 0.0037 0.0038 0.0039 Columns 137 through 1400.0040 0.0041 0.0042 0.0043 Columns 141 through 1440.0044 0.0045 0.0046 0.0047 Columns 145 through 1480.0049 0.0050 0.0051 0.0053 Columns 149 through 1520.0054 0.0056 0.0057 0.0059 Columns 153 through 1560.0060 0.0062 0.0064 0.0066 Columns 157 through 1600.0068 0.0069 0.0071 0.0074 Columns 161 through 1640.0076 0.0078 0.0080 0.0083 Columns 165 through 1680.0085 0.0088 0.0090 0.0093 Columns 169 through 1720.0096 0.0099 0.0102 0.0105 Columns 173 through 1760.0108 0.0112 0.0115 0.0119 Columns 177 through 1800.0123 0.0127 0.0131 0.0135 Columns 181 through 1840.0139 0.0144 0.0148 0.0153 Columns 185 through 1880.0158 0.0163 0.0168 0.0174 Columns 189 through 1920.0180 0.0185 0.0191 0.0198 Columns 193 through 1960.0204 0.0211 0.0218 0.0225 Columns 197 through 2000.0233 0.0241 0.0249 0.0257 Columns 201 through 2040.0265 0.0274 0.0284 0.0293 Columns 205 through 2080.0303 0.0313 0.0324 0.0335 Columns 209 through 2120.0346 0.0358 0.0370 0.0382 Columns 213 through 2160.0395 0.0409 0.0423 0.0437 Columns 217 through 2200.0452 0.0467 0.0483 0.0500 Columns 221 through 2240.0517 0.0534 0.0552 0.0571 Columns 225 through 2280.0591 0.0611 0.0632 0.0654 Columns 229 through 2320.0676 0.0699 0.0723 0.0748 Columns 233 through 2360.0773 0.0800 0.0827 0.0856 Columns 237 through 2400.0885 0.0915 0.0947 0.0979 Columns 241 through 2440.1013 0.1047 0.1083 0.1121 Columns 245 through 2480.1159 0.1199 0.1240 0.1282 Columns 249 through 2520.1326 0.1372 0.1419 0.1467 Columns 253 through 2560.1518 0.1570 0.1624 0.1679 Columns 257 through 2600.1737 0.1796 0.1858 0.1921 Columns 261 through 2640.1987 0.2055 0.2125 0.2198 Columns 265 through 2680.2273 0.2351 0.2431 0.2514 Columns 269 through 2720.2600 0.2689 0.2781 0.2876 Columns 273 through 2760.2974 0.3076 0.3180 0.3289 Columns 277 through 2800.3401 0.3517 0.3637 0.3761 Columns 281 through 2840.3889 0.4021 0.4158 0.4299 Columns 285 through 2880.4446 0.4597 0.4753 0.4915 Columns 289 through 2920.5082 0.5254 0.5433 0.5617 Columns 293 through 2960.5807 0.6004 0.6208 0.6418 Columns 297 through 3000.6636 0.6861 0.7093 0.7333 Columns 301 through 3040.7581 0.7838 0.8103 0.8377 Columns 305 through 3080.8660 0.8952 0.9254 0.9567 Columns 309 through 3120.9890 1.0223 1.0568 1.0924 Columns 313 through 3161.1292 1.1673 1.2066 1.2472Columns 317 through 3201.2892 1.3326 1.3774 1.4237Columns 321 through 3241.4715 1.5210 1.5721 1.6249Columns 325 through 3281.6794 1.7357 1.7940 1.8541Columns 329 through 3311.9163 1.98052.0468（4）⎪⎩⎪⎨⎧=<≤<≤<≤+--=5.2:5.0:0,322110,121,2224t t t t t t t t z 其中提示：用逻辑表达式求分段函数值。
MATLAB上机实验实验报告
MATLAB上机实验一一、实验目的初步熟悉MATLAB 工作环境，熟悉命令窗口，学会使用帮助窗口查找帮助信息。
命令窗口二、实验内容(1) 熟悉MATLAB 平台的工作环境。
(2) 熟悉MATLAB 的5 个工作窗口。
(3) MATLAB 的优先搜索顺序。
三、实验步骤1. 熟悉MATLAB 的5 个基本窗口①Command Window (命令窗口)②Workspace (工作空间窗口)—③Command History (命令历史记录窗口)④Current Directory (当前目录窗口)⑤Help Window (帮助窗口)(1) 命令窗口(Command Window)。
在命令窗口中依次输入以下命令：>>x=1>> y=[1 2 34 5 67 8 9];>> z1=[1:10],z2=[1:2:5];>> w=linspace(1,10,10);>> t1=ones(3),t2=ones(1,3),t3=ones(3,1)>> t4=ones(3),t4=eye(4)x =1z1 =1 2 3 4 5 6 7 8 9 10 t1 =1 1 11 1 11 1t2 =1 1 1t3 =111t4 =1 1 11 1 11 1 1t4 =1 0 0 00 1 0 00 0 1 00 0 0 1思考题：①变量如何声明，变量名须遵守什么规则、是否区分大小写。
答：（1）变量声明1.局部变量每个函数都有自己的局部变量，这些变量只能在定义它的函数内部使用。
当函数运行时，局部变量保存在函数的工作空间中，一旦函数退出，这些局部变量将不复存在。
脚本（没有输入输出参数，由一系列MATLAB命令组成的M文件）没有单独的工作空间，只能共享调用者的工作空间。
当从命令行调用，脚本变量存在基本工作空间中；当从函数调用，脚本变量存在函数空间中。
2.全局变量在函数或基本工作空间内，用global声明的变量为全局变量。
matlab上机习题详解试题答案
百度文库 - 让每个人平等地提升自我P 第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
MATLAB第一次 07111014
数学软件与实验第一次上机作业上机时间：2013-4-10 地点：E204班级:071111 学号:07111014 姓名:曹红兴1.计算三角形三边分别为a，b，c中c边对应内角的角度>> a = 3; b = 3; c = 3;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =60.0000>> a = 3; b = 4; c = 5;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =90>> a = 3; b = 4; c = 20;>> acos((a^2+b^2-c^2)/2/a/b) *180/pians =1.8000e+002 -1.9715e+002i2.试分别生成5 阶的单位阵、8 阶均匀分布的随机矩阵及其下三角矩阵，要求矩阵元素为介于10~99之间整数>> M=eye(5)M =1 0 0 0 00 1 0 0 00 0 1 0 00 0 0 1 00 0 0 0 1>> N=randsrc(8,8,[10:99])N =83 96 47 71 34 49 73 9691 96 92 78 14 44 77 4021 24 81 76 18 78 34 6292 97 96 45 84 81 71 3066 96 69 68 72 26 68 7718 53 13 25 38 54 24 3235 82 86 73 95 50 20 5559 22 94 12 13 68 54 72>> Z=tril(N)Z =83 0 0 0 0 0 0 091 96 0 0 0 0 0 021 24 81 0 0 0 0 092 97 96 45 0 0 0 066 96 69 68 72 0 0 018 53 13 25 38 54 0 035 82 86 73 95 50 20 059 22 94 12 13 68 54 723.生产列向量x=[1, 3, 10, 15, 20, 25, 30, 35, 40]>> x=[1;3;10;15;20;25;30;35;40]x =13101520253035404.生成以x 的元素为对角线的矩阵A，并输出A 的行数>> A=diag(x)a=length(A)A =1 0 0 0 0 0 0 0 00 3 0 0 0 0 0 0 00 0 10 0 0 0 0 0 00 0 0 15 0 0 0 0 00 0 0 0 20 0 0 0 00 0 0 0 0 25 0 0 00 0 0 0 0 0 30 0 00 0 0 0 0 0 0 35 00 0 0 0 0 0 0 0 40a =95.生成一个与A 同阶的正态分布的随机矩阵B, 要求B矩阵元素为介于0~99之间整数>> B=round(0+rand(9)*(99-0))B =75 10 79 65 25 83 47 38 4674 82 19 16 50 25 35 56 139 69 48 12 69 81 82 8 3365 31 44 49 88 24 58 5 1617 94 64 95 95 92 54 53 7970 3 70 34 54 35 91 77 313 43 75 58 14 19 28 92 5227 38 27 22 15 25 75 13 165 76 67 74 25 61 75 56 606.输出A 与B 的乘积矩阵C>> C=A*BC =75 10 79 65 25 83 47 38 46222 246 57 48 150 75 105 168 3390 690 480 120 690 810 820 80 330975 465 660 735 1320 360 870 75 240340 1880 1280 1900 1900 1840 1080 1060 15801750 75 1750 850 1350 875 2275 1925 77590 1290 2250 1740 420 570 840 2760 1560945 1330 945 770 525 875 2625 455 560200 3040 2680 2960 1000 2440 3000 2240 24007.生成由A 与B 点乘得到的矩阵D>> D=A.*BD =75 0 0 0 0 0 0 0 00 246 0 0 0 0 0 0 00 0 480 0 0 0 0 0 00 0 0 735 0 0 0 0 00 0 0 0 1900 0 0 0 00 0 0 0 0 875 0 0 00 0 0 0 0 0 840 0 00 0 0 0 0 0 0 455 00 0 0 0 0 0 0 0 24008.生成一个由D 的第8、4、2、5 行和第7、1、6、9、2 列组成的子矩阵 E>> E=D([8 4 2 5],[7 1 6 9 2])E =0 0 0 0 00 0 0 0 00 0 0 0 2460 0 0 0 09.求出矩阵E 的最大元素>> m=max(max(E))m =24610.计算的值>> 1.369^2+sin(7*pi/10)*sqrt(26.48)/2.9ans =3.309711.产生一个5阶魔术方阵，并执行如下操作：a)将矩阵的第2行3列元素赋值给变量cb)将由矩阵第2，3，4行第3，5列构成的子矩阵赋值给变量d>> Q=magic(5)c=Q(2,3)d=[Q(2:4,3),Q(2:4,5)]Q =17 24 1 8 1523 5 7 14 164 6 13 20 2210 12 19 21 311 18 25 2 9c =7d =7 1613 2219 312.给出区间[0，1]上的6个等分点数据。
Matlab上机实验(第四次和第五次)
Matlab 上机实验（第四次和第五次）要求：每次上机要签到，以下题目在第四次上机时以电子版形式交到共享文件夹中，共享文件夹位置见公共邮箱。
上交电子文档的文件名按照专业班级学号姓名顺序给出，比如应用电子0701班李伟文（学号3106002803）上交的文件名为：应电0701班3106002803李伟文。
上机地点：信息学院机房404，406。
1.设,求>> syms x;f=x*exp(2*x);diff(f,x,20)ans = 10485760*exp(2*x)+1048576*x*exp(2*x)2.设,求>> syms x;f=sin(x).^4+cos(x).^4;diff(f,x,10)ans =1572864*sin(x)^2*cos(x)^2-262144*cos(x)^4-262144*sin(x)^43.设,求>> syms x y;z=y^4*sin(x)^4+x^4*cos(y)^4;diff(z,x)ans =4*y^4*sin(x)^3*cos(x)+4*x^3*cos(y)^44．⎰++dx x x x ))(1(122>> syms x;f=((x^2+1)*(x^2+x))\1;int(f,x)ans = log(x)-1/2*log(x+1)-1/4*log(x^2+1)-1/2*atan(x)5．dx x x ⎰-632)9(>> syms x;f=x^6\(9-x^2)^(2\3);int(f,x)ans =-1/45/x^5*(9-x^2)^(5/2)6．⎰+213x x dx >> syms x;f=(x+x^3)\1;int(f,x,1,2)ans = 3/2*log(2)-1/2*log(5)7.求级数的和: ∑∞=+122n n n （利用帮助，查找symsum ）>> syms n;symsum(2^n\(n+2),1,inf)ans =48.利用函数int 计算二重不定积分 >> syms x y;z=(x+y)*exp(-x*y);t=int(int(z,x),y)t =x*Ei(1,x*y)+x*(1/x/y*exp(-x*y)-Ei(1,x*y))+1/x*exp(-x*y)9．设方程（1）0)cos(3=+a x ，（2）0)cos(3=+x x 及（3）0)cos(3=+ax x ，用符号运算工具箱函数solve 分别求x 的解。
Matlab上机作业部分参考答案
参考答案
2. 求方程组的根
解： >> a=[1 4 -3;2 5 -1;1 6 1]; >> b=[2;11;12]; >> x=a\b 2x1+5x2-x3=11 x1+6x2+x3=12
3. 解方程组:
2 9 0 13 3 4 11 x 6 2 2 6 6
[答案] >> a=[2 9 0;3 4 11;2 2 6]; >> b=[13 6 6]'; >> x=a\b x= 7.4000 -0.2000 -1.4000
4. 试判定下面的线性代数方程是否有解。
3 13 16 2 1 5 11 10 8 3 X 9 4 7 6 12 4 14 15 1 7
【求解】可以将方程写成矩阵形式，得出的两列向量为方程 的基础解系。 >> A=[6,1,4,-7,-3; -2,-7,-8,6,0; -4,5,1,-6,8; -34,36,9,-21,49; -26,-12,-27,27,17]; A=sym(A); rank(A) ans = 3 >> null(A) ans = [ 191/34, 95/17] [ 0, 1] [ 1, 0] [ 109/34, 103/34] [ 173/34, 151/34]
A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ]; b=[4 3 -1 8 -5 2]'; B=[A b]; C=rref(B) A=[1 2 4 6 -3 2;2 4 -4 5 1 -5;3 6 2 0 5 -9;2 3 0 4 0 1;0 -4 -5 2 1 4 ;5 5 -3 6 6 -4 ]; b=[4 3 -1 8 -5 2]'; B=[A b]; n=6; R_A=rank(A) R_B=rank(B) format rat if R_A==R_B&R_A==n X=A\b elseif R_A==R_B&R_A<n X=A\b C=null(A,'r') else X='equation has no solves' end syms k1 k2 X=k1*C+X pretty(X)
Matlab上机作业部分参考答案
3. 有一组学生的考试成绩如下表，根据规定，成绩在100分时 为满分，成绩在90~99之间时为优秀，成绩在80~89分之间时 为良好，成绩在60~79分之间为及格，成绩在60分以下时为不 及格，编制一个根据成绩划分等级的程序。
学生 姓名 王 张 刘 李 陈 杨 于 黄 郭 赵成 绩Βιβλιοθήκη 728356
94
100
z1 x2 2y 2 , 2. 列出求下列空间曲面
（提示find 命令）
z2 2 x 3 y 交线的程序
clc clear x=[-5:0.5:5]; [X,Y]=meshgrid(x); z1=X.^2-2*Y.^2; z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z')
surf(X,Y,z1) hold on surf(X,Y,z2) k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k) z3=x1.^2-2*y1.^2 hold on
88
96
68
54
65
4. 编制写一个函数实现，输入一个3x4矩阵，求其中最大元素 和最小元素（禁止直接用Matlab命令） 5. 编程，可输入有10个实数的数组，调用函数f求他们的平均值 V，然后输出其中与V之差的绝对值为最小的那个数组元素。 6. 编制写一个函数实现，一个给定二维（3x3）数组转置。 （禁止直接用Matlab命令）
7.
编制写一个函数实现，输入16进制，输出10进制。（禁止直 接用Matlab命令）
8. 编程序实现功能：在一个元素值按升序存放的整型数组中插 入一个数，使得插入后的数组元素仍然有序。 9. 编制写一个函数实现，用“冒泡法” 排列输入的10个数字。
MATLAB上机题.doc
西北农林科技大学本科课程试卷2011—2012学年第1学期《MATLAB与科学计算》课程A卷（上机考试）专业班级：09农水班命题教师：王增红审题教师：马孝义学生姓名学号：考试成绩：考试须知：1、所有考生请在规定时间内，完成以下操作内容，并将.m文件和操作结果粘贴到试卷中。
2、考试过程中，请注意保存文件，因个人操作造成的文件丢失，按0分计。
3、编程题请注意编程规范，否则酌情扣分。
4、请将姓名和学号填入试卷中指定位置，word文件名请按“s机号-学号-姓名”规则保存。
5、考试结束后停止操作，学生和教师确认文件拷贝无误后，方可离开考场，否则后果自负。
y2=[2 5 4;9 10 5]plot(x1,y1)%x1的列与y1的行相等，与y1的每列对应一条直线，两条直线plot(x1,x2)%x1的列与x2的列相等，与x2的每行对应一条直线，两条直线plot(x1,y2)%x1的列与y2的列相等，与y2的每行对应一条直线，两条直线plot(y1,x1)%y1的行与x1的列相等，与x1的每行对应一条直线，两条直线plot(x2,x1)%x2的列与x1的列相等，与x1的每行对应一条直线，两条直线plot(y2,x1)%y2的列与x1的列相等，与x1的每行对应一条直线，两条直线plot(x2,y2)%x2与y2的行列相等，与y2的每列对应一条直线，三条直线plot(y2,x2)%y2与x2的行列相等，与x2的每列对应一条直线，三条直线其图形如下所示：二、绘图题（根据题意，对下列各题进行绘图，共20分）得分：分（本大题中的3个小题编辑M文件，运行查看图形，并将M文件粘贴到试卷中）1、画一二维饼形图，至少包含三个元素。
（5分）答：x=[100 200 300 400]explode=[1 1 1 1]pie(x,explode,{'一','二','三','四'})2、在同一窗口绘制6个子图，已知t=::2*pi，每个子图分别绘制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t,y5=2ln(t),y6=e t）（10分）答：t=::2*pisubplot(2,3,1)y1=sin(t)plot(t,y1)title('y1=sin(t)')subplot(2,3,2)y2=cos(2*t)plot(t,y2)title('y2=cos(2t)')subplot(2,3,3)y3=sin(t)+cos(t)plot(t,y3)title('y3=sin(t)+cos(t)')subplot(2,3,4)y4=3*tplot(t,y4)title('y4=3t') subplot(2,3,5) y5=2*log(t) plot(t,y5)title('y5=2ln(t)') subplot(2,3,6) y6=exp(t) plot(t,y6) title('y6=e^t')3、在matlab 的图形窗口中写出标题为表达式⎰⎰±=m mzdx ydx qt 00，字体大小为14号，如右图所示。
Matlab上机作业部分参考答案
上机练习二 参考答案
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并 且按照从大到小的顺序排列好！ 【求解】 a=10*rand(1,10)-5; b=sort(a,'descend')
上机练习二 参考答案
2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵A将得出什么 结果？
Matlab 上机课作业
吴梅红 2012.10.15
上机练习一
上机练习一 参考答案
上机练习一 参考答案
上机练习一 参考答案
上机练习二
1. 产生一个1x10的随机矩阵，大小位于（-5 5），并且按 照从大到小的顺序排列好！ 2、用MATLAB 语句输入矩阵A 和B
前面给出的是4 ×4 矩阵，如果给出A(5，6) = 5 命令，矩阵 A将得出什么结果？ 3、假设已知矩阵A ，试给出相应的MATLAB 命令，将其全 部偶数行提取出来，赋给B 矩阵，用A =magic(8) 命令生成A 矩阵，用上述的命令检验一下结果是不是正确。
【求解】用课程介绍的方法可以直接输入这两个矩阵 >> A=[1 2 3 4; 4 3 2 1; 2 3 4 1; 3 2 4 1] A= 1234 4321 2341 3241 若给出A(5,6)=5 命令，虽然这时的行和列数均大于A矩阵当前的维数， 但仍然可以执行该语句，得出 >> A(5,6)=5 A= 123400 432100 234100 324100 000005 复数矩阵也可以用直观的语句输入 3+2i 4+1i; 4+1i 3+2i 2+3i 1+4i; 2+3i 3+2i 4+1i 1+4i; 3+2i 2+3i 4+1i 1+4i]; B= 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i
matlab上机习题详解 试题答案
第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
同济大学 Matlab 第五次实验报告
-4
10
-3
10
-2
Inverse Tenperature K-1
Inverse Tenperature K-1
1 γ 铁的扩散率与温度倒数 T 的曲线图：
1.4 1.2
2 Diffusivity cm /s 2 Diffusivity cm /s
x 10
-6
Diffusivity of C in Fe 线线 1.4 1.2 1 0.8 0.6 0.4 0.2 1 1.5 2 2.5 3 3.5 x 10
1 α 铁的扩散率与温度倒数 T 的曲线图：
对 数 '),xlabel('Inverse Tenperature
对 数 '),xlabel('Inverse Tenperature
数 '),xlabel('Inverse
Tenperature
性 '),xlabel('Inverse
Tenperature
实 验 报 告
班级 学号 姓名 Matlab 第五次试验 1． 自学 SIMULINK 的封装和自定义功能模块的建立 2． 对于本领域常用工具学会建立自定义功能模块以 供使用方便 实 验 目 的 3． 将建立过程的重要步骤截图表示即可 4． 半对数坐标图的绘制
5． ch 结构的应用，结构数组的应用
-3
x 10
-6
Diffusivity of C in Fe x轴轴轴轴
1 0.8 0.6 0.4 0.2 0 0.5
Inverse Tenperature K-1 Diffusivity of C in Fe x轴轴轴轴
0 -4 10
10
Matlab上机练习参考答案
Matlab 上机练习二班级 学号 姓名按要求完成题目，并写下指令和运行结果。
（不需要画图）1、 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ans =- -- ++ -- ++ -2、计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8];>> a.*bans =12 36 38 42 403、 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;>> X=A\BX =4、 ⎥⎦⎤⎢⎣⎡-=463521a ，⎥⎦⎤⎢⎣⎡-=263478b ，观察a 与b 之间的六种关系运算的结果。
>> a=[1 2 3;4 5 6];>> b=[8 –7 4;3 6 2];>> a>bans =0 1 01 0 1 >> a>=bans =0 1 01 0 1 >> a<bans =1 0 1 0 1 0 >> a<=bans =1 0 1 0 1 0 >> a==bans =0 0 0 0 0 0>> a~=bans =1 1 11 1 15、[]7.0=-a，在进行逻辑运算时，a相当于什么样的逻辑量。
82.05-相当于a=[1 1 0 1 1]。
通信工程MATLAB上机实验12345
MATLAB实验指导书王海燕石油大学计通学院通信工程系2008年9月实验一熟悉MATLAB环境[实验目的]1．熟悉MATLAB主界面，学会常用的窗口的操作2．学会简单的矩阵输入3．熟悉MATLAB基本语句的使用、矩阵与向量的基本运算[实验原理]1. 用户第一次使用MATLAB时，可通过demo命令，启动MATLAB的演试程序，领略MATLAB所提供的强大的运算与绘图功能。
也可以键入help进行进一步了解。
2.MATLAB 是以矩阵作为基本编程单元的一种程序设计语言。
它提供了各种矩阵的运算与操作，并有较强的绘图功能。
在MATLAB中输入矩阵最方便的的方式是直接输入矩阵的元素，其方法为：1）用中括号[ ]把所有矩阵元素括起来；2）同一行的不同数据元素之间用空格或逗号间隔；3）用分号（；）指定一行结束；4）也可以分成几行进行输入，用回车符代替分号；5）矩阵元素可以是任何matlab表达式（系统将自动计算结果），可以是实数，也可以是复数，复数可用特殊函数i，j 输入。
3.MATLAB语言最基本的赋值语句结构为：变量=表达式。
表达式由操作符或其它字符，函数和变量名组成，表达式的结果为一个矩阵，显示在屏幕上，同时输送到一个变量(ans为默认)中并存放于workspace工作空间中以备调用。
[实验内容]1．（1）用空格或逗号输入矩阵A=[1 2 3；4 5 6；7 8 9]或A=[1,2,3;4,5,6;7,8,9]（2）试用回车代替A中的分号，观察输出结果。
要求：总结矩阵的输入2．输入矩阵B=[9，8，7；6，5，4；3，2，1]；C=[4，5，6；7，8，9；1，2，3]；要求：（1）说明；的作用（2）给出查看B及C的结果的方法？3．用edit或选择File|new菜单中的M-file进入编辑器，输入:Base=[9,8,7;6,5,4;3,2,1];Base_info=Base+1保存为Base.m文件，并通过Debug中的run命令（或F5）运行之，观察有何异常现象。
现代控制理论MATLAB上机实验指导1
《现代控制理论》MATLAB实践指导书刘红军编写《现代控制理论》MATLAB实践指导书1 MA TLAB概述MA TLAB是MA Trix LABoratory的缩写，早期主要用于现代控制中复杂的矩阵、向量的各种运算。
由于MA TLAB提供了强大的矩阵处理和绘图功能，很多专家因此在自己擅长的领域用它编写了许多专门的MA TLAB工具包（toolbox），如控制系统工具包（control systems toolbox）；系统辨识工具包（system identification toolbox）；信号处理工具包（signal processing toolbox）；鲁棒控制工具包（robust control toolbox）；最优化工具包(optimization toolbox）等等。
由于MA TLAB功能的不断扩展，所以现在的MA TLAB已不仅仅局限与现代控制系统分析和综合应用，它已是一种包罗众多学科的功能强大的“技术计算语言（The Language of Technical Computing）”。
MathWorks 公司于1992 年推出了具有划时代意义的MA TLAB 4.0 版本，并推出了交互式模型输入与仿真系统SIMULINK，它使得控制系统的仿真与CAD应用更加方便、快捷，用户可以方便地在计算机上建模和仿真实验。
1997 年MathWorks推出的MA TLAB 5.0 版允许了更多的数据结构，1999 年初推出的MA TLAB 5.3 版在很多方面又进一步改进了MA TLAB 语言的功能。
2000 年底推出的MA TLAB 6.0。
最新版本是MA TLAB7.0。
MA TLAB以矩阵作为基本编程单元，它提供了各种矩阵的运算与操作，并有较强的绘图功能。
MA TLAB集科学计算、图像处理、声音处理于一身，是一个高度的集成系统，有良好的用户界面，并有良好的帮助功能。
MA TLAB不仅流行于控制界，在机械工程、生物工程、语音处理、图像处理、信号分析、计算机技术等各行各业中都有极广泛的应用。
Matlab上机练习一答案
Matlab上机练习⼀答案Matlab 上机练习⼀班级学号姓名按要求完成题⽬，并写下指令和运⾏结果。
（不需要画图）1、求下列联⽴⽅程的解=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59942、设 ??++=)1(sin 35.0cos 2x x x y ，把x=0-2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')3、产⽣8×6阶的正态分布随机数矩阵R1, 求其各列的平均值、和、乘积。
并求该矩阵全体数的平均值。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)或者u=reshape(a,1,48);p1=mean(u)p2=var(u)4、绘制)(222y x e x z +-=在定义域x=[-2,2],y=[-2,2]内的曲⾯。
（利⽤meshgrid ）x=-2:2;y=x;[X,Y]= meshgrid(x,y);Z=X^2*exp(-(X^2+Y^2));mesh(X,Y,Z)5、求代数⽅程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。
（利⽤roots 函数）p=[3 4 7 2 9 12];roots(p)6、把1开五次⽅，并求其全部五个根。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第五次上机作业时间：2013-6-25 E-207（周二）晚上6：00~10：00班级：071111学号：07111014姓名：曹红兴符号运算：1. 求符号矩阵111213212223313233a a a A a a a a a a ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦的行列式值和逆，所得结果应采用“子表达式置换”简洁化。
程序：A = sym ( ' [a11 a12 a13 ; a21 a22 a23 ; a31 a32 a33] ' );B = det(A);C = inv(A);[ Bs,e] = subexpr( B,'e')[ Cs,d] = subexpr( C,'d' )运行结果：Bs =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 e =[ empty sym ]Cs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)][ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)][ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)]d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) >>2. 对函数0()00ka k f k k ⎧≥=⎨<⎩，当a 为正实数，求0()k k f k z ∞-=∑。
程序：syms k a zf = a^k * z^(-k);S = symsum(f,k,0,inf)运行结果：S =piecewise([a = z, Inf], [a <> 0 and z = 0, 0], [a <> z and z <> 0, limit((a^k*z)/(z^k*(a - z)), k = Inf) - z/(a - z)])>>3. 对0x >，求21021211k k x k x +∞=-⎛⎫ ⎪++⎝⎭∑。
程序：syms k xf = 2*((x-1)/(x+1))^(2*k+1)/(2*k+1)S = symsum(f,k,0,inf)运行结果：f =(2*((x - 1)/(x + 1))^(2*k + 1))/(2*k + 1)S =piecewise([abs(x + 1)*(abs(x - 1) - abs(x + 1)) < 0, 2*atanh((x - 1)/(x + 1))])>>4. 求()|sin |y t t =的导数dy dt ，0t dy dt -=，2t dy dt π= 程序：syms xf = abs(sin(x))g = diff(f,x)g1 = subs(g,x,-eps)g2 = subs(g,x,pi/2) 运行结果：f =abs(sin(x))g =sign(sin(x))*cos(x)g1 =-1g2 =6.1232e-017>>5. 计算二重积分222211()x x y dydx +⎰⎰程序：syms x y f = int(int((x^2+y^2),y,1,x^2),x,1,2)vpa(f)运行结果：f =1006/105ans =9.5809523809523809523809523809524>>6. 求||1||()0||t A t f t t τττ⎧⎛⎫-≤⎪ ⎪=⎝⎭⎨⎪>⎩的Fourier 变换，并画出2,2A τ==时的频谱程序；syms A t w tao ;ft = A*(1-abs(t)/tao)*(heaviside(t+tao)-heaviside(t-tao));Fw = fourier(ft,t,w)h = subs(Fw,A,2);h = subs(h,tao,2);x = linspace(-20,20,1000)；sn = subs(h,w,x);plot(x,sn)运行结果：Fw =A*((sign(tao)/w^2 + 1/w^2 + pi*dirac(-w, 1)*i - (exp(tao*w*i)*sign(tao))/w^2 + (tao*exp(tao*w*i)*sign(tao)*i)/w)/tao - (1/w^2 - sign(tao)/w^2 + pi*dirac(-w, 1)*i + ((1/exp(tao*w*i))*sign(tao))/w^2 + (tao*(1/exp(tao*w*i))*sign(tao)*i)/w)/tao + (1/exp(tao*w*i))*(- pi*dirac(-w) + i/w) - exp(tao*w*i)*(- pi*dirac(-w) + i/w))Warning: Imaginary parts of complex X and/or Y arguments ignored> In workfive at 10>>8．用solve 求方程组2212x y xy ⎧+=⎨=⎩的解。
程序：syms x yS = solve(x^2+y^2-1,x*y-2,x,y)x1 = S.xx2 = S.y运行结果：S =x: [4x1 sym]y: [4x1 sym]x1 =(1/2 + (15^(1/2)*i)/2)^(1/2)/2 - (1/2 + (15^(1/2)*i)/2)^(3/2)/2 - (1/2 + (15^(1/2)*i)/2)^(1/2)/2 + (1/2 + (15^(1/2)*i)/2)^(3/2)/2 (1/2 - (15^(1/2)*i)/2)^(1/2)/2 - (1/2 - (15^(1/2)*i)/2)^(3/2)/2 - (1/2 - (15^(1/2)*i)/2)^(1/2)/2 + (1/2 - (15^(1/2)*i)/2)^(3/2)/2 x2 =(1/2 + (15^(1/2)*i)/2)^(1/2)-((15^(1/2)*i)/2 + 1/2)^(1/2)(1/2 - (15^(1/2)*i)/2)^(1/2)-(1/2 - (15^(1/2)*i)/2)^(1/2)>>9. 求一阶微分方程2',(0)2x at bt x =+=的解程序：syms a b t xS = dsolve('Dx = a*t^2+b*t','x(0) = 2','t')运行程序：S =(t^2*(3*b + 2*a*t))/6 + 2>>10. 求边值问题34,43,(0)0,(0)1df dg f g f g f g dx dx=+=-+==的解，并绘制解的图形程序： S = dsolve('Df = 3*f + 4*g ','Dg = -4*f+3*g','f(0) = 0','g(0) = 1' ,'x') f1 = S.fg1 = S.gclf;hold onf11 = ezplot(f1,[20,60]);set(f11,'Color','k')g11 = ezplot(g1,[20,60]);set(g11,'Color','r')grid on运行结果：S =f: [1x1 sym]g: [1x1 sym]f1 =sin(4*x)*exp(3*x)g1 =cos(4*x)*exp(3*x)>>●“证明”“哥德巴赫猜想”：1742年6月7日哥德巴赫写信给当时的大数学家欧拉，正式提出了以下的猜想：a.任何一个大于6的偶数都可以表示成两个素数之和。
b.任何一个大于9的奇数都可以表示成三个素数之和。
陈景润于1966年证明：任何充份大的偶数都是一个质数与一个自然数之和，而后者可表示为两个质数的乘积。
任务：用计算机证明106以内的数满足歌德巴赫猜想程序：a(1) = 2;a(2) = 3;k = 3;for i = 1:10^6%求出10^6以内的所有素数，记录在数组a中for j=2:floor(sqrt(i))if mod(i,j) == 0break;endif j == floor(sqrt(i))a(k) = i;k = k + 1;endendendflag = 0; %标记对于整数n，是否找到要求的俩个或三个素数。
如果找到，赋值为1；否则赋值为0for n = 1:10^6if mod(n,2)==0 && n > 6%大于6的偶数for i = 1:length(a)if flag==1break;endfor j = i:length(a)if n==a(i) + a(j) %找到了这个偶数所对应的俩个素数flag = 1;break;endendendif flag == 0disp('error!'); %如果没有找到要求的三个素数，输出：error！endendflag = 0;if mod(n,2)==1&& n > 9%大于9的奇数for i = 1:length(a)if flag == 1break;endfor j = i:length(a)if flag == 1break;endfor k = j:length(a)if n==a(i) + a(j) +a(k) %找到了这个奇数所对应的三个素数flag = 1;break;endendendendif flag == 0disp('error!') %如果没有找到要求的三个素数，输出：error！endendflag = 0;end注：10^6运行时间太久，只验证了10^4，运行结果正确。