matlab第四次作业
物流管理定量分析方法形成性考核册(第3版)答案
第一次作业物资调运方案优化的表上作业法1.若某物资的总供应量大于总需求量，则可增设一个（A )，其需求量取总供应量与总需求量的差额，并取各产地到该销地的单位运价为0，可将不平衡运输问题化为平衡运输问题。
（A）虚销地（B)虚产地（C）需求量 (D）供应量2。
将下列某物资的供求不平衡运输问题（供应量、需求量单位:吨;运价单位：元/吨)化为供求平衡运输问题：量与总供应量的差额，并取该产地到各销地的单位运价为0，并将供不应求运输问题化为供求平衡运输问题。
(A）大于 (B）小于 (C）等于（D)大于等于4．将下列某物资的供求不平衡运输问题（供应量、需求量单位:吨；运价单位：元/吨)化为供求平衡运输问题：5A,B，C，D四个仓库中收存，四仓库收进的数量分别为100吨、1500吨、400吨和1100吨，仓库和发货点之间的单位运价如下表所示：/吨）试用最小元素法确定一个初始调运方案，再调整寻求最优调运方案，使运输总费用最小。
解：6。
某物资要从产地A1,A2,A3调往销地B1，B2，B3，运输平衡表（单位：吨）和运价表（单位：元/吨）如下表所示：试用最小元素法编制初始调运方案，并求最优调运方案。
所有检验数全为正，此调运方案最优.最低运输总费用：（元)7。
设某物资要从产地A1，A2,A3调往销地B1，B2,B3，B4，运输平衡表(单位：吨）和运价表(单位：百元/吨）如下表所示：试问应怎样调运才能使总运费最省？所有检验数全为正,初始调运方案就是最优调运方案.最小运输总费用为（元）8．有一运输问题，涉及3个起始点A1，A2，A3和4个目的点B1，B2，B3，B4，3个起始点的供应量分别为50吨、50吨、75吨，4个目的点的需求量分别为40吨、55吨、60吨、20吨.运输平衡表及各起始点与目的点之间的距离(公里）如下表所示：假设每次装车的额外费用不计，运输成本与所行驶的距离成正比。
试求最优的调运方案，并求最小吨公里数。
检验数全为正，达到最优调运方案。
matlab第4次作业
《数学实验》报告实验名称线性代数相关运算及数值方法计算定积分学院专业班级姓名学号2013年 5月一、【实验目的】1、掌握矩阵的基本运算、特征值、特征向量和线性方程组的求解；2、能熟练运用数值方法求定积分。
二、【实验任务】P114--12、14;P115--21 (1) (2);P167--17(2)、18。
三、【实验程序】P114—12：>> A=rand(6,6)>> A'>> det(A)>> rank(A)>> rref(A)P114—14：>> p=poly(A)>> poly2str(p,'x')>> [V,D]=eig(A)P115--21 (1)：>> A=[1 1 2 -1;-1 1 3 0;2 -3 4 -1]>> rref(A)P115--21 (2)：>> B=[1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -1/2]>> rref(B)P167--17(2)：>> x=0:pi/100:pi;>> y=x.*sin(x)./(1.+(cos(x)).^2);>> s1=sum(y(1:100))*pi/100s2=sum(y(2:101))*pi/100s3=trapz(x,y)s4=quad('x.*sin(x)./(1.+(cos(x)).^2)',0,pi)P167—18：>> h=pi/50;x=0:h:pi/4;>> y=1./(1-sin(x));>> s1=sum(y(1:length(x)-1))*h;>> s2=sum(y(2:length(x)))*h;>> s3=trapz(x,y);>> ff=inline('1./(1-sin(x))','x');>> s4=quad(ff,0,pi/4)a1=s1-2^0.5a2=s2-2^0.5a3=s3-2^0.5a4=s4-2^0.5四、【实验结果】P114—12：P114—14：P115--21 (1)：结果分析：系数矩阵A的秩为3，小于未知量个数4，有无穷多解，原方程组对应的通解方程组为：X1-0.56X4=0X2-0.2X4=0X3-0.12X4=0取X4=1解得方程组的基础解系为：0.560.2ξ1=0.121所以方程组的通解为：X1 0.56X2 0.2=KX3 0.12，其中K为任意实数。
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。
矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。
4 作业-MatLab练习题-绘图
MatLab 练习题-绘图学号：班级：姓名：注意: 1、请按要求作题（即写出程序代码+结果）.2、整理好格式，排版好，不留过多的空格。
1,在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注“y=sin(x)”, “y=cos(x)” ,x轴,y轴,标题为“正弦余弦函数图象”.答案：x=linspace(0,4*pi,40);y=sin(x);z=cos(x);plot(x,y,x,z,'ro')xlabel('x')ylabel('y')title('ÕýÏÒÓàÏÒº¯ÊýͼÏñ')gtext('sinx'),gtext('cosx')xy2, 画出椭圆221 169x y+=。
答案：a=linspace(0,2*pi,40); x=cos(a)*4; y=sin(a)*2; plot(x,y); xlabel('x') ylabel('y') title('ÍÖÔ²')-4-3-2-101234-2-1.5-1-0.500.511.52xy椭圆3,绘制螺旋线. (12cos )(12sin )3x u y u z u =+⎧⎪=+⎨⎪=⎩[0,10]u ∈答案：u=0:0.1:10;plot3(1+2*cos(u),1+2*sin(u),3*u) xlabel('x') ylabel('y') title('ÂÝÐýÏß')x螺旋线y4,绘制四叶玫瑰线sin 2ρθ=。
matlab1-8章课后作业
MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。
Matlab第二章1．创建double的变量，并进行计算。
(1)a=87，b=190，计算a+b、a-b、a*b。
(2)创建uint8 类型的变量，数值与(1)中相同，进行相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u)ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。
数学建模作业题+答案
数学建模MATLAB 语言及应用上机作业11. 在matlab 中建立一个矩阵135792468101234501234A ⎡⎤⎢⎥⎢⎥=⎢⎥-----⎢⎥⎣⎦答案：A = [1,3,5,7,9;2,4,6,8,10;-1,-2,-3,-4,-5;0,1,2,3,4]2. 试着利用matlab 求解出下列方程的解（线性代数22页例14）123412423412342583692254760x x x x x x x x x x x x x x +-+=⎧⎪--=⎪⎨-+=-⎪⎪+-+=⎩ 答案：A=[2 ,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6]; B=[8;9;-5;0]; X=A\B 或A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6] b=[8,9,-5,0]' X=inv(A)*b3. 生成一个5阶服从标准正态分布的随机方阵，并计算出其行列式的值，逆矩阵以及转置矩阵。
答案:A=randn(5) det(A) inv(A) A'4. 利用matlab 求解出110430002A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦的特征值和特征向量。
答案：A=[-1,1,0;-4,3,0;0,0,2] [V,D]=eig(A)5.画出衰减振荡曲线3sin3t y et -=在[0,4]π上的图像。
要求，画线颜色调整为黑色，画布底面为白色。
（在实际中，很多打印机时黑白的，因此大多数作图要考虑黑白打印机的效果。
） 给出恰当的x ，y 坐标轴标题，图像x 轴的最大值为4π。
6. 生成一个0-1分布的具有10个元素的随机向量，试着编写程序挑选出向量中大于0.5的元素。
数学建模和Matlab 上机作业2（2016-9-20）跟老师做（不用整合进作业中）：上机演示讲解：函数，递归的两个例子的写法。
附：1. Fibonacci Sequence （斐波那契数列）在数学上，费波那西数列是以递归的方法来定义： F1= 1；F2= 1；F （n ）=F （n-1）+F （n-2） 2. 阶乘举例：数学描述：n!=1×2×……×n ；计算机描述：n!=n*(n-1)!自己做（需要整合进作业中，提交到系统中）：1. 写一个m 文件完成分值百分制到5分制的转换（即输入一个百分制，转换后输出一个5级对应的得分，联系条件控制语句）。
MATLAB第四次作业
《工程应用数学》第四次作业化学工程张少石622085216041 一、目的熟悉迭代法求解线性方程组的过程，通过MATLAB的m文件详细了解各种不同迭代法构造的迭代函数的差别进而巩固迭代求解线性方程组的思想。
尝试使用Jacobi迭代法计算线性方程组。
二、思路Jacobi文件如下：function x = Jacobi(A, c, x0, tol, trace)if nargin < 4 | isempty(tol)tol = 1e-6;endif nargin >= 4 & tol == 0tol = 1e-6;endif nargin < 5 | isempty(trace)trace = 0;endif tracefprintf('\n Initial guess :\n')fprintf('%8.6g ',x0)endc = (c(:).')'; % Make sure it's a column vectorx0 = (x0(:).')'; % Make sure it's a column vectorn = length(c);[nr nc] = size(A);% Check coefficient matrix, vector of constants and% vector of unknownsif nr ~= ncerror('Coefficient matrix is not square.')endif nr ~= nerror('Coefficient matrix and vector of constants do not have the same length.')endif length(x0) ~= nerror('Vector of unknowns and vector of constants do not have the same length.')end% Check if the coefficient matrix is singularif det(A) == 0fprintf('\n Rank = %7.3g\n',rank(A))error('The coefficient matrix is singular.')end% Building modified coefficient matrix and modified% vector of coefficientsD = diag(diag(A)); % The diagonal matrixa0 = inv(D)*A - eye(n); % Modified matrix of coefficientsc0 = inv(D)*c; % Modified vector of constantsx = x0;x0 = x + 2 * tol;iter = 0;% Substitution procedurewhile max(abs(x - x0)) >= tolx0 = x;x = c0 - a0 * x0;if traceiter = iter + 1;fprintf('\n Iteration no. %3d\n',iter)fprintf('%8.6g ',x)endend三、计算过程要计算的线性方程组如下：4x1+3x3=243x1+3x2−x3=30−x2+4x3=−24在MATLAB命令窗口运行如下：>> A=[4 3 0;3 3 -1;0 -1 4];>> b=[24;30;-24];>> x0=[1 1 1];>> x=Jacobi(A,b,x0)x =-3.000012.0000-3.0000四、结果与讨论虽然迭代求解线性方程组只能得到近似解，但在提高精度要求后可以得到十分可靠地近似解。
matlab1-8章课后作业
MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。
Matlab第二章1．创建double的变量，并进行计算。
(1)a=87，b=190，计算a+b、a-b、a*b。
(2)创建uint8 类型的变量，数值与(1)中相同，进行相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =ans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans =15.4189>> sqrt(u-3*v)/(u*v) ans =0 + 0.4410i 4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。
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第四次模板
数学软件实验报告
实验序号：04 实验项目名称：MATLAB在线性代数中的应用
7、运行在MATLAB中迭代求解线性代数方程组的函数文件
二、实验设备（环境）及要求
针对所给习题，编写程序。实验过程中，务必分析实验结果，按要求写
出实验报告。（建议同时网上提交电子版实验报告）
三、实验内容与步骤
1.输入矩阵，求得矩阵的维数
2.求矩阵的范数、向量的范数
3.求矩阵的秩
4.求矩阵的条件数
5.求矩阵的行列式
>> rank(a),rank(b)
ans =
2
ans =
2
4. 矩阵的条件数
>> a=[1 2 3 4];b=3*ones(5);c=[1:5;3:7;5:9;(10*rand(5));
>> e=diag(a,1);
>> f=e+b+c+d;
>> cond(f,1),cond(f),cond(f,inf)
6.求非奇异阵的逆矩阵
7.利用‘\’‘/’解线性代数方程组
8.求方阵的特征值与特征向量
9.求矩阵的奇异值分解
10.求矩阵的QR分解
11.求矩阵的LU分解
12.求矩阵的Jordan分解
13.用迭代法求解线性代数方程组和利用LU分解法求解线性代数方程组
四、实验程序、结果与数据处理
1. 输入矩阵，并求得矩阵的维数
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第四讲上机实验
实验内容（注：1，2无需写入报告）
1．用命令或图标进入simulink模块库，仔细观察各个模块组并分别打开，熟悉其模块的构成。
回忆连续模块组中有哪些常用模块？Sum模块在哪个模块组中？
2．用3种不同的方法，分别打开一空白模型窗口，熟悉有关的菜单及选项，并向其中拖拽一模块，如Step模块，对其进行大小的调整、方向的旋转、标签名字的改变及隐藏、增加阴影、参数改变等操作，熟悉右键的使用。
3．在Simulink中构建模型：对一个正弦信号进行积分处理，然后将原始正弦信号和积分后的信号送到示波器中同时显示出来。
4． 一系统可由微分方程322()y y y u t ++=&
&&来描述，系统初始状态为零，求取系统响应。
5. 建立Simulink模型，用信号发生器产生幅值4.5V、频率为1Hz的正弦波，并叠加一个均值0.1V, 方差0.5V的白噪声（正态分布）信号，将叠加后的信号显示在示波器上, 并同时传送至matlab工作空间。