MATLAB 练习(含代码)
MATLAB编程练习(含答案很好的)
001双峰曲线图：
z=peaks(40);
mesh(z);
surf(z)
002解方程：
A=[3,4,-2;6,2,-3;45,5,4];
>> B=[14;4;23];
>> root=inv(A)*B
003傅里叶变换
load mtlb ;
subplot(2,1,1);
plot(mtlb);
>> title('原始语音信息');
>> y=fft(mtlb);
>> subplot(2,1,2);
>> yy=abs(y);
>> plot(yy);
>> title('傅里叶变换')
004输入函数：
a=input('How many apples\n','s')
005输出函数
a=[1 2 3 4 ;5 6 7 8;12 23 34 45;34 435 23 34]
a =
1 2 3 4
5 6 7 8
12 23 34 45
34 435 23 34
disp(a)
a =
1 2 3 4
5 6 7 8
12 23 34 45
34 435 23 34
b=input('how many people\n' ,'s')
how many people
two people
b =
two people
>> disp(b)
two people
>>
006求一元二次方程的根
a=1;b=2;c=3;
d=sqrt(b^2-4*a*c);
x1=(-b+d)/(2*a)
x1 =
-1.0000 + 1.4142i
>> x2=(-b-d)/(2*a)
matlab编程实例100例
1-32是：图形应用篇
33-66是：界面设计篇
67-84是：图形处理篇
85-100是：数值分析篇
实例1：三角函数曲线（1）function shili01
h0=figure('toolbar','none',...
'position',[198 56 350 300],...
'name','实例01');
h1=axes('parent',h0,...
'visible','off');
x=-pi:0.05:pi;
y=sin(x);
plot(x,y);
xlabel('自变量X');
ylabel('函数值Y');
title('SIN( )函数曲线');
grid on
实例2：三角函数曲线（2）function shili02
h0=figure('toolbar','none',...
'position',[200 150 450 350],...
'name','实例02');
x=-pi:0.05:pi;
y=sin(x)+cos(x);
plot(x,y,'-*r','linewidth',1);
grid on
xlabel('自变量X');
ylabel('函数值Y');
title('三角函数');
实例3：图形的叠加
function shili03
h0=figure('toolbar','none',...
'position',[200 150 450 350],...
'name','实例03');
x=-pi:0.05:pi;
y1=sin(x);
y2=cos(x);
plot(x,y1,...
'-*r',...
MATLAB--练习(含代码)
MATLAB 练习
实验一常见分布的概率密度、分布函数生成
[实验目的]
1. 会利用MATLAB软件计算离散型随机变量的概率，连续型随机变量概率密度值。
2.会利用MATLAB软件计算分布函数值，或计算形如事件。
3.会求上分位点以及分布函数的反函数值。
[实验要求]
1.掌握常见分布的分布律和概率密度的产生命令，如binopdf,normpdf
2. 掌握常见分布的分布函数命令，如binocdf,normcdf
3. 掌握常见分布的分布函数反函数命令，如binoinv,norminv
[实验内容]
1 事件A在每次试验中发生的概率是0.3，计算
（1）在10次试验中A恰好发生6次的概率；
（2）在10次试验中A至多发生6次的概率.
binopdf(6,10,0.3)
>> binocdf(6,10,0.3)
2设随机变量X服从参数是3的泊松分布，求概率
poisspdf(6,3)
3设随机变量X服从区间[2,6]上的均匀分布，求
（1）X=4时的概率密度值；
（2）.
unifpdf(4,2,6)
>> unifcdf(5,2,6)
4设随机变量X服从参数是6的指数分布，求
（1）X=0,1,2,3,4,5,6时的概率密度值;
(2).
exppdf(0:6,6)
>> expcdf(5,6)
5设随机变量X服从均值是6，标准差是2的正态分布，求
(1) X=3,4,5,6，7,8,9时的概率密度值;
（2）X=3,4,5,6，7,8,9时的分布函数值；
（3）若=0.345,求x;
（4）求标准正态分布的上0.05分位数。
normpdf(3:9,6,2)
matlab练习题
matlab练习题
Matlab是一种高级的技术计算软件，被广泛应用于科学与工程领域。通过编写Matlab代码，我们可以解决各种数学问题、进行数据分析和
可视化、进行仿真等。本文将提供一些针对初学者的Matlab练习题，
帮助读者熟悉Matlab的基本语法和常用操作。
练习一：矩阵运算
在Matlab中，矩阵运算是一项重要的操作。请编写Matlab代码完
成以下任务：
1. 创建两个3×3的矩阵A和B，元素的值自行选择。
2. 计算矩阵A和B的和，并将结果存储在变量C中。
3. 计算矩阵A和B的乘积，并将结果存储在变量D中。
4. 计算矩阵A的转置，并将结果存储在变量E中。
提示：使用Matlab的矩阵运算符和内置函数可以简化代码编写过程。
练习二：函数编写与图形绘制
Matlab是一个功能强大的数学软件，可以用来编写自定义的函数并
进行图形绘制。
1. 编写一个Matlab函数，该函数可以计算任意输入整数n的阶乘，并将结果输出。
2. 调用该函数，计算并输出5的阶乘。
3. 使用Matlab绘制函数y=2x+1的图像，其中x的取值范围为-10
到10。
提示：使用Matlab的plot函数可以实现图形的绘制，并可使用grid
函数添加刻度线。
练习三：数据分析
Matlab提供了丰富的数据分析工具和函数，可以用来处理和分析实
验数据。
1. 编写一个Matlab脚本，实现对给定实验数据进行统计分析的功能。
2. 实验数据保存在文本文件data.txt中，每行包含一个数据点。读
取数据并将其存储在一个向量中。
3. 计算数据的均值、方差和标准差，并将结果输出。
matlab编程训练
matlab训练题目1、建立一个脚本文件，用于两个变量值的互相交换。2、比较从键盘输入的两个数的大小，输出其中较大的数。（扩展，输出两个同型矩阵中对应位置较大的数）3、编写一个函数求一元二次方程的根，并编写测试脚本文件测试结果。4、编写一个函数求圆的面积与周长，并编写测试脚本文件测试结果。5、编写一个函数用于求数n的阶乘（分别用循环与递归），再编写测试脚本文件求1！+2！+。。。+10！。6、编写脚本测试分段函数：--2x-1
7、有函数，编
写脚本文件测试该函数，x的值从键盘输入。 8、编写脚本文件，功能：从键盘输入一个字符，如果是大写字母，输出其后继字母；如果是小写字母，输出其前导字母，（‘z’的后继为’a’，‘a’的前导为’z’）；如果是数字字符，输出其ASCII码；其他字符原样输出。 9、某商场对顾客购买的商品总价打折，规则是500以下不打折，500～1000的95折，1000～3000的9折，3000～5000的85折，5000以上8折。从键盘输入总价，分别用if
语句和switch语句计算折扣，输出应付价。10、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：price<200 没有折扣200≤price<500 3%折扣500≤price<1000 5%折扣1000≤price<2500 8%折扣2500≤price<5000 10%折扣5000≤price 14%折扣输入所售商品的价格，求其实际
销售价格。 11、给出学生的百分制成绩，
matlab练习题
练习一：基本操作和简单语句输入
【学习要点】 认识MATLAB 的命令窗口，在命令窗口中输入简单语句，认识
“，”和“；”结束语句的区别，注意变量大小写的区别。
1．打开MATLAB 的命令窗口，键入help general 等帮助命令，
2．键入 demo 演示命令，
6．用who 命令观察工作空间中的变量，用clear 命令清除变量。
7．已知：x=π/4, y=2, 创建矩阵2345sin()9A x y ⎡⎤=⎢⎥⎣⎦
。 8．任意生成一个向量V ，计算V ×V ’和V ’×V 。
注意：三角函数用弧度计算；默认变量ans
matlab简单编程21个题目及答案
1、设
⎥
⎦
⎤
⎢
⎣
⎡
+
+
=
)
1(
sin
3
5.0
cos
2
x
x
x
y
，把x=0~2π间分为101点，画出以x为横坐
标，y为纵坐标的曲线。
第一题的matlab源程序：
①考虑cos（x）为一个整体，然后乘以中括号里面的全部
x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x).*(0.5+3*sin(x)./(1+x.^2)); %y的表达式
plot(x,y)%画出图形
图如下：
②考虑对整体求解cos，先求x乘以括号中的部分
x=0:2*pi/100:2*pi; %x的步长以及范围从0到2*pi y=cos(x.*(0.5+3*sin(x)./(1+x.^2))); %y的表达式
plot(x,y) %画出图形
图如下：
2、产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。
第二题的matlab源程序如下：
R1=randn(8,6) %产生正态分布随机矩阵
R1 =
1.0933 -0.7697 1.5442 -0.1924 1.4193 0.2157
1.1093 0.3714 0.0859 0.8886 0.2916 -1.1658
-0.8637 -0.2256 -1.4916 -0.7648 0.1978 -1.1480
0.0774 1.1174 -0.7423 -1.4023 1.5877 0.1049
-1.2141 -1.0891 -1.0616 -1.4224 -0.8045 0.7223
-1.1135 0.0326 2.3505 0.4882 0.6966 2.5855
matlab 代码例子
matlab 代码例子
标题：Matlab代码例子
1. 矩阵乘法
在Matlab中，可以使用'*'运算符来进行矩阵乘法。例如，下面的代码演示了如何计算两个矩阵的乘积：
```matlab
A = [1 2; 3 4];
B = [5 6; 7 8];
C = A * B;
disp(C);
```
2. 求解线性方程组
Matlab提供了求解线性方程组的函数`linsolve`，可以方便地求解形如Ax=b的线性方程组。下面的代码演示了如何使用`linsolve`函数求解一个线性方程组：
```matlab
A = [1 2; 3 4];
b = [5; 6];
x = linsolve(A, b);
disp(x);
```
3. 绘制函数图像
Matlab的绘图功能非常强大，可以绘制各种函数的图像。下面的代码演示了如何绘制一个简单的函数图像：
```matlab
x = linspace(0, 2*pi, 100);
y = sin(x);
plot(x, y);
```
4. 随机数生成
Matlab提供了多种生成随机数的函数，下面的代码演示了如何生成一个服从正态分布的随机数向量：
```matlab
mu = 0;
sigma = 1;
n = 100;
x = mu + sigma * randn(n, 1);
disp(x);
```
5. 求解最优化问题
Matlab的优化工具箱提供了多种求解最优化问题的函数，下面的代码演示了如何使用`fmincon`函数求解一个无约束的最优化问题：
```matlab
fun = @(x) x(1)^2 + x(2)^2;
matlab练习题及答案详解(附程序)
Matlab练习题及答案详解
一、基础题
1、下表为
1)
2) 把任意两点距离的表中，自己到自己的距离，转变成无穷大
clc
clear
close all
data=xlsread('ti1.xls');%从excel里导入数据
x=data(1,:);%将data里第一行的所有列赋值给x
y=data(2,:);%将data里第二行的所有列赋值给x
n=length(x);%表x的长度
dis=zeros(n,n);%对距离表进行初始化
for i=1:n
for j=1:n
dis(i,j)=sqrt((x(i)-x(j))^2+(y(i)-y(j))^2);%平面坐标距离公式
end
end
dis(dis==0)=inf;%将距离表中为0的数值，变成无穷大
redis=dis;
2、按要求步骤，操作下表
要求：
1) 对每列进行[0,1]区间化'(min())/(max min())
=--；
x x x x
2) 把区间化后的表的右边3列，连接在左边3列下边，构成16行3列的表；
3) 把16行3列进行转置，得到3行16列的表。然后对列进行操作，把3行
16列的表变成一列，后转置成一行，得到最终的结果。
clc
clear
close all
x=load('ti2.txt');%从txt里导入数据
big=max(x);small=min(x);%算出每列的最大最小值
[m,n]=size(x);%表x的行和列的大小
newx1=zeros(m,n);%对newx1进行初始化，大小和初始x一样大
for i=1:n
newx1(:,i)=(x(:,i)-small(i))./(big(i)-small(i));
Matlab例题
例15 在同一坐标中，绘制3个同心圆，并加坐标控制。 在同一坐标中，绘制3 程序如下： function circle(R) theta=0:0.01:2*pi; x=R*sin(theta); y=R*cos(theta); plot(x,y),grid on,hold on; axis equal clear;clc for R=1:2:20 circle(R) end
例14 在0≤x≤2p区间内，绘制曲线y1=2e-0.5x和 0≤x≤2p区间内，绘制曲线y1=2e-0.5x和 y2=cos(4πx)，并给图形添加图形标注。 y2=cos(4πx)，并给图形添加图形标注。 程序如下： x=0:pi/100:2*pi; y1=2*exp(y1=2*exp(-0.5*x); y2=cos(4*pi*x); plot(x,y1,x,y2),grid on title('x from 0 to 2{\pi}'); 2{\ %加图形标题 %加图形标题 xlabel('Variable X'); %加 %加X轴说明 ylabel('Variable Y'); %加 %加Y轴说明 text(0.8,1.5,'曲线y1=2e^{text(0.8,1.5,'曲线y1=2e^{-0.5x}'); %在指定位置添加 %在指定位置添加 图形说明 text(2.5,1.1,'曲线y2=cos(4{\ text(2.5,1.1,'曲线y2=cos(4{\pi}x)'); legend(‘y1’ y2’ legend(‘y1’,’y2’) %加图例 %加图例
matlab程序练习
hjl 笔记
1.将存于E 盘根目录的图片按顺序运动 for i=1:3
k=int2str(i);
k1=strcat('e:\',k,'.jpg'); #读取存于E 盘根目录的后缀名为.jpg 的图片按顺序运动# a1=imread(k1); image(a1);
m(:,i)=getframe; end
movie(m,10)
2.随机生成一个6行10列矩阵，矩阵元素在0与1之间。如果那个元素大于等于0.5，则把这个元素变为1；如果哪个元素小于0.5，则把这个元素变为0。最后把结果显示出来 Ra=rand(6,10); R=Ra; for i=1:6
for j=1:10
if Ra(i,j)<0.5 Ra(i,j)=0; else
Ra(i,j)=1; end end end R Ra
3.使用程序绘制颜色渐变的图形 for j=1:30 for i=1:20 a(i,j)=j end end
image(a)
5
10
15
20
25
30
2468101214161820
for i=1:20 for j=1:30 a(i,j)=i end end
image(a)
5
10
15
20
25
30
2468101214161820
r=1;
for t=0:0.01:2*pi
for m=0:0.01:2*pi
x=r*cos(m)*cos(t); y=r*sin(m)*cos(t); z=r*sin(t); hl=plot3(x,y,z);
set(hl,'color',[z,z,z]); end hold on; end
lighting none
matlab编程实例100例
1-32是：图形应用篇
33-66是：界面设计篇
67-84是：图形处理篇
85-100是：数值分析篇
实例2：三角函数曲线（2）
function shili02
h0=figure('toolbar','none',...
'position',[200 150 450 350],...
'name','实例02');
x=-pi:0.05:pi;
y=sin(x)+cos(x);
plot(x,y,'-*r','linewidth',1);
grid on
xlabel('自变量X');
ylabel('函数值Y');
title('三角函数');
实例3：图形的叠加
function shili03
h0=figure('toolbar','none',...
'position',[200 150 450 350],...
'name','实例03');
x=-pi:0.05:pi;
y1=sin(x);
y2=cos(x);
plot(x,y1,...
'-*r',...
x,y2,...
'--og');
grid on
xlabel('自变量X');
ylabel('函数值Y');
title('三角函数');
实例4：双y轴图形的绘制function shili04
h0=figure('toolbar','none',...
'position',[200 150 450 250],...
'name','实例04');
x=0:900;a=1000;b=0.005;
y1=2*x;
y2=cos(b*x);
[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot'); axes(haxes(1))
matlab总练习题(完整版)
>> vpa('pi',20)
ans =
3.85
>> vpa('exp(1)',20)
ans =
2.54
>> x=linspace(-pi,pi,21);
>> y=sin(x)
y =
Columns 1 through 6
-0.0000 -0.3090 -0.5878 -0.8090 -0.9511 -1.0000 Columns 7 through 12
-0.9511 -0.8090 -0.5878 -0.3090 0 0.3090
Columns 13 through 18
0.5878 0.8090 0.9511 1.0000 0.9511 0.8090 Columns 19 through 21
0.5878 0.3090 0.0000
>> k=1:1000;
>> kk=1./k;
>> kkk=kk./k;
>> res=sum(kkk)-(pi^2)/6
res =
-9.9950e-04
%承接上题
>> sum(kk)-log(1000)
ans =
0.5777
>> power(1+eps,1/eps) ans =
2.7183
>> a=rand(2,3)
a =
0.8147 0.1270 0.6324
0.9058 0.9134 0.0975 >> x=a(1,:)
x =
0.8147 0.1270 0.6324 >> y=a(2,:)
y =
MATLAB部分课后习题代码
习题1.6.2 第5题
clear
clc
A=[1 4 7 10;2 5 8 11;3 6 9 12]
B=flipud(A') %将A左旋90度C=fliplr(A') %将A右旋90度
习题2.6 第6题
clear
clc
x=input('输入10个两位整数')
b=mean(x);
for i=1:10
if x(i) < b & rem(x(i),2)~=0
x(i)
end
end
习题2.6 第7题
clear
clc
x=input('请输入一个正整数')
for i=2:x
if rem(x,i)==0
break
end
end
if i < x & i > 1
disp('不是素数')
else
disp('是素数')
end
3.4.2 自己练第6题
clear
clc
x=-3:0.01:3;
y=exp(-x.^2);
plot(x,y)
3.4.2 自己练第8题
clear
clc
x=-2:0.1:2;
y=-1:0.05:1;
[x,y]=meshgrid(x,y);
z=sqrt(5-x.^3/3-y.^2/7);
surf(x,y,z)
grid on
3.5习题第1题
（1）
clear
clc
x=-1:0.01:1;
y=log((1-x)./(1+x));
plot(x,y)
（2）
clear
clc
x=-pi/2:pi/100:pi/2;
y=sqrt(cos(x));
plot(x,y)
（3）
clear
clc
x=-1:0.03:1;
y=sin(1./x);
plot(x,y)
（4）
clear
clc
x=-0.5:0.01:0.5;
matlab练习题及答案详解附程序
Matlab练习题及答案详解1、下表为8组坐标
x 4.22
■
7.92
6.
56 8.4 9 6.79 '7.43 6.55 7. 06
y 9.16 ! 9.59
0. 36 9.3 4 7.58 3.92 1.71 0. 32
1）算出任意两点间的空间距离，得到任意两点距离表；
2）把任意两点距离的表中，自己到自己的距离，转变成无穷大
clc clear
close all
data=xlsread( 'ti1.xls' );% 从excel 里导入数据
x=data(1,:);%将data里第一行的所有列赋值给x
y=data(2,:);%将data里第二行的所有列赋值给x
n=length(x);% 表x 的长度
dis=zeros( n,n);% 对距离表进行初始化
for i=1: n
for j=1: n
dis(i,j)=sqrt((x(i)-x(j))A2+(y(i)-y(j))A2);% 平面坐标距离公式
end
end
dis(dis==0)=inf;% 将距离表中为0的数值，变成无穷大
redis二dis;
2、按要求步骤，操作下表
基础题
要求：
1)对每列进行[0,1]区间化x' (x min(x))/(max min(x))；
2)把区间化后的表的右边3列，连接在左边3列下边，构成16行3列的表；
3)把16行3列进行转置，得到3行16列的表。然后对列进行操作，把3行16列的表变成一列，后转置成一行，得到最终的结果。
clc
clear
close all
x=load('ti2.txt');% 从txt 里导入数据
MATLAB练习
MATLAB 练习
练习1 计算36sin 5e ++π的结果
练习2 输入矩阵
并删除矩阵A 的第一行。
练习 3 输入矩阵 ,分析C=(A
练习4 随便输入一个七阶矩阵A ，
(1)生成由它的2,4,6行和3，4，5列产生的子矩阵B
(2)生成由它的2,4,7行和4，1列产生的子矩阵C
(3)编程求出A 矩阵的最大元素
练习 5 设可逆方阵为A ，编写同时求 的M 函数文件。
练习6 编写一函数的M 文件,此函数表达式为：
并计算
练习7 分别用for 和while 语句编写一M 文件，计算1！+2！+⋯+50！的值
练习8 有一函数y xy x y x f 2sin ),(2++=，编写一函数M 文件输入自变量的值 输出相应函数值
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=162ln 973sin 56231
A ⎥⎦
⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=2211,2121B A ⎪⎪⎩
⎪⎪⎨⎧>+≤<-≤≤-<=30ln sin 3010113101121)(x x x x x x x x x x f ',,,12A A A A -)10(),30(),2(),2.0(πf f f f
练习9编一个转换成绩等级的函数文件，其中成绩等级转换标准为考试成绩分数在
[90，100]分显示优秀；在[80，90)分显示良好；在[60，80)分显示及格；在
[0，60)分显示不及格。
练习10在同一个图中画出衰减振荡曲线及它的包络线
t的取值范围是
]
4,0[π，给出标题，两坐标轴及曲线说明。
练习11 画出椭圆
1
2
52
2
2
2
=
+
y
x
的曲线图。
练习12在同一平面的两个窗口中画出函数在-3 ≤ x ,y≤ 3
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
[Ex,Dx]=normstat(1,sqrt(9)); >> [Ey,Dy]=normstat(0,sqrt(16)); >> rxy=-0.5; >> syms x y z >> z=x/3+y/2; >> covxy=rxy*sqrt(Dx)*sqrt(Dy); >> Ez=Ex/3+Ey/2 >> Dz=Dx/9+Dy/4+2*1/3*1/2*(covxy) >> covxz=Dx/3+covxy/2 实验四 统计中的样本数字特征 23 随机生成 4 组 10 个整数数据，求每组数据的平均值。 > X=fix(20*rand(10,4)) >> M=mean(X) 24 随机生成服从标准正态分布的 6 组 10 个数据，求每组数据的极差、样本方差、样本标准 差。 X=normrnd(0,1,10,6) M1=range(X) M2=var(X) M3=std(X) 实验五 两个正态总体均值差，方差比的区间估计 [实验目的] 1 掌握两个正态总体均值差，方差比的区间估计方法 2 会用 MATLAB 求两个正态总体均值差，方差比的区间估计
MATLAB 练习 实验一 常见分布的概率密度 、分布函数生成 [实验目的] 1. 会利用 MATLAB 软件计算离散型随机变量的概率，连续型随机变量概率密度值。 2.会利用 MATLAB 软件计算分布函数值，或计算形如事件。 3.会求上分位点以及分布函数的反函数值。 [实验要求] 1.掌握常见分布的分布律和概率密度的产生命令，如 binopdf,normpdf 2. 掌握常见分布的分布函数命令，如 binocdf,normcdf 3. 掌握常见分布的分布函数反函数命令，如 binoinv,norminv [实验内容] 1 事件 A 在每次试验中发生的概率是 0.3，计算 （1）在 10 次试验中 A 恰好发生 6 次的概率； （2）在 10 次试验中 A 至多发生 6 次的概率. binopdf(6,10,0.3) >> binocdf(6,10,0.3) 2 设随机变量 X 服从参数是 3 的泊松分布，求概率 poisspdf(6,3) 3 设随机变量 X 服从区间[2,6]上的均匀分布，求 （1）X=4 时的概率密度值； （2）. unifpdf(4,2,6) >> unifcdf(5,2,6) 4 设随机变量 X 服从参数是 6 的指数分布，求 （1）X=0,1,2,3,4,5,6 时的概率密度值; (2). exppdf(0:6,6) >> expcdf(5,6) 5 设随机变量 X 服从均值是 6，标准差是 2 的正态分布，求 (1) X=3,4,5,6，7,8,9 时的概率密度值; （2）X=3,4,5,6，7,8,9 时的分布函数值； （3）若=0.345,求 x; （4）求标准正态分布的上 0.05 分位数。 normpdf(3:9,6,2) normcdf(3:9,6,2) norminv(0.345,6,2) norminv(0.95,0,1) 6 设随机变量 X 服从自由度是 6 的 t 分布 ,求
[实验要求] 两个正态总体的区间估计理论知识 [实验内容] 25 为比较甲乙两种型号步枪子弹的枪口速度， 随机地取甲种型号子弹 10 发， 得到枪口速度 的平均值为 500（m/s）,标准差 1.10（m/s） ；乙种型号子弹 20 发，得到枪口速度的平均值 为 496（m/s）,标准差 1.20（m/s） 。两个总体都近似服从正态分布，且方差相等。求两个总 体均值差的置信水平为 0.95 的置信区间。 t=tinv(0.975,28); >> s=sqrt((9*1.10^2+1.20^2)/28); >> d1=(500-496)-t*s*sqrt(1/10+1/20) >> d2=(500-496)+t*s*sqrt(1/10+1/20) 26 从甲乙两个蓄电池厂生产的产品中，分别抽取 10 个产品，测得它们的电容量（单位： Ah）为： 甲厂：146 141 138 142 140 143 138 137 142 137 乙厂：141 143 139 139 140 141 138 140 142 136 若蓄电池的电容量服从正态分布， 求两个工厂生产的蓄电池的电容量方差之比的置信水平为 0.90 的置信区间。 x1=[146 141 138 142 140 143 138 137 142 137]; >> x2=[141 143 139 139 140 141 138 140 142 136]; >> s1=var(x1); >> s2=var(x2); >> f1=finv(0.95,9,9); >> f2=finv(0.05,9,9); >> d1=s1/s1*(1/f1) >> d2=s1/s2*(1/f2) 实验六 假设检验 [实验目的] 1 会用 MATLAB 进行单个正态总体均值及方差的假设检验 2 会用 MATLAB 进行两个正态总体均值差及方差比的假设检验 [实验要求] 熟悉 MATLAB 进行假设检验的基本命令与操作 [实验内容] 27 某工厂生产 10
>> plot(x,y) 13 设随机变量 X 服从均值是 6，标准差是 2 的正态分布。 （1） 画出 X 的概率密度图形 x=-10:0.01:10; >> y=normpdf(x,6,2); >> plot(x,y) （2） 画出 X 的分布函数图形 > x=-10:0.01:10; >> y=normcdf(x,6,2); >> plot(x,y) （3） 在同一个坐标系中画出均值是 6，标准差是 1，,3 的正态分布概率密度图形 x=0:0.01:13; >> y1=normpdf(x,6,1); >> y2=normpdf(x,6,2); >> y3=normpdf(x,6,3); >> plot(x,y1,x,y2,x,y3) 14 设随机变量 X 服从自由度是 6 的 t 分布 （1）画出 X 的概率密度图形 x=-10:0.01:10; >> y=tpdf(x,6); >> plot(x,y) （2）画出 X 的分布函数图形 x=-10:0.01:10; >> y=tcdf(x,6); >> plot(x,y) 15 设随机变量 X 服从自由度是 6 的 （1） 画出 X 的概率密度图形 x=0:0.01:10; >> y=chi2pdf(x,6); >> plot(x,y) （2）画出 X 的分布函数图形 x=0:0.01:10; >> y=chi2cdf(x,6); >> plot(x,y) 16 设随机变量 X 服从第一自由度是 2，第,二自由度是 6 的 F 分布 （1）画出 X 的概率密度图形 x=0:0.001:10; >> y=fpdf(x,2,6); >> plot(x,y) （2）画出 X 的分布函数图形 x=0:0.001:10; >> y=fcdf(x,2,6);
[实验要求] 1.掌握 MATLAB 画图命令 plot 2.掌握常见分布的概率密度图像和分布函数图像的画法 [实验内容] 9 事件 A 在每次试验中发生的概率是 0.3，记 10 次试验中 A 发生的次数为 X. (1)画出 X 的分布律图形； x=0:10; >> y=binopdf(x,10,0.3); >> plot(x,y,'.') >> （2）画出 X 的分布函数图形； x=0:0.01:10; >> y=binocdf(x,10,0.3); >> plot(x,y) 10 设随机变量 X 服从参数是 6 的指数分布， （1）画出 X 的概率密度图形 x=0:0.01:10; >> y=exppdf(x,6); >> plot(x,y) （2）画出 X 的分布函数图形 x=-1:0.01:10; >> y=expcdf(x,6); >> plot(x,y) 11 设随机变量 X 服从参数是 3 的泊松分布。 (1)画出 X 的分布律图形； x=0:10; >> y=poisspdf(x,3); >> plot(x,y,'.') （2）画出 X 的分布函数图形； >> x=0:0.01:10; >> y=poisscdf(x,3); >> plot(x,y) 12 设随机变量 X 服从区间[2,6]上的均匀分布。 （1）画出 X 的概率密度图形 x=0:0.01:10; >> y=unifpdf(x,2,6); >> plot(x,y,'*') （2）画出 X 的分布函数图形 >> x=0:0.01:10; >> y=unifcdf(x,2,6);
（1）X=-3,-2,-1,0,1,2,3 时的概率密度值; （2）X=-3,-2,-1,0,1,2,3 时分布函数值; （3）若=0.345,求 x; （4）求 t 分布的上 0.05 分位数. tpdf(-3:3,6) tcdf(-3:3,6) tinv(0.345,6) tinv(0.95,6) 7 设随机变量 X 服从自由度是 6 的 (1) X=0,1,2,3,4,5,6 时的概率密度值; (2) X=0,1,2,3,4,5,6 时的分布函数值; (3) 若=0.345,求 x; (4) 求
>> pro(2)=poisscdf(1,0.8); >> pro(3)=poisscdf(4,0.8)-pro(2); >> pro(1)=1-pro(2)-pro(3) >> Ey=pro*price' 22 设随机变量 Xห้องสมุดไป่ตู้N(1,9),Y~N(0,16) ， 且 X 与 Y 的相关系数为， 令 Z=X/3+Y/2.求 （1） E(Z) ,D(Z) ; (2) 求 X 与 Z 的相关系数 。 解 根据题意，有 E(X)=1,D(X)=9,E(y)=0 ,D(y)=16 由 E(Z)=E(X/3)+E(y/2)=E(X)/3+E(Y)/2 得 E(Z)=1/3+0=1/3 由 ，cov(X,Y)= 有 Cov(X,Y)=0.5=6 D(Z)=D(X/3+Y/2)=D(X/3)+D(Y/2)+2cov(X/3,Y/2) 得到 D(Z)=D(X)/9+D(Y)/4+2cov((X,Y) =9/9+16/4+2=3 得 cov（X,Z）=cov(X,X/3)+cov(X,Y/2) =1/3cov(X,X)+1/2cov(X,Y) =9/3+(-6)=0
2
分布 ,求
2
分布的上 0.05 分位数.
chi2pdf(0:6,6) chi2cdf(0:6,6) chi2inv(0.345,6) chi2inv(0.95,6) 8 设随机变量 X 服从第一自由度是 2，第,二自由度是 6 的 F 分布 ,求 (1) X=0,1,2,3,4,5,6 时的概率密度值; (2) X=0,1,2,3,4,5,6 时的分布函数值; (3) 若=0.345,求 x; (4) 求 F 分布的上 0.05 分位数. fpdf(0:6,2,6) fcdf(0:6,2,6) finv(0.345,2,6) finv(0.95,2,6) 实验二 概率作图 [实验目的] 1.熟练掌握 MATLAB 软件的关于概率分布作图的基本操作 2.会进行常用的概率密度函数和分布函数的作图 3.会画出分布律图形
2
分布
>> plot(x,y) 实验三 数字特征 [实验目的] 1 加深对数学期望，方差的理解 2 理解数学期望，方差的意义，以及具体的应用 3 加深对协方差，相关系数的理解 4 了解协方差，相关系数的具体的应用 [实验要求] 1 概率与频率的理论知识，MATLAB 软件 2 协方差，相关系数的理论知识，MATLAB 命令 cov,corrcoef [实验内容] 17 若 ~B(20,0.3), 求 E(X),D(X). [M,,V]=binostat(20,0.3) 18 随机变量 X 的概率密度为，求 E(X),D(X). Syms X f1=x; f2=2-x; Ex=int(xf1,0,1)+int(xf2,1,2); Ex2=int(x^2f1,0,1)+int(x^2f2,1,2); Dx=Ex2-Ex^2 19 设(X,Y)的概率密度为, 求 EX,EY. Syms x y fxy=8xy; Ex=int(int(fxyx,y,0,x),x,0,1) Ey=int(int(fxyy,y,0,x),x,0,1) 20（续 19）求 cov（X,Y） Syms x y fxy=8xy; Ex=; Ey=; Cxy=int(int(fxy(x-Ex)(y-Ey),y,0,x),x,0,1) 21 某种商品每件表面上的疵点数 X 服从泊松分布，平均每件上有 0.8 个疵点。若规定表面 不超过一个疵点的为一等品，价值 10 元，表面疵点数大于一个不多于 4 个的为二等品，价 值 8 元。表面疵点数多于 4 个则为废品，求产品价值的均值。 解 设 X 表示产品表面上的疵点数，由已知，EX=0.8,且 X 服从泊松分布，故 EX=0.8, = 设 Y 表示产品价值，则 Y 有分布律： ==0.8088 =0.1898 =0.0014 故有 EY=0(元) > pro=[ ]; >> price=[0 10 8];