插值-matlab实验课件

插值-matlab实验课件
插值-matlab实验课件

Matlab数学实验报告一

数学软件课程设计 题目非线性方程求解 班级数学081 姓名曹曼伦

实验目的:用二分法与Newton迭代法求解非线性方程的根; 用Matlab函数solve、fzero、fsolve求解非线性方程(组)的解。 编程实现二分法及Newton迭代法; 学会使用Matlab函数solve、fzero、fsolve求解非线性方程(组)的解。 通过实例分别用二分法及迭代法解非线性方程组并观察收敛速度。 实验内容: 比较求exp(x)+10*x-2的根的计算量。(要求误差不超过十的五次方) (1)在区间(0,1)内用二分法; (2)用迭代法x=(2-exp(x))/10,取初值x=0 。 试验程序 (1)二分法: format long syms x s=exp(x)+10*x-2 a=0; b=1; A=subs(s,a) B=subs(s,b) f=A*B %若f<0,则为由根区间 n=0; stop=1.0e-5; while f<0&abs(a-b)>=stop&n<=100; Xk=(a+b)/2; %二分 M= subs(s, Xk); if M* A<0 symbol=1 %若M= subs(s, Xk)为正,则与a二分 b= Xk else symbol=0 % 若M= subs(s, Xk)为负,则与b二分 a= Xk end n=n+1 end Xk n (2)牛顿迭代法; format long

syms x s= (2-exp(x))/10; %迭代公式 f=diff(s); x=0; %迭代初值 a=subs(f,x); %判断收敛性(a是否小于1) s=(2-exp(x))/10; stop=1.0e-5; %迭代的精度 n=0; while a<1&abs(s-x)>=stop&n<=100; x=s %迭代 s=(2-exp(x))/10; n=n+1 end 实验结果: (1)二分法: symbol =1 b =0.50000000000000 n =1 symbol =1 b =0.25000000000000 n =2 symbol =1 b =0.12500000000000 n =3 symbol =0 a =0.06250000000000 n =4 symbol =1 b =0.09375000000000 n =5 symbol =0 a =0.07812500000000 n =6 symbol =1 b =0.09054565429688 n =15 symbol =1 b =0.09053039550781 n =16 symbol =0 a =0.09052276611328 n =17 Xk =0.09052276611328 n =17 (2)迭代法 由x =0.10000000000000 n =1 x =0.08948290819244 n =2 x =0.09063913585958 n =3 x =0.09051261667437 n =4 x =0.09052646805264 n =5 试验结果可见用二分法需要算17次,而用迭代法求得同样精度的解仅用5次,但由于迭代法一般只具有局部收敛性,因此通常不用二分法来求得非线性方程的精确解,而只用它求得根的一个近似解,再用收敛速度较快的迭代法求得其精确解。

数学软件MATLAB实验作业

数学软件与数学实验作业 一.《数学软件》练习题(任选12题,其中19-24题至少选2题): 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)

MATLAB练习题七

实验六 观看课件第三章,学习数学模型的连接(课件80页始) 注意:了解最后结果形式优先级,而传函优先级最低 1、串联: (1)两个系统、环节串联连接时,s3=series(s1,s2),s3结果的形式由优先原则确定。 (2)多个(包含2个)系统、环节串联连接时,可用s6=s2*s3*s4*s5*-- 2、并联: (1)两个系统、环节并联连接时,s3=paralell(s1,s2),s3结果的形式由优先原则确定。 (2)多个(包含2个)系统、环节并联连接时,可用s6=s2+s3+s4+s5+--- 3、反馈:sys=feedback (s1,s2) [ num,den]=cloop(num,den,sign) 习题一、求串联环节的传递函数 G1=tf([1],[1 1]); G2=tf([1], [1 0]); G3=series(G1,G2); G4=tf([1],[1 2]); Y=series(G3,G4) Transfer function: R(s)

1 ----------------- s^3 + 3 s^2 + 2 s 习题二、求并联环节的传递函数 G1=tf([1],[1 1]); G2=tf([1],[1]); Y=parallel(G1,G2) Transfer function: s + 2 ----- s + 1 习题三:已知 6.01,141)(,2314)(,212)(43221=+=+++=+=G s s G s s s s G s s s G ,其连接结构如下图所示,编程求取系统闭环传函。 G2=tf([4 1],[1 3 2]); G3=tf([1],[4 1]);

MATLAB数学实验第二版答案(胡良剑)

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值,b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明:编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明:编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1,但数据类型分别为数值,字符,逻辑,注意a与c相等,但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12; >> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点 f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500 >> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点 >> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500

MATLAB实验报告

实验一 MATLAB 环境的熟悉与基本运算 一、实验目的及要求 1.熟悉MATLAB 的开发环境; 2.掌握MATLAB 的一些常用命令; 3.掌握矩阵、变量、表达式的输入方法及各种基本运算。 二、实验内容 1.熟悉MATLAB 的开发环境: ① MATLAB 的各种窗口: 命令窗口、命令历史窗口、工作空间窗口、当前路径窗口。 ②路径的设置: 建立自己的文件夹,加入到MATLAB 路径中,并保存。 设置当前路径,以方便文件管理。 2.学习使用clc 、clear ,了解其功能和作用。 3.矩阵运算: 已知:A=[1 2;3 4]; B=[5 5;7 8]; 求:A*B 、A.*B ,并比较结果。 4.使用冒号选出指定元素: 已知:A=[1 2 3;4 5 6;7 8 9]; 求:A 中第3列前2个元素;A 中所有列第2,3行的元素; 5.在MATLAB 的命令窗口计算: 1) )2sin(π 2) 5.4)4.05589(÷?+ 6.关系及逻辑运算 1)已知:a=[5:1:15]; b=[1 2 8 8 7 10 12 11 13 14 15],求: y=a==b ,并分析结果 2)已知:X=[0 1;1 0]; Y=[0 0;1 0],求: x&y+x>y ,并分析结果 7.文件操作 1)将0到1000的所有整数,写入到D 盘下的文件 2)读入D 盘下的文件,并赋给变量num

8.符号运算 1)对表达式f=x 3 -1 进行因式分解 2)对表达式f=(2x 2*(x+3)-10)*t ,分别将自变量x 和t 的同类项合并 3)求 3(1)x dz z +? 三、实验报告要求 完成实验内容的3、4、5、6、7、8,写出相应的程序、结果

MATLAB软件与基础数学实验

软件与基础数学实验 实验1 基本特性与基本运算 例1-1求[12+2×(7-4)]÷32的算术运算结果。 >> >> (12+2*(7-4))/3^2 s = 2 例1-2计算5!,并把运算结果赋给变量y 1; 1:5 *i; y 例1-3计算2开平方 >> 2^(0.5) s = 1.4142 >> 例1-4 计算2开平方并赋值给变量x(不显示) 查看x的赋值情况 2; ^(0.5); x 例1-4设 75 , 24= - =b a,计算|) tan(| |) | | sin(| b a b a + + 的值。 (-24)/180*; 75/180*; a1(a); b1(b); ();

(a11)/((c))^(0.5) 例1-5 设三角形三边长为2,3,4===c b a ,求此三角形的面积。 432; ()/2; (p*()*()*())^(0.5) 例1-7 设 ??????????=101654321A ,?? ??? ?????-=112311021B ,计算||,,A AB B A +,1-A 。 [1,2,3;4,5,6;1,0,1]; [-1,2,0;1,1,3;2,1,1]; ; *b; (a); (a); 例1-8 显示上例中矩阵A 的第2行第3列元素,并对其进行修改. [1,2,3;4,5,6;1,0,1]; (2,3); a(2,3)(' ') 例1-9 分别画出函数x x y cos 2 =和x x z sin = 在区间[-6π,6π]上的图形。 1; 1/6*:0.01:1/6*; (x.*x).*(x); (x); (); 例1-10 试求方程组??????????=????????? ?--432201624121X 的解。 [1,2,1;4,261,0,2]; [2;3;4]; (a)*b 例1-11 试求矩阵方程??????=????? ?????--111321201624121X 的解。 [1,2,1;4,261,0,2]; [1,2,3;1,1,1]; *(a)

MATLAB数学实验100例题解

一元函数微分学 实验1 一元函数的图形(基础实验) 实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧. 初等函数的图形 2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势. 解:程序代码: >> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps); plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象: 程序代码: >> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps); plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象: cot(x) 4在区间]1,1[-画出函数x y 1 sin =的图形. 解:程序代码: >> x=linspace(-1,1,10000); y=sin(1./x); plot(x,y); axis([-1,1,-2,2]) 图象:

二维参数方程作图 6画出参数方程???==t t t y t t t x 3cos sin )(5cos cos )(的图形: 解:程序代码: >> t=linspace(0,2*pi,100); plot(cos(t).*cos(5*t),sin(t).*cos(3*t)); 图象: 极坐标方程作图 8 作出极坐标方程为10/t e r =的对数螺线的图形. 解:程序代码: >> t=0:0.01:2*pi; r=exp(t/10); polar(log(t+eps),log(r+eps)); 图象: 90270 分段函数作图 10 作出符号函数x y sgn =的图形. 解:

数学软件MATLAB实验报告 实验八

实验八:概率论与数理统计的MATLAB 实现 实验目的与要求: 能运用MATLAB 提供的针对概率统计课程的工具箱。 实验内容: 1、用normpdf函数计算正态概率密度函数。 该函数的调用格式为:Y=normpdf(X,MU,SIGMA) 2、用normpdf函数计算正态分布的分布函数。 该函数的调用格式为:F=normcdf(X,MU,SIGMA) 3、用chi2inv函数计算卡方分布的分布函数的逆函数。 分布函数的逆函数及其调用格式:x=chi2inv(P,v) 4、随机取8只活塞环,测得他们直径为(以mm计): 74.001 74.005 74.003 74.001 74.000 73.998 74.006 74.002 。 设环直径的测量值服从正态分布,现估计总体的方差2 程序代码: x=[74.001 74.005 74.003 74.001 74.000 73.998 74.006 74.002]; p=mle('norm',x); sigma2hatmle=p(2)^2 5、从一批灯泡中随机的取5只做寿命试验,测得寿命(以小时计)为: 1050 1100 1120 1250 1280 设灯泡寿命服从正态分布,求灯泡寿命平均值的95%置信区间。 程序代码: x=[1050 1100 1120 1250 1280]; [p,ci]=mle('norm',x,0.05) 6、下面列出的是某工厂随机选取的20只部件的装配时间(分): 9.8 10.4 10.6 9.6 9.7 9.9 10.9 11.1 9.6 10.2

10.3 9.6 9.9 11.2 10.6 9.8 10.5 10.1 10.5 9.7 设装配时间的总体服从正态分布,标准差为0.4,是否可以认为装配时间的均值在0.05的水平上不小于10. 0H :10<μ vs 1H :10≥μ 程序: %正态总体的方差已知时的均值检验 x1=[9.8 10.4 10.6 9.6 9.7 9.9 10.9 11.1 9.6 10.2]; x2=[10.3 9.6 9.9 11.2 10.6 9.8 10.5 10.1 10.5 9.7]; x=[x1 x2]'; m=10;sigma=0.4;a=0.05; [h,sig,muci]=ztest(x,m,sigma,a,1) 因此,在0.05的水平下,可以认为装配时间的均值不小于10。 7、某种电子元件的寿命x (以小时计)服从正态分布,2 δμ和均未知。现测得16只元件的寿命如下: 159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170 问是否有理由认为元件的平均寿命大于225(小时)? 0H :225≤μ vs 1H :225>μ 程序: %正态总体的方差求知时的均值检验 x=[159 280 101 212 224 379 179 264 222 362 168 250 149 260 485 170]; m=225;a=0.05; [h,sig,muci]=ttest(x,m,a,1)

实验一 Matlab软件的使用..

实验一 Matlab软件的使用 1、实验目的和要求 ㈠掌握Matlab软件的使用方法; ㈡常用离散时间信号的产生、显示和运算。 2、实验内容 ㈠Matlab软件常用命令和工具箱的使用,学会简单的编程; ㈡编程实现常用离散时间信号; ㈢离散时间信号的叠加、移位、线性卷积等基本运算。 3、实验原理 利用软件生成数字信号处理系统中所涉及的信号及信号的加、移位、卷积等运算。 4、主要仪器设备 计算机,Matlab软件。 实验步骤: 程序1: n=-5:30; x=[zeros(1,5),1,zeros(1,30)]; stem(n,x,'fill'); grid on 程序2: n=-5:30; x=[zeros(1,5),ones(1,8),zeros(1,23)]; stem(n,x,'fill'); grid on; 程序3: n=-15:20; x=[zeros(1,15),1,2,3,4,5,zeros(1,16)]; stem(n,x,'fill'); grid on; 程序4: N=10; M=10; L=N+M-1; x=[1,2,3,4,5,zeros(1,5)]; h=[1,2,1,2,zeros(1,6)]; y=conv(x,h); n=0:L-1; stem(n,y,'*k'); grid on ;

实验结果:1、 2、 3、 4、

实验二 DTFT和Z变换的Matlab实现 1、实验目的和要求 ㈠掌握DTFT和Z变换的Matlab实现的方法; ㈡进一步了解离散时域系统的时频特性。 2、实验内容 ㈠ DTFT和z变换的计算和基本性质; ㈡离散时域系统的仿真; ㈢利用函数实现离散时域系统的频率响应,分析零、极点的分布,进一步了解离散时域系统的时域特性。 3、实验原理 利用Matlab软件计算出系统函数的零极点分布、单位脉冲响应和系统频率响应等的图像并于笔算结果进行比较,找出异同。编译合适程序能计算取值范围不同的离散卷积。 4、主要仪器设备 计算机,Matlab软件。 实验步骤: 程序1: n=0:30; x=zeros(1,length(n)); h=zeros(1,length(n)); x([find((n>=0)&(n<=4))])=1; h([find((n>=0)&(n<=8))])=0.5; figure(1) subplot(3,1,1); stem(n,x); axis([0,30,0,2]); title('输入序列'); xlabel('n'); ylabel('x(n)'); subplot(3,1,2); stem(n,h); axis([0,30,0,2]);

南邮MATLAB数学实验答案(全)

第一次练习 教学要求:熟练掌握Matlab 软件的基本命令和操作,会作二维、三维几何图形,能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字,教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令,画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位(1-9班)或4位(10班以上) 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =,求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB程序设计实验报告

MATLAB 程序设计实验报告 一、实验目的 1. 通过实验熟悉MATLAB 仿真软件的使用方法; 2. 掌握用MATLAB 对连续信号时域分析、频域分析和s 域分析的方法,利用绘图命令绘制出典型信号的波形,了解这些信号的基本特征; 3. 掌握用MATLAB 对离散信号时域分析、频域分析和z 域分析的方法,利用绘图命令绘制出典型信号的波形,了解这些信号的基本特征; 4. 通过绘制信号运算结果的波形,了解这些信号运算对信号所起的作用。 二、实验设备 1. 计算机 : 2. MATLAB R2007a 仿真软件 三、实验原理 对系统的时域分析 信号的时域运算包括信号的相加、相乘,信号的时域变换包括信号的平移、反折、倒相及信号的尺度变换。 (1)信号的相加和相乘:已知信号)(1t f 和)(2t f ,信号相加和相乘记为 )()(1t f t f =)(2t f +;)()(1 t f t f =)(2t f *。 (2)信号的微分和积分:对于连续时间信号,其微分运算是用diff 函数来完成的,其语句格式为:diff(function,’variable’,n),其中function 表示需要进行求导运算的信号,或者被赋值的符号表达式;variable 为求导运算的独立变量;n 为求导的阶数,默认值为求一阶导数。连续信号的积分运算用int 函数来完成,语句格式为:diff(function,’variable’,a,b),其中function 表示需要进行被积信号,或者被赋值的符号表达式;variable 为求导运算的独立变量;a,b 为积分上、下限,a 和b 省略时为求不定积分。 (3)信号的平移、翻转和尺度变换 信号的平移包含信号的左移与右移,信号的翻转包含信号的倒相与折叠,平移和翻转信号不会改变信号)(t f 的面积和能量。信号的尺度变换是对信号)(t f 在时间轴上的变化,可使信号压缩或扩展。)(at f 将原波形压缩a 倍,)/(a t f 将原波形扩大a 倍。 ¥ 对系统频率特性的分析

Matlab 实验

一、实验目的 熟悉MATLAB 软件的用户环境;了解MATLAB 软件的一般命令;掌握MATLAB 向量、数组、矩阵操作与运算函数;掌握MATLAB 软件的基本绘图命令;掌握MATLAB 语言的几种循环、条件和开关选择结构,及其编程规范。 通过该实验的学习,使学生能灵活应用MATLAB 软件解决一些简单问题,能借助MATLAB 软件的绘图功能,对函数的特性进行探讨,广泛联想,大胆猜想,发现进而证实其中的规律。 二、实验仪器、设备或软件: 电脑,MATLAB 软件 三、实验内容 1.MATLAB 软件的数组操作及运算练习; 2.直接使用MATLAB 软件进行作图练习; 3.用MATLAB 语言编写命令M 文件和函数M 文件。 四、实验步骤 1.在D 盘建立一个自己的文件夹; 2.开启软件平台——MATLAB ,将你建立的文件夹加入到MATLAB 的搜索路径中; 3.利用帮助了解函数max, min, sum, mean, sort, length ,rand, size 和diag 的功能和用法; 4.开启MATLAB 编辑窗口,键入你编写的M 文件(命令文件或函数文件); 5.保存文件(注意将文件存入你自己的文件夹)并运行; 6.若出现错误,修改、运行直到输出正确结果; 7.写出实验报告,并浅谈学习心得体会。 五、实验要求与任务 根据实验内容和步骤,完成以下具体实验,要求写出实验报告(实验目的→问题→算法与编程→计算结果或图形→心得体会) 1. 已知矩阵??????????=321212113A , ???? ??????--=101012111B 要求:(1)屏幕输出A 与B ;(2)A 的转置A′;(3)求A+B 的值;(4)求A-B 的值;(5)求4A ;(6)求A×B ;(7)求A -1.

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令,通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πx)(λ为非负实数)进行了分岔与混沌的研究,试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150

plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];

for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图

2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏,长满了草。他要将一头牛拴在牛栏边界的桩栏上,但只让牛吃到一半草,问拴牛鼻子的绳子应为多长? 2.2.2问题分析 如图所示,E为圆ABD的圆心,AB为拴牛的绳子,圆ABD为草场,区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半,可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积,2a÷π×πx2=2aπ2,再求AB的面积,用扇形ABE的面积减去三角形ABE的面积即可。

实验1答案 Matlab软件环境的基本使用

实验1 Matlab的基本运算(1) 实验目的: 1、了解Matlab的发展和主要功能; 2、熟悉Matlab工作环境的各个窗口; 3、掌握建立矩阵的方法; 4、掌握Matlab各种表达式的书写规则以及常用函数的使用。 实验内容: 1.在命令窗口依次输入下列命令,根据执行结果分析其功能: help 在命令窗口输入后,会看到帮助的目录信息,如下图所示。 lookfor inv:输入后会在命令窗口显示帮助中含有字母inv的函数和函数功能的列表。 help inv:查找函数inv,并显示其帮助信息。

which inv:显示函数保存的位置信息,和类型。 2.建立自己的工作目录,再将自己的工作目录设置到Matlab搜索路径下。用help命令查 询自己的工作目录;

步骤:File菜单——Set Path——Add Folder——选择自建文件夹——确定——点击save ——点击close。 3.利用Matlab的帮助功能,分别查询inv、plot、max、round等函数的功能及用法; help inv

help plot help max help round

4. 先求下列表达式的值,然后显示Matlab 工作空间的使用情况并保存全部变量。 (1) sin z e = +12 2851 (2 ) ln(z x =+212,其中.i x +??=??-?? 2120455 (3)...sin(.)ln ,.,.,.,...,.,.,.a a e e a z a a --+= ++=---0303303 0330292828293022 提示:利用冒号表达式生成a 向量,求各点的函数值时用点乘运算。

数学软件MATLAB实验报告 实验一

实验一:了解数学软件MATLAB 实验目的与要求: 了解MATLAB的基本特点,掌握MATLAB界面上的主要窗口,熟练掌握MATLAB的帮助系统。 实验内容: 1启动按钮 打开MATLAB主界面以后,单击“Start”按钮,显示一个菜单,利用“START”菜单机器子菜单中的选项,打开MATLAB的有关工具。 2命令窗口 命令窗口(Command Window)是用于输入数据,运行MA TLAB函数和脚本并显示结果的主要工具之一。命令窗口没有打开时,从“Desktop”菜单中选择“Command Window”选项可以打开它。 在命令窗口中键入命令并执行:a=[123;456;789] 在上述语句末尾加分号“;”,结果是什么?请与不加分号的情况作比较。

功能。 命令历史窗口(command history)显示命令窗口中最近输入的所有语句。先关闭历史窗口,再分别用“Desktop”菜单打开它和用command history命令打开它。 (1)将命令历史窗口中的语句复制到命令窗口中; (2)直接双击命令历史窗口中的语句。 4工作空间窗口 清空工作空间的命令是:clear

清空命令窗口的命令是:clc 在命令窗口中键入: t=0:pi/4:2*pi y=sin(t) 在命令窗口中键入:who,看运行结果;

在命令窗口中键入:whos,看运行结果; 在命令窗口中键入:whos y,看运行结果。

退出MATLAB时,工作空间中的内容随之清除。可以将当前空间中的部分或全部变量保存到一个MA T文件中,它是一个二进制文件,扩展名为.mat。然后可以在以后使用它时载入它,请练习这一操作过程。 5帮助浏览器 使用帮助浏览器可以搜索和查询所有Math Works产品的文档和演示。帮助浏览器是集成到MA TLAB桌面的一个HTML查看器。请打开帮助浏览器,熟悉它。 分别用Help函数和doc函数获取format函数的帮助,进而说明format函数的功能

matlab数学实验复习题(有答案)

复习题 1、写出3个常用的绘图函数命令 2、inv (A )表示A 的逆矩阵; 3、在命令窗口健入clc 4、在命令窗口健入 clear 5、在命令窗口健入6、x=-1:0.2:17、det (A )表示计算A 的行列式的值;8、三种插值方法:拉格朗日多项式插值,分段线性插值,三次样条插值。 9、若A=123456789?? ???????? ,则fliplr (A )= 321654987?? ???????? A-3=210123456--??????????A .^2=149162536496481?????? ???? tril (A )=100450789?? ???????? triu (A ,-1)=123456089??????????diag (A )=100050009?? ???? ???? A(:,2),=258A(3,:)=369 10、normcdf (1,1,2)=0.5%正态分布mu=1,sigma=2,x=1处的概率 [t,x]=ode45(f,[a,b],x0),中参数的涵义是fun 是求解方程的函数M 文件,[a,b]是输入向量即自变量的围a 为初值,x0为函数的初值,t 为输出指定function 开头;17、二种数值积分的库函数名为:quad;quadl

4 3,4 21、设x )的功能是作出将X 十等分的直方图 22、interp1([1,2,3],[3,4,5],2.5) Ans=4.5 23、建立一阶微分方程组???+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来) 二、写出运行结果: 1、>>eye(3,4)=1000 01000010 2、>>size([1,2,3])=1;3 3、设b=round (unifrnd (-5,5,1,4)),则=3 5 2 -5 >>[x,m]=min(b);x=-5;m=4 ,[x,n]=sort(b) -5 2 3 5 4 3 1 2 mean(b)=1.25,median (b )=2.5,range (b )=10 4、向量b 如上题,则 >>any(b),all(b<2),all(b<6) Ans=1 0 1 5、>>[5 6;7 8]>[7 8;5 6]=00 11 6、若1234B ?? =?? ??,则 7、>>diag(diag(B))= 10 04 8、>>[4:-2:1].*[-1,6]=-4 12 9、>>acos(0.5),atan(1) ans= 1.6598 ans= 0.7448 10、>>norm([1,2,3]) Ans=3.3941 11、>>length ([1,3,-1])=3

数学实验三 MATLAB软件入门(绘图)

西安理工大学 学生实验报告 数学实验 实验课程名 称: 实验名称:实验三MATLAB软件入门(绘图)学院:自动化与信息工程学院学生姓名: 班级: 学号:

一、实验目的及意义 [1]掌握MATLAB软件的基本绘图命令; [2]掌握MATLAB语言的几种循环、条件和开关选择结构。 通过该实验的学习,使学生能灵活应用MATLAB软件解决一些简单问题,能借助MATLAB软件的绘图功能,对函数的特性进行探讨,广泛联想,大胆猜想,发现进而证实其中的规律。 二、实验内容 [1]使用MATLAB进行作图练习; [2]用MATLAB语言编写命令M文件。 三、实验心得体会 经过腾讯课堂视频教学与展示,以及多次练习,已经能够熟练掌握所学内容,通过MATLAB各种绘图函数的调用,解决平面及三维绘图,着实感受到MATLAB的方便与强大。 四、实验任务 1. 学习plot命令的使用 (1)采用plot命令绘制y=sin(x) -2π

(2)采用plot命令绘制y=cos(x) -2π

(4)采用一句plot命令将两个sin与cos的图形绘制在一张图上要求sin 采用蓝色实线; cos 采用红色点划线 Legend命令在图上给出图标 采用axis命令,是x坐标显示范围[-1,1] ,y坐标显示范围[-2,2]

(4)采用help命令学习plot命令的用法 >> help plot plot - 二维线图 此MATLAB 函数创建Y 中数据对X 中对应值的二维线图。如果X 和Y 都是向量,则它们的长度必须相同。plot 函数绘制Y 对X 的图。 如果X 和Y 均为矩阵,则它们的大小必须相同。plot 函数绘制Y 的列对X 的列的图。如果X 或Y 中的一个是向量而另一个是矩阵,则矩阵的各维中必须有一维与向量的长度相等。如果矩阵的行数等于向量长度,则plot

数学实验与数学软件(Mathmaticandmatlab)

数学软件与数学实验2013-2014学年度秋季学期期末试卷 专业:统计学 班级:11级2班 学号:20110723 姓名:晏静

一、按要求计算出下列表达式的值 (1)318, 3 162 53 ?? + ? ?? , 21 eπ+, 2.5 tg, 2 log15; (2)给出π的9位和e的10位近似值; (3)求658和4102的最大公约数及35和25的最小公倍数; (4)产生10个0与10之间随机数的一个表; (5)求虚数1453 i i i i +- -的实部,虚部,模,共轭,辐角。 (6)自己运用Table建立两个表,并进行表运算,如连接、并集、交、排序等操作。

二、因式分解 22212321332112322 1 22(1)()()()4;(2)21;x x x x x x x x x x x x x x x +++++---- 解: 三、解方程(组) 1234234124234-2+344-+-3(1)+31-73+3 x x x x x x x x x x x x x -=??=? ? +=??+=-? 65432(2)5232002000.x x x x x x -+--++= 四、求极限 () 20 (1)1sin ;(2);(3)56! ctg x n x n n n Lim x Lim n n →→∞ →∞++

(1) (2) (3) 五、求导数 32 22(1)()=ln(x+1+);(2)()=cos 2,; (3)=log (),Z . x f x x f f x e y x y Z xy x y y ???求的导数已知求求关于的二阶导 (1) (2) (3) 六、求下列定积分与不定积分: ()()()12201+sin ln 1+(1);(2);(3)sin (1+cos ) +1(1+)(2+-) x x dx dx x x x x x x ? ? ?2 2-(4)=0,=1,==.y D D x y y x I x e d σ??设是由直线围成的区域,计算的值 (1) (2)

实验一 Matlab软件初步(终稿)讲解

实验一 MATLAB软件初步 1.1实验目的 在计算机上进行数学实验需要一个功能完备的计算平台,本数学实验系列课程均在MATLAB软件环境中进行数学实验。而本实验的目的就是熟悉MATLAB软件使用环境,并掌握其基本应用,为后续实验课程奠定必要的计算机基础。 1.2 实验要求 本实验要求学生掌握MATLAB软件的基本功能,会使用MATLAB进行简单的数值运算,会使用plot函数绘制简单二维图形,会编写简单的MATLAB程序,会使用常用MATLAB帮助命令。 1.3知识概要 1.3.1 导言 MATLAB是 MATrix LABoratory 的缩写,是由美国MathWorks公司开发的工程计算软件。MATLAB语言是以数组为基本数据单位,包括控制流程语句、函数、数据结构、输入输出及面向对象等特点的高级语言,具有以下主要特点:1)运算符和库函数极其丰富,语言简洁,编程效率高。MATLAB除了提供和C语言一样的运算符外,还提供广泛的矩阵和向量运算符。 2)既具有结构化的控制语句(如for循环、while循环、break语句、if语句和switch语句),又有面向对象的编程特性。 3)图形功能强大。它既包括对二维和三维数据可视化、图像处理、动画制作等高层次的绘图命令,也包括可以完全修改图形局部及编制完整图形界面的、低层次的绘图命令。 1.3.2 MATLAB软件环境 启动MATLAB后,软件界面如下图所示,该界面主要包括命令窗口、工作区窗口、历史命令窗口、当前目录窗口。

图 1 MATLAB 主窗口 1 命令窗口(Command Window) 命令窗口是用户与MATLAB 进行交互的主要场所,命令窗口的空白区域,用于输入和显示计算结果,可以在该区域键入各种MATLAB 命令进行各种操作,键入数学表达式进行计算。 例如,当键入变量赋值命令:123x =+?,并回车, 将在命令行的下面显示: 再输入求三角正弦函数值的表达式sin()4 y =并回车,将显示:

matlab数学软件实验测试题

数学软件实验测试题 Matlab 作业电子版 姓名:** 学号:**

一、选择语句和循环语句编程 1. 用if 语句实现以下的计算,其中,a b 的值从键盘输入。 2lg cos cos ,52sin ln ,219,14x b ax b x x a y a b b x x e x ?+-≤<-? ? =++-≤=-2&x<1)y=a*sin(b)+log(abs(b+a/x)); else if ( -5<=x&x<-2)y=a*x^2+b*cos(x); else if (1<=x&x<=4) y=9^(log10(x))+exp(cos(b)); end end end 2. function y=f(n)

y=0; for i=1:n y=y+1/i.^4; end 二、函数和调用函数 22(,)(1,1)(2,2)(3,3)(100,100)f a b a b y f f f f =+=+++1. 已知,求的值。L (,)f a b y 要求:定义为函数,利用循环语句求。 代码如下: 1. function m=f(a,b) m1=0;m2=0; for a=1:100 for b=1:100 m2=m2+b*b end m1=m1+m2+a*a; end 三、画图和拟合 1. 已知1cos(4),2sin(2)1,11y x y x y x ==?-≤≤,完成下列操作: (1)在同一坐标系下用不同的颜色和线型绘制两条曲线。 (2)以子图形式绘制两条曲线。(以一行两列形式绘图) 2. 已知函数f (x )在[1,101]区间上10个整数采样点的函数值如实验表1所示: 实验表1 10个采样点的函数值 x 1 11 21 31 41 51 61 71 81 91 101

相关文档
最新文档