1Matlab基础汇总.pptVIP

%求素数的程序 function a=avg(x) a=sum(x(:))/length(x(:)); 1999 B题 钻井布局 已知旧井位坐标如下, 求网格状的新井位的位置, 使得旧井的利用率最高. i 1 2 3 4 5 6 ai? 0.50 1.41 3.00 3.37 3.40 4.72 bi? 2.00 3.50 1.50 3.51 5.50 2.00 i 7 8 9 10 11 12 ai? 4.72 5.43 7.57 8.38 8.98 9.50 bi? 6.24 4.10 2.01 4.50 3.41 0.80 将a, b去整数部分后的图形如下 接下来就可以在[0,1]× [0,1]区域内寻找一个位置P, 使得P点与这些点的距离小于给定的ε的个数最多. 一个比较简单的方法是使用离散穷举的方法来搜索, 即在[0,1]× [0,1]区域内均匀地取n2个点, 计算每个点能够利用的旧井数, 求出最大的即可. n越大, 结果越准确, 但计算量也越大; 反之亦然. clear clc a=[0.5 1.41 3.00 3.37 3.4 4.72 4.72 5.43 7.57 8.38 8.89 9.50]; b=[2.00 3.50 1.50 3.51 5.50 2.00 6.24 4.10 2.01 4.50 3.41 0.80]; e=0.05; A=a-floor(a); %取小数部分 B=b-floor(b); %取小数部分 n=0; x1=0; y1=0; for x=0:0.01:1 for y=0:0.01:1 n1=sum(abs(A-x)e abs(B-y)e); %计算与点(x,y)距离小于e的点的个数 if(n1n) %如果点的个数比原来的更多, 则更新点数及位置 n=n1; x1=x; y1=y; end end end n [x1 y1] plot(A,B,*) %用点画出数据图 hold on plot(x1,y1,+r); %用+画出搜索到的位置 plot([x1-e,x1+e,x1+e,x1-e,x1-e],[y1-e,y1-e,y1+e,y1+e,y1-e],r); axis([0 1 0 1]); %设置坐标范围 hold off 运行结果: n = 4 ans = 0.3600 0.4600 在2008年国赛的A题中, 有一个问题是如果确定下图中的各黑点的中心. 一个近似的方法是用“重心”来代替中心. 首先将图片的黑点独立出来, 然后将黑色像素的坐标进行平均, 这样得到的平均值就是黑点的“重心”. A=imread(GrayeImage.bmp); %读取灰度图像, 保存在当前目录下 A=im2bw(A,0.5); %变为二值图像 A=~A; %将图像数据取反，结果为1表示黑色，0表示白色 [B,n]=bwlabeln(A); %将图像分成n块 CG=zeros(2,n); %作用是分配内存 for i=1:n [Y,X]=find(B==i); %获取第i块的坐标 CG(1,i)=mean(X); %计算第i块的重心x坐标 CG(2,i)=mean(Y); %计算第i块的重心y坐标 end CG imshow(~A); %作图 hold on plot(CG(1,:),CG(2,:),+r); hold off 运行结果 作函数 的图像. 画出一个中心为(0,0)、长半轴为2、短半轴为1的椭圆及其外切矩形. 产生一个5×10的矩阵, 其元素为0或1的随机数. 矩阵初等变换. (1) 交换矩阵两行; (2) 某一行乘以一个常数; (3) 某行乘以一个常数加到另一行. 5. 用蒙特卡罗法求多圆体的周长. 已知n个圆的圆心及半径, 求由这n个圆组成的复合体的周长. xi 0.4173 0 0.9027 0.9448 0.4909 yi 0.4893 0.3377 0.9001 0.3692 -0.3 ri 0.7803 0.3897 0.2417 0.4039 0.45 《数学实验》主要内容分3部分； 第1部分是软件操作实验， 并对MATLAB命令的使用格式作了较为详细的说明； 第2部分是应用篇，由12个各自独立的实验组成， 其目的是让学生学会利用MATLAB软件解决一些 简单的实际问题； 第3部分是数学建模实验，由12个独立的综合建模实验 组成，内容涉及生活、经济、管理、科技等多个学科的建模 问题，题材广泛，具有一定的实用价值和趣味性。 最后，附录中对Lingo数学