matlab实验九、十、十一
Matlab第十一讲--数据的统计分析
输入：name——概率分布类型、P（概率）——数据向 量、A——分布参数 输出：Y——逆累积分布向量
name同前
计算均值为0，方差为1的正态分布， x = icdf(‘norm',0.1:0.2:0.9,0,1)
Matlab相关命令介绍
Matlab相关命令介绍
mle 系列函数：参数估计
[phat,pci]=mle(‘name’,X,alpha)
load 从matlab数据文件中载入数据
S=load('数据文件名') 如果数据格式是XXXX.mat ，可以直接 load XXXX; 如果文本格式XXXX.txt，也可以用load载入，load 'XXXX.txt'; 另外文本格式也可以通过Import data转换成.mat格式， matlab默认处理.mat格式数据！
name+stat 系列函数：均值与方差函数
数字特征的相关函数
var 方差
1 n 2 2 var( X ) s ( xi X ) n 1 i 1
%若X为向量，则返回向量的样本方差。A为矩阵，返回A的列向 量的样本方差构成的行向量。 std 标准差
n n
1 1 2 2 ( xi X ) 或者 n ( xi X ) n 1 i 1 i 1
Matlab相关命令介绍
name+rnd 系列函数：随机数发生函数 random(‘name’,A1,A2,A3,m,n) %(通用函数）
normrnd(1:6,1./(1:6)) normrnd([1 2 3;4 5 6],0.1,2,3) %mu为均值矩阵 normrnd(10,0.5,[2,3]) %mu为10，sigma为0.5的2 行3列个正态随机数 y=random('norm',2,0.3,3,4) %产生12（3行4列） 个均值为2，标准差为0.3的正态分布随机数
matlab数学实验练习题
Matlab 数学实验实验一 插值与拟合实验内容：预备知识：编制计算拉格朗日插值的M 文件。
1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。
通过数值和图形输出，将三种插值结果与精确值进行比较。
适当增加n ，再做比较，由此作初步分析。
下列函数任选一种。
（1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s10≤≤-=x x y（4）、22),exp(2≤≤--=x x y2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为)(0)()(t eV V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。
试由下面一组t ，V 数据确定0V 和τ。
实验二 常微分方程数值解试验实验目的：1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法；2. 掌握用微分方程模型解决简化的实际问题。
实验内容：实验三地图问题1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量：以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。
根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土的近似面积，并与它的精确值41288km2比较。
实验四狼追兔问题狼猎兔问题是欧洲文艺复兴时代的著名人物达.芬奇提出的一个数学问题。
当一个兔子正在它的洞穴南面60码处觅食时，一只恶狼出现在兔子正东的100码处。
当两只动物同时发现对方以后，兔子奔向自己的洞穴，狼以快于兔子一倍的速度紧追兔子不放。
狼在追赶过程中所形成的轨迹就是追击曲线。
高等数学：MATLAB实验
MATLAB实验
2.fplot绘图命令 fplot绘图命令专门用于绘制一元函数曲线,格式为:
fplot('fun',[a,b]) 用于绘制区间[a,b]上的函数y=fun的图像.
MATLAB实验 【实验内容】
MATLAB实验
由此可知,函数在点x=3处的二阶导数为6,所以f(3)=3为 极小值;函数在点x= 1处的二阶导数为-6,所以f(1)=7为极大值.
MATLAB实验
例12-10 假设某种商品的需求量q 是单价p(单位:元)的函 数q=12000-80p,商 品的总成本C 是需求量q 的函数 C=25000+50q.每单位商品需要纳税2元,试求使销售 利润达 到最大的商品单价和最大利润额.
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验 实验九 用 MATLAB求解二重积分
【实验目的】 熟悉LAB中的int命令,会用int命令求解简单的二重积分.
MATLAB实验
【实验M步A骤T】 由于二重积分可以化成二次积分来进行计算,因此只要
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
MATLAB实验
实验七 应用 MATLAB绘制三维曲线图
【实验目的】 (1)熟悉 MATLAB软件的绘图功能; (2)熟悉常见空间曲线的作图方法.
【实验要求】 (1)掌握 MATLAB中绘图命令plot3和 mesh的使用; (2)会用plot3和 mesh函数绘制出某区间的三维曲线,线型
Matlab程序设计与应用第二版刘卫国课后实验答案
实验一：T1：%%第一小题z1=2*sin(85*pi/180)/(1+exp(2))%%第二小题x=[2,1+2i;,5];z2=1/2*log(x+sqrt(1+x.^2));z2%%第三小题a=::;z3=1/2*(exp*a)-exp*a)).*sin(a++log(+a )/2)%%第四题t=0::z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t. ^2-1)+(t>=2&t<3).*(t.^2-2*t+1)T2：A=[12,34,-4;34,7,87;3,65,7]B=[1,3,-1;2,0,3;3,-2,7]disp ('A+6*B=');disp(A+6*B);disp('A-B+I=');disp(A-B+eye(3));disp('A*B=');disp(A*B);disp('A.*B=');disp(A.*B);disp('A^3=');disp(A^3);disp('A.^3=');disp(A.^3);disp('A/B=');disp(A/B);disp('B\A=');disp(B\A);disp('[A,B]=');disp([A,B]);disp('[A([1,3],:);B^2]=');disp([A([1,3],:);B^2]);T3:z=1:25;A=reshape(z,5,5)';B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];C=A*BD=C(3:5,2:3)T4-1:a=100:999;b=find(rem(a,21)==0);c=length(b)T4-2:a=input('请输入一个字符串：','s');b=find(a>='A'&a<='Z');a(b)=[];disp(a);实验二：T1：E=eye(3),R=rand(3,2),O=zeros(2,3),S=d iag([1,2]);A=[E,R;O,S]disp('A^2=');disp(A^2);disp('[E,R+RS;O,S^2]');B=[E,R+R*S;O,S^2]T2：H=hilb(5)P=pascal(5)Hh=det(H)Hp=det(P)Th=cond(H)Tp=cond(P)a=abs(Th-1);b=abs(Tp-1);if a>bdisp('帕萨卡矩阵P性能更好'); elseif a<bdisp('希尔伯特矩阵H性能更好'); elsedisp('两个矩阵性能相同');endT3:a=1:25;A=reshape(a,5,5)disp('行列式的值：');disp(det(A));disp('矩阵的秩：');disp(rank(A));disp('矩阵的迹：');disp(trace(A));disp('矩阵的范数：');disp(norm(A));T4:A=[-29,6,18;20,5,12;-8,8,5][V,D]=eig(A)T5:A=[1/2,1/3,1/4;1/3,1/4,1/5;1/4,1/5,1/ 6]B=[,,]'X1=A\BB(3)=X2=A\Bdisp('系数矩阵A的条件数：');disp(cond(A));T6:a=1:25;A=reshape(a,5,5)disp('B1=sqrtm(A)');B1=sqrtm(A)disp('B2=sqrt(A)');B2=sqrt(A)disp('B1*B1');B1*B1disp('B2.*B2');B2.*B2实验三：T1:x=:2:for i=1:length(x);if (x(i)<0)&(x(i)~=-3)y(i)=x(i)^2+x(i)-6;elseif(x(i)>=0)&(x(i)<5)&(x(i)~=2)&(x(i)~=3 );y(i)=x(i)^2-5*x(i)+6;elsey(i)=x(i)^2-x(i)-1;endendyT2-if:s=input('please enter the score:');while (s<0||s>100)disp('the score is not reasonable'); s=input('please enter the score:'); endif (s>=90&s<=100);disp('A');elseif(s>=80&s<90);disp('B');elseif(s>=70&s<80);disp('C');elseif(s>=60&s<70);disp('D');elsedisp('E');endT2-switch:s=input('please enter the score:'); while (s<0||s>100)disp('the score is not reasonable'); s=input('please enter the score:'); endswitch fix(s/10)case{9,10}disp('A');case{8}disp('B');case{7}disp('C');case{6}disp('D');case{0,1,2,3,4,5}disp('E');endT3:t=input('请输入工时time=');if t>120w=120*84+(t-120)**84;elseif t<60w=t*84-700;elsew=84*t;enddisp('应发工资为：');disp(w);T4:a=10+floor(rand(1)*89)b=10+floor(rand(1)*89)s=input('请输入+或-或*或/','s');while(s~='+'&&s~='-'&&s~='*'&&s~='/') disp('输入的符号错误，请重新输入'); s=input('请输入+或-或*或/','s'); endswitch scase{'+'}c=a+b;case{'-'}c=a-b;case{'*'}c=a*b;case{'/'}c=a/b;endcT5:A=rand(5,6)n=input('请输入n:');while (n<1)disp('输入的n有误，请重新输入'); n=input('请输入n:');endif n>5B=A(5,:)elseB=A(n,:)End实验四：T1:n=input('请输入n:');m=0;for i=1:n;m=m+1/i^2;endpii=sqrt(6*m)T1-2:n=input('请输入n:');m=1./(1:n).^2;s=sum(m);pii=sqrt(6*s)T2:n=1000;y=0;for i=1:n;y=y+1/(2*i-1);if y>=3disp('最大n值为：');i-1disp('对应的y值为：');y-1/(2*i-1)break;endendT3:format longa=input('请输入a:');b=input('请输入b:');x(1)=a/+b);i=1;x(2)=a/(x(i)+b);while (abs(x(i+1)-x(i))>&i<500)i=i+1;x(i+1)=a/(x(i)+b);endr1=(-b+sqrt(b^2+4*a))/2;r2=(-b-sqrt(b^2+4*a))/2;disp(['x',num2str(i+1),'=',num2str(x( i+1),8)]);disp(['x',num2str(i),'=',num2str(x(i) ,8)]);disp(['r1=',num2str(r1,8),',r2=',num2 str(r2,8)]);format shortT4:f(1)=1;f(2)=0;f(3)=1;i=4;while(i<=100)f(i)=f(i-1)-2*f(i-2)+f(i-3);i=i+1;endfmax=max(f)fmin=min(f)fsum=sum(f)fp=length(find(f>0))fn=length(find(f<0))f0=length(find(f==0))T5:j=0;s=0;for i=3:50m=i*(i-1)-1;if isprime(m)==1ms=s+m;j=j+1;endenddisp(['亲密数的对数：',num2str(j)]); disp(['所有亲密素数之和：',num2str(s)]);实验五：fn1:function f=fn1(n)f=n+10*log(n.^2+5);fn2:function f=fn2(n)f=0;for i=1:nf=f+i*(i+1);endfx:function f=fx(x)f=1./((x-2).^2++1./((x-3).^4+;T1:x=input('请输入一个复数：');fe=exp(x)flg=log(x)fsin=sin(x)fcos=cos(x)T2:m1=input('请输入m1:');m2=input('请输入m2:');t=input('请输入theta:');A=[m1*cos(t),-m1,-sin(t),0;m1*sin(t), 0,...cos(t),0;0,m2,-sin(t),0;0,0,-cos(t),1 ]B=[0,m1*,0,m2*]'disp('[a1,a2,N1,N2]=');disp([A\B]');T3:j=1;n=1;k=1;for i=10:99 % 挑选出10~99中所有的素数if isprime(i)==1m(j)=i;j=j+1;endendfor t=1:length(m);s(n)=10*rem(m(t),10)+fix(m(t)/10); %挑选出的素数进行位置交换if isprime(s(n))==1 %再挑选出交换绝对素数r(k)=m(t);k=k+1;endn=n+1;enddisp(['所有两位绝对素数：',num2str(r)])T5:y1=fn1(40)/(fn1(30)+fn1(20)) %函数调用y2=fn2(40)/(fn2(30)+fn2(20))实验六：T1:x=linspace(0,2*pi,101);y=+3*sin(x)./(1+x.^2)).*cos(x);plot(x,y);T2:x=linspace(0,2*pi,100);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(2,2,1);plot(x,y1,'r:',x,y2,'y-',x,y3,'b-.'); subplot(2,2,2);plot(x,y1,'r:');subplot(2,2,3);plot(x,y2,'y-');subplot(2,2,4);plot(x,y3,'b-.');T2-3:x=linspace(0,2*pi,10);y1=x.^2;y2=cos(2*x);y3=y1.*y2;subplot(3,4,1);bar(x,y1,'r');axis([0,7,0,40]); subplot(3,4,2);stairs(x,y1,'y');axis([0,7,0,40]); subplot(3,4,3);stem(x,y1,'b');axis([0,7,0,40]); subplot(3,4,4);fill(x,y1,'g');axis([0,7,0,40]); subplot(3,4,5);bar(x,y2,'r');axis([-1,7,-1,1]); subplot(3,4,6);stairs(x,y2,'y');axis([-1,7,-1,1]); subplot(3,4,7);stem(x,y2,'b');axis([-1,7,-1,1]); subplot(3,4,8);fill(x,y2,'g');axis([-1,7,-1,1]); subplot(3,4,9);bar(x,y3,'r');axis([0,7,-30,50]); subplot(3,4,10);stairs(x,y3,'y');axis([0,7,-30,50]); subplot(3,4,11);stem(x,y3,'b');axis([0,7,-30,50]); subplot(3,4,12);fill(x,y3,'g');axis([0,7,-30,50]);T3:x=linspace(-5,5,100);y=[];for x0=x;if x0<=0;y=[y,(x0+sqrt(pi))/exp(2)]; elsey=[y,1/2*log(x0+sqrt(1+x0^2))];endendsubplot(1,2,1);plot(x,y);for i=1:length(x)if x(i)<=0y(i)=(x(i)+sqrt(pi))/exp(2); elsey(i)=1/2*log(x(i)+sqrt(1+x(i)^2)); endendsubplot(1,2,2);plot(x,y);T4:a=input('请输入a:');b=input('请输入b:');n=input('请输入n:');theta=0::2*pi;y=a.*sin(b+n.*theta);polar(theta,y);T5:x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-sqrt(x.^2+y.^2 )/4);subplot(1,2,1);surf(x,y,z);subplot(1,2,2);contour3(x,y,z);T6:s=0:pi/2;t=0:3*pi/2;[s,t]=meshgrid(s,t);x=cos(s).*cos(t);y=cos(s).*sin(t);z=sin(s);surf(x,y,z);shading interp;实验八：T1:A=rand(1,30000);disp('均值为：');mean(A)disp('标准方差为：');std(A)disp('最大元素为：');max(A)disp('最小元素为：');min(A)k=find(A>;b=length(k);disp('大于的随机数个数占总数的百分比：');sprintf('%f%%', b/300)T2:P=45+50*rand(100,5)disp('最高分及序号：');[Y,U]=max(P,[],1)disp('最低分及序号：');[Y1,U1]=min(P,[],1)disp('每门课平均分：');mean(P,1)disp('每门课的标准方差：');std(P,0,1)s=sum(P,2);disp('总分最高分：');[Y2,U2]=max(s,[],1)disp('总分最低分：');[Y3,U3]=min(s,[],1)[zcj,xsxh]=sort(s,1,'descend') %desc end降序T3:h=6:2:18;t1=[,,,,,,];t2=[,,,,,,];t=:2:;T1=interp1(h,t1,t,'spline')T2=interp1(h,t2,t,'spline')T4:x=1:10:101;y=log10(x);p=polyfit(x,y,5)x1=1::101;y1=log10(x1);y2=polyval(p,x1);plot(x1,y1,'r-',x1,y2,'b-.');T5:P1=[1,2,4,0,5];P2=[0,0,0,1,2];P3=[0,0,1,2,3];P4=conv(P2,P3);sp1=length(P1);sp2=length(P4);P=[zeros(1,sp2-sp1),P1]+P4x=roots(P)A=[-1,,;,2,;0,5,];Ax=polyval(P,A)AA=polyvalm(P,A)实验九：T1:f=inline('det([x x^2 x^3;1 2*x 3*x^2;0 2 6*x])');i=1;x=1;while x<=g(i)=f(x);i=i+1;x=x+;enddx=diff(g)/;dx(1)dx(101)dx(length(g)-1)T2:f1=inline('sqrt(cos(t.^2)+4*sin(2*t). ^2+1)')f2=inline('log(1+x)./(1+x.^2)')I1=quad(f1,0,2*pi)I2=quad(f2,0,2*pi)T3:A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0 ,2];B=[-4,13,1,11]';x1=A\B[L,U]=lu(A);x2=U\(L\B)[Q,R]=qr(A);x3=R\(Q\B)fun:function F=fun(X);x=X(1);y=X(2);z=X(3);F(1)=sin(x)+y^2+log(z)-7;F(2)=3*x+2^y-z^3+1;F(3)=x+y+z-5;T5:f=inline('3*x+sin(x)-exp(x)');fzero(f,X=fsolve('fun',[1 1 1]',optimset('Display','off'))实验十：T1:x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y))T2:syms x y;t=sym('5135');factor(x^4-y^4)factor(t)T3:syms beta1 beta2 x;simplify(sin(beta1)*cos(beta2)-cos(be ta1)*sin(beta2))simplify((4*x^2+8*x+3)/(2*x+1))T4:syms a b c d e f g h k;p1=[0,1,0;1,0,0;0,0,1];p2=[1,0,0;0,1,0;1,0,1];A=[a,b,c;d,e,f;g,h,k];B=p1*p2*AB1=inv(B)B1*B;tril(B)det(B)T5:syms x t a yf1=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/sin(x)^3;limit(f1,x,0) %(1)f2=(sqrt(pi)-sqrt(acos(x)))/sqrt(x+1) ;limit(f2,x,-1,'right') %(2)f3=(1-cos(2*x))/x;diff(f3,x,1) %(3)diff(f3,x,2)A=[a^x,t^3;t*cos(x),log(x)];diff(A,x,1) &(4)diff(A,t,2)diff(diff(A,x,1),t,1)f5=(x^2-2*x)*exp(-x^2-y^2-x*y);yx=-diff(f5,x,1)/diff(f5,y,1) %(5) fxy=diff(diff(f5,x,1),y,1)x=sym('0');y=sym('1');eval(fxy)T6:x=sym('x');f1=1/(1+x^4+x^8);int(f1,x)f2=1/(asin(x))^2/sqrt(1-x^2);int(f2,x)f3=(x^2+1)/(x^4+1);int(f3,x,0,inf)f4=exp(x)*(1+exp(x))^2;int(f4,x,0,log(2))实验十一：T1：syms n x ;S1=symsum(1/(2*n-1),n,1,10)S2=symsum(n^2*x^(n-1),n,1,inf)S3=symsum(n^2/5^n,n,1,inf)T2：x=sym('x');f=log(x);taylor(f,x,6,1)T3：x1=solve('log(1+x)-5/(1+sin(x))=2')x2=solve('x^2+9*sqrt(x+1)-1=0')x3=solve('3*x*exp(x)+5*sin(x)=0')[x4,y4]=solve('sqrt(x^2+y^2)-100=0',' 3*x+5*y-8=0')T4：dsolve('D2y+4*Dy+29*y=0','y(0)=0','Dy (0)=15','x')T5：[x,y,z]=dsolve('Dx=2*x-3*y+3*z','Dy=4 *x-5*y+3*z','Dz=4*x-4*y+2*z','t')。
MATLAB实验答案
实验一1.（1）z1=2*sin(85*pi/180)/(1+exp(2))（2）>> x=[2,1+2i;-0.45,5];>> z2=(1/2)*log(x+sqrt(1+x*x))(3)>> a=-3.0:0.1:3.0;>>z3=((exp(0.3*a)-exp(-0.3*a))/2).*sin(a+0.3)+log((0.3+a)/2)(4)t=0:0.5:2.5;>>z4=(t.*t).*(t>=0&t<1)+(t.*t-1).*(t>=1&t<2)+(t.*t-2.*t+1).*(t>=2& t<3)2.(1)>> A=[12,34,-4;34,7,87;3,65,7];>> B=[1,3,-1;2,0,3;3,-2,7];>> A+6*B(2)>> A*B(3)>> A^3(4)>> A/B(5)>> [A,B]3.(1)>>A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,2 4,25];>> B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];>> C=A*B(2)>> D=C(3:end,2:end)(3)Name valueclassA <5×5double> doubleB <5×3double> doubleC <5×3double> doubleD [520 397;705 557;…] double4.(1)>> X=100:999;>> length(find(~rem(X,21)))(2)>> y='ahHGJjskjJjsdsaSa';>> M=y(find(y<='z'&y>='a'))实验三1 x=input('请输入x的值')if x<0&x~=-3y=x^2+x-6;elseif x>=0&x<5&x~=2&x~=3 y=x^2-5*x+6;elsey=x^2-x-1;endy2（1）x=input('请输入成绩')if x>=90&x<=100disp('等级为A') elseif x>=80&x<=89disp('等级为B') elseif x>=70&x<=79disp('等级为C') elseif x>=60&x<=69disp('等级为D')elseif x>=0&x<=59disp('等级为E')elsedisp('成绩有误')end(2)x=input('请输入成绩')switch(x)case{90,100}disp('等级为A') case{80,89}disp('等级为B') case{70,79}disp('等级为C') case{60,69}disp('等级为D') case{0,59}disp('等级为E')otherwisedisp('成绩错误')end>>3Num=input('输入员工的工号:');Time=input('该工号员工的工作时数:'); salary=input('该员工底薪:');if Time>=120salary=salary*(1+0.15);elseif Time<=60salary=salary-700;elseif Time>60&Time<120salary=Time*84;endsalary=salary4:x=input('请输入x');y=input('请输入y');n=input('请输入运算符','s');if n=='+'m=x+y;elseif n=='-'m=x-y;elseif n=='*'m=x*y;elsem=x/y;endm5:a=rand(5,6)n=input('请输入要输出的第二行数');if n>5a(5,:)elsea(n,:)end实验四1：程序：n=input('输入n:');for a=1:n;b(a)=1/(a^2); endpi=sqrt(6*sum(b))2：.程序： sum=0;for n=1:200;b(n)=1/(2*n-1);sum=sum+b(n);if(sum>=3)break;endendn-14.程序：f(1)=1;max=0;f(2)=0;min=0;f(3)=1;c=0;zhengshu=0;fushu=0;zero=0;for n=4:100f(n)=f(n-1)-2*f(n-2)+f(n-3);endfor m=1:100if(f(m)>max)max=f(m);endif(f(m)<min)min=f(m);endc=sum(f);if(f>0)zhengshu=zhengshu+1; endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;endif(f(m)>0)zhengshu=zhengshu+1; endif(f(m)<0)fushu=fushu+1;endif(f(m)==0)zero=zero+1;enddisp('最大值')maxdisp('最小值')mindisp('和')cdisp('正数的个数')zhengshudisp('负数的个数')fushudisp('零的个数')zero5p=0;i=0;h=0;for m=2:49a=m;b=m+1;c=a*b-1;for k=1:cif rem(c,k)==0i=i+1;endendif i==2disp(['亲密数对是','(',num2str(a),',',num2str(b),')'])p=p+1;h=[h,c];endi=0;endsum=sum(h);disp(['亲密数对的个数是：',num2str(p)])disp(['亲密数对的和是：',num2str(sum)])实验五1，函数文件：文件名：Yunsuan程序：function[z,d,s,c]=yunsuan(a,b)z=exp(a+b*i);d=log(a+b*i);s=sin(a+b*i);c=cos(a+b*i);disp(['复数',num2str(a),'+',num2str(b),'i的指数为:',num2str(z),'，对数为：',num2str(d),'，正弦为：',num2str(s),'，余弦为：',num2str(c)])函数调用：程序：clearx=input('请输入复数实部的值');y=input('请输入复数虚部的值');[x,y]=yunsuan(x,y);2函数文件：文件名：zhu程序：function[p,q,w]=zhu(p,q,w)A=[p*cos(w),-p,-sin(w),0;p*sin(w),0,cos(w),0;0,q,-sin(w),0;0 ,0,-cos(w),1];B=[0;9.8*p;0;9.8*q];X=A\B函数调用：clearm1=input('m1=');m2=input('m2=');th=input('theta=');theta=pi*th/180;[m1,m2,theta]=zhu(m1,m2,theta); 5(1)function[f]=fn(n);f=n+10*log(n*n+5);[f40]=fn(40);[f30]=fn(30);[f20]=fn(20);y=f40/(f30+f20)(2)function[f]=fn(n);f=0;for i=1:nf=f+i*(i+1);end[f40]=fn(40);[f30]=fn(30);[f20]=fn(20);y=f40/(f30+f20)2_1clear;x=linspace(0,2*pi,101);y1=x.^2;y2=cos(2.*x);y3=y1.*y2;plot(x,y1,'b-',x,y2,'r:',x,y3,'g-.'); %y1蓝色实线，y2红色虚线，y3绿色点画线2_2subplot(2,2,1);%分四个子图(先画2行2列第1块) plot(x,y1);subplot(2,2,2);plot(x,y2),subplot(2,2,3);plot(x,y3);2_3 ()subplot(3,4,1); %y1的四种图形bar(x,y1);subplot(3,4,2);stairs(x,y1),subplot(3,4,3);stem(x,y1);subplot(3,4,4);fill(x,y1,'b');subplot(3,4,5); %y2bar(x,y2); %条形图subplot(3,4,6);stairs(x,y2), %阶梯图subplot(3,4,7);stem(x,y2); %杆图subplot(3,4,8);fill(x,y2,'b'); %填充图，注意必须加填充颜色subplot(3,4,9); %y3bar(x,y3);subplot(3,4,10);stairs(x,y3),subplot(3,4,11);stem(x,y3);subplot(3,4,12);fill(x,y3,'b');3clear;x=-5:0.1:5;if x<=0y=(x+sqrt(pi)/exp(2));elsey=0.5.*log(x+sqrt(1+x.^2));endplot(x,y);4M文件，假设文件名为Untitled6 a=input('a=');b=input('b='); %b单位为pi/4 b=b*pi/4;n=input('n=');q=linspace(-2*pi,2*pi,100);p=a*sin(b+n*q);plot(q,p);hold on; %保持图形命令窗口调用情况>> Untitled6a=1b=1n=1>> Untitled6a=2b=2n=25.程序：x=linspace(-5,5,21);y=linspace(0,10,31);[x,y]=meshgrid(x,y);z=cos(x).*cos(y).*exp(-0.25*sqrt(x.^2+y.^2)); subplot(1,2,1);title('surf(x,y,z)');surf(x,y,z);subplot(1,2,2);title('surfc(x,y,z)');surfc(x,y,z);实验七701clear allhf=figure('Color',[1,0,0],'KeyPressFcn','disp(''Left Button Pressed'')');702clear allx=linspace(-2*pi,2*pi,500);y=x.^2.*exp(2*x);h=plot(x,y);set(h,'Color','r','LineStyle',':','LineWidth',5)title('下图是y=x^2e^{2x}曲线的图像');703clear alla=linspace(-2*pi,2*pi,40);b=linspace(-2*pi,2*pi,40);[x,t]=meshgrid(a,b);v=10*exp(-0.01*x).*sin(2000*pi*t-0.2*x+pi);axes('view',[-37.5,30]);h=surface(x,t,v,'FaceColor','w','EdgeColor','flat');grid on;title('函数图像如下');set(h,'FaceColor','flat');705%实验七底层绘图操作clear all;[x,y,z]=cylinder(3,500);%cylinder是生成柱体的函数surf(x,y,z)title('实验七第五题圆柱体的光照和材料处理');xlabel('X-axis');ylabel('Y-axis');zlabel('Z-axis');axis([-5,5,-5,5,0,1])grid off;light('Color','r','Position',[-4,0,0],'Style','infinite');shading interp;material shiny;view(0,10);lighting phong;axis off;实验九程序：I1=quad('sqrt(cos(t.^2)+4*sin(2*t).^2+1)',0,2*pi) I2=quad('(log(1+x))./(1+x.*x)',0,1)程序：A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2]b=[-4;13;1;11]x=inv(A)*b实验十程序：x=sym('6');y=sym('5');z=(x+1)/(sqrt(3+x)-sqrt(y))分解因式（1）程序：syms x y;A=x^4-y^4;factor(A)（2）程序：factor(sym('5135'))3、化简表达式（1）程序：syms beta1 beta2y=sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2)simple(y)（2）程序：syms xy=(4*x^2+8*x+3)/(2*x+1)simple(y)4.(1)p1=sym(‘[0 1 0;1 0 0;0 01]’);p2= sym(‘[1 0 0;0 1 0;1 0 1]’);A= sym(‘[a b c;d e f;g h l]’);B=P1*P2*A(2)inv(B) (3)tril(B) (4)determ(B)5、用符号方法求下列极限或导数（1）程序：syms xf=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(sin(x)) limit(f)（2）程序：syms xy=(sqrt(pi)-sqrt(acos(x)))/(sqrt(x+1));limit(f,x,-1,'right')（3）程序：syms xy=(1-cos(2*x))/x;y1=diff(y)y2=diff(y,x,2)6、用符号方法求下列积分（1）程序：syms xf=1/(1+x^4+x^8)int(f)（2）程序：syms xf=1/(((asin(x))^2)*sqrt(1-x^2)) int(f)（3）程序：syms xf=(x^2+1)/(x^4+1)int(f,x,0,inf)（4）程序：syms xf=exp(x)*(1+exp(x))^2 y=int(f,x,0,log(2)) double(y)。
Matlab实验报告（1）
《Matlab语言与应用》课程实验报告*名：**班级：电信114学号：************指导老师：***二〇一三年十一月二十一日Matlab实验报告实验一一、实验问题求[12 + sin(2)×( 22 −4)]÷3^2的算术运算结果。
二、问题的分析该题主要熟悉Matlab环境下的对数学运算的熟悉，如何输入数据、建立函数输出结果。
三、上机实验结果如图四、实验的总结与体会通过本次实验，我学会了用Matlab来计算数学运算中的复杂技术。
我们也可以自己编写一个可以调用的函数，首先我们要了解Matlab语言函数的基本结构，结构如下：Function [返回变量列表]=函数名（输入变量列表）注释说明语句段，由%引导输入、返回变量格式的检测函数体语句例如：输入变量为k，返回的变量为m和s,其中s为前m项的和Function [m,s]=findsum(k)s=0;m=0;while(s<=k),m=m+1;s=s+m;end编写了函数，就可以将其存为findsum.m文件，这样就可以在Matlab环境中对不同的k值调用该函数了。
这样就可以灵活的实现我们想要的数据。
实验二一、实验问题二、问题分析输入矩阵时，空格或逗号表示间隔，分号表示换行，比如上面的矩阵A应写为A=[1,2,3;4,5,6;7,8,9]三、上机实验结果如图四、实验的总结与体会通过对本次上机实验了解到在Matlab中对矩阵的代数运算矩阵转置、矩阵的加减法运算、矩阵乘法、矩阵的左、右除、矩阵翻转、矩阵乘方运算、点运算等。
实验三一、实验要求画图，理解plot函数用法二、代码如下clear; clf;t=0:pi/20:2*pi;R=5;x=R*sin(t); y=R*cos(t);plot(x,y,'b:'), gridhold onrrr=[x;y;x+y];plot(rrr(1,:),'.','MarkerSize',10,'Color','r')plot(rrr(2,:),rrr(3,:),'o','MarkerSize',15,'Color','b'); axis([-8,20,-8,8]), % axis squarehold off三、Matlab运行结果如图实验四一、实验要求二、代码如下t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b') Grid三、Matlab运行结果如图实验五一、实验要求傅里叶频谱分析二、代码及分析(1)首先生成数据，包含50Hz和120Hz频率的正弦波x >>t = 0:.001:.25;>>x = sin(2*pi*50*t) + sin(2*pi*120*t);(2)再生成噪音信号yy = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')(3)对y进行快速傅里叶变换Y = fft(y,256);(4)计算功率谱Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('Frequency (Hz)')(5)只查看200Hz以下频率段plot(f(1:50),Pyy(1:50))title('Power spectral density')xlabel('Frequency (Hz)')三、上机结果如下图实验六一、实验要求FIR数字滤波器设计代码如下clear;close allf=[0,0.6,0.6,1]; m=[0,0,1,1]; % 设定预期幅频响应b=fir2(30,f,m); n=0:30; % 设计FIR 数字滤波器系数subplot(3,2,1),stem(n,b,'.')xlabel('n'); ylabel('h(n)');axis([0,30,-0.4,0.5]),line([0,30],[0,0])[h,w]=freqz(b,1,256);subplot(3,2,2),plot(w/pi,20*log10(abs(h)));gridaxis([0,1,-80,0]),xlabel('w/pi'); ylabel('幅度(dB)');二、上级结果如图实验七二、实验要求用guide实验一个简单的加减乘除计算器二、实验步骤在命令行输入guide命令，进入guide界面新建一个空白guide文件在空白文件中设置好功能模块如图模块建立好后，就要把编写好的加减乘除代码加入到各自的回调函数中，见下图两个被加数代码如下图加模块代码如下图减模块代码如下图乘模块代码如下图除模块代码如下图各模块的回调函数加完后就可以运行了，运行结果如下图总结：Matlab一个高级的距阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。
MATLAB与科学计算教学大纲
《MATLAB与科学计算》教学大纲一、课程的教学目标与任务数值计算与MA TLAB语言是信息科学与技术平台必修课程，注重锻炼学生的数学建模、分析能力等所需的基础知识和基本能力。
MA TLAB是一种以数值计算和数据图示为主的计算机软件，并包含适应多个学科的专业软件包，以及完善程序开发功能。
本课程要求学生掌握MATLAB的数据类型、矩阵输入和操作方法、语法结构、函数的使用以及二维、三维绘图功能，并能够熟练地将MATLAB应用于学习中，解决相关课程中的复杂的数学计算问题。
上机操作是本课程重要的教学环节，学生只有通过上机实习，才能领会MATLAB中众多功能，才能达到熟练应用的程度。
二、本课程与其它课程的联系MATLAB与科学计算是汽车检测与维修专业的公共基础课程，并能够熟练地将MA TLAB 应用于学习中，解决相关课程中的复杂的数学计算问题。
三、课程内容及基本要求（一）认识MATLAB （8学时）基本内容：MATLAB软件的发展历史，MATLAB的基本情况，MA TLAB启动和退出MATLAB，功能介绍，命令窗口，MATLAB帮助。
1．基本要求：了解MATLAB软件的发展历史，MATLAB的基本情况，以及学习的意义。
熟练掌握启动和退出MA TLAB的方法。
熟练掌握MATLAB的各种功能介绍。
熟练掌握命令窗口的使用。
熟练掌握MATLAB帮助。
2．重点、难点：重点：掌握MA TLAB软件基本操作。
难点：MA TLAB命令窗口的使用。
（二）MATLAB程序设计（4学时）基本内容：M文件的建立，调试。
程序流程语句，函数文件及编程技巧。
1．基本要求：了解M文件的建立，调试。
掌握程序流程语句，函数文件及编程技巧。
并能运用流程语句做简单程序设计。
2．重点、难点：重点：流程语句掌握。
难点：运用流程语句编写程序。
（三）MALAB绘图（6学时）基本内容：图形窗口；二维平面图形与坐标系；三维绘图。
1．基本要求：了解图形窗口。
matlab概率论部分数学实验指导书
1.9
0.8
1.1
0.1
0.1
4.4 5.5 1.6 4.6 3.4
0.7 -1.6 -0.2 -1.2 -0.1 3.4 3.7 0.8 0.0 2.0
试就下列两种情况分析这两种药物的疗效有无显示性的差异。 （ α = 0.05 ） 。 ① X 与 Y 的方差相同；② X 与 Y 的方差不同。 （7） 、 已知某一试验， 其温度服从正态分布， 现在测量了温度的五个值为： 1250， 1265，1245，1260，1275。问是否可以认为 µ = 1277 （8） 、其它教材上的题目或自己感兴趣的题目。 （ α = 0.05 ） 。 ？
A =[16 25 19 20 25 33 24 23 20 24 25 17 15 21 22 26 15 23 22
20 14 16 11 14 28 18 13 27 31 25 24 16 19 23 26 17 14 30 21 18 16 18 19 20 22 19 22 18 26 26 13 21 13 11 19 23 18 24 28 13 11 25 15 17 18 22 16 13 12 13 11 09 15 18 21 15 12 17 13 14 12 16 10 08 23 18 11 16 28 13 21 22 12 08 15 21 18 16 16
实验四、样本的统计与计算 实验目的： 熟练使用 matlab 对样本进行基本统计，包括样本的位置统计、分散性统计、样 本中心矩、分布的形状统计。求样本均值、中位数、样本方差，偏度、峰度、 样本分位数和其它数字特征，并能做出频率直方图和经验分布函数。 实验内容： 来自总体的样本观察值如下，计算样本的样本均值、中位数、样本方差、极差, 偏度、峰度、画出频率直方图，经验分布函数图。
matlab在控制方面的示例
一、简介MATLAB（Matrix Laboratory）是一种用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。
MATLAB被广泛应用于科学和工程领域，特别是在控制系统设计和模拟方面具有重要的作用。
在控制方面，MATLAB提供了丰富的工具和函数，可用于设计、分析和实现各种类型的控制系统，并且提供了许多示例来帮助用户更好地理解控制系统。
二、控制系统的建模和仿真1. 实例一：DC电机控制假设我们希望设计一个用于控制直流电机的系统。
我们可以使用MATLAB来建立直流电机的数学模型，并使用Simulink进行仿真。
通过编写方程或使用Simulink的模块化建模工具，我们可以描述电机的动态行为和控制器的工作原理，从而获得一个完整的控制系统模型。
我们可以通过仿真来评估不同的控制策略，优化系统性能，并进行实验验证。
2. 实例二：PID控制器设计在控制系统中，PID（Proportional-Integral-Derivative）控制器是一种常用的控制器类型。
使用MATLAB中的Control System Toolbox，我们可以设计和调试PID控制器。
我们可以通过输入系统的传递函数或状态空间模型来创建控制系统对象。
可以利用Control System Toolbox提供的自动调整功能，根据系统的要求和性能指标，自动调整PID控制器的参数来实现系统稳定和性能优化。
三、控制系统分析和优化1. 实例三：系统频域分析在设计控制系统时，频域分析是一种重要的方法。
MATLAB提供了许多函数和工具，可用于进行频域分析。
我们可以使用bode函数来绘制系统的频率响应曲线，了解系统的增益和相位裕度，并进行稳定性分析。
MATLAB还提供了工具来进行奈奎斯特图和极点分析等分析方法，帮助用户更好地理解系统的动态特性。
2. 实例四：多目标优化在实际控制系统设计中，通常需要同时满足多个设计指标，例如稳定性、快速响应和抑制干扰等。
matlab数学实验报告
MATLAB数学实验报告指导老师：班级：小组成员：时间：201_/_/_Matlab第二次实验报告小组成员：1题目：实验四;MATLAB选择结构与应用实验目的：掌握if选择结构与程序流程控制;重点掌握break;return;pause语句的应用..问题：问题1：验证“哥德巴赫猜想”;即：任何一个正偶数n>=6均可表示为两个质数的和..要求编制一个函数程序;输入一个正偶数;返回两个质数的和..问题分析：由用户输入一个大于6的偶数;由input语句实现..由if判断语句判断是否输入的数据符合条件..再引用质数判断函数来找出两个质数;再向屏幕输出两个质数即可..编程：function z1;z2=geden;n=input'please input n'if n<6disp'data error';returnendif modn;2==0for i=2:n/2k=0;for j=2:sqrtiif modi;j==0k=k+1;endendfor j=2:sqrtn-iif modn-i;j==0k=k+1;endendif k==0fprintf'two numbers are'fprintf'%.0f;%.0f';i;n-ibreakendendend结果分析：如上图;用户输入了大于6的偶数返回两个质数5和31;通过不断试验;即可验证哥德巴赫猜想..纪录：if判断语句与for循环语句联合嵌套使用可使程序结构更加明晰;更快的解决问题..2题目：实验四;MATLAB选择结构与应用实验目的：用matlab联系生活实际;解决一些生活中常见的实际问题..问题：问题四：在一边长为1的四个顶点上各站有一个人;他们同时开始以等速顺时针沿跑道追逐下一人;在追击过程中;每个人时刻对准目标;试模拟追击路线;并讨论.. （1）四个人能否追到一起（2）若能追到一起;每个人跑过多少路程（3）追到一起所需要的时间设速率为1问题分析：由正方形的几何对称性和四个人运动的对称性可知;只需研究2个人的运动即可解决此问题..编程：hold onaxis0 1 0 1;a=0;0;b=0;1;k=0;dt=0.001;v=1;while k<10000d=norma-b;k=k+1;plota1;a2;'r.';'markersize';15;plotb1;b2;'b.';'markersize';15;fprintf'k=%.0f b%.3f;%.3f a%.3f;%.3f d=%.3f\n';k;b1;b2;a1;a2;da=a+b1-a1/d*dt;b2-a2/d*dt;b=b+b2-a2/d*dt;-b1-a1/d*dt;if d<=0.001breakendendfprintf'每个人所走的路程为：%.3f';k*v*dtfprintf'追到一起所需要的时间为%.3f';k*dt结果分析：上图为2人的模拟运动路线;有对称性可解决所提问题..-上图为运算过程和运算结果..四个人可以追到一起;走过的路程为1.003;时间也为1.003.纪录：此题利用正方形和运动的对称性可以简便运算..3题目：实验八;河流流量估计与数据插值目的：由一些测量数据经过计算处理;解决一些生活实际问题..问题：实验八上机练习题第三题：瑞士地图如图所示;为了算出他的国土面积;做以下测量;由西向东为x轴;由南向北为y轴;从西边界点到东边界点划分为若干区域;测出每个分点的南北边界点y1和y2;得到以下数据mm..已知比例尺1:2222;计算瑞士国土面积;精确值为41288平方公里..测量数据如下：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158 ;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;问题分析：先由题目给定的数据作出瑞士地图的草图;再根据梯形法;使用trapz语句;来估算瑞士国土的面积..编程：x=7.0 10.5 13.0 17.5 34 40.5 44.5 48 56 61 68.5 76.5 80.5 91 96 101 104 106 111.5 118 123.5 136.5 142 146 150 157 158;y1=44 45 47 50 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68;y2=44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121 121 122 116 83 81 82 86 85 68;plotx;y1;'r.';'markersize';15;plotx;y2;'r.';'markersize';15;axis0 160 0 135grid;hold ont=7:158;u1=splinex;y1;t;u2=splinex;y2;t;plott;u1plott;u2s1=trapzt;u1;s2=trapzt;u2;s=s2-s1*2222*22222/10000000;fprintf'S=%.0f';s结果分析：上图为由所给数据绘制出的瑞士地图..上图为运算结果;计算出瑞士的国土面积为42472平方公里;与准确值41288较为接近..纪录：使用梯形分割的方法;trapz语句可以方便计算不规则图形面积;但存在一定误差..4题目：实验七：圆周率的计算与数值积分目的：将数值积分最基本的原理应用于matlab之中;解决一些与积分有关的问题..问题：实验七上机练习题第一题：排洪量某河床的横断面如图7.3所示;为了计算最大排洪量;需要计算其断面积;试根据所给数据m用梯形法计算其断面积..问题分析：河床断面可近似分割成若干曲边梯形;近似处理把它们当做梯形来计算面积可使问题得到简化..编程：clc;clear;x=0 4 10 12 15 22 28 34 40;y=0 1 3 6 8 9 5 3 0;y1=10-y;plotx;y1;'k.';'markersize';15;axis0 40 0 10;grid;hold ont=0:40;u=splinex;y1;t;plott;u;s=40*10-trapzt;u;fprintf's=%.2f\n';s结果分析：上图为河床的断面图..上图为计算结果面积约为180.70平方米..纪录：使用梯形法计算不规则图形面积十分简便易行..5题目：实验七：圆周率的计算与数值积分目的：使用matlab计算解决一些有关积分的问题..问题：实验七上机练习题第三题：从地面发射一枚火箭;在最初100秒内记录其加速度如下;试求火箭在100秒时的速度..Ts=0 10 20 30 40 50 60 70 80 90 100;Am/s*s=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.67 54.01 57.23;问题分析：加速度为速度的微分;已知微分求积分;类似于面积问题;可使用梯形法来计算..编程：clc;clear;x=0 10 20 30 40 50 60 70 80 90 100;y=30.00 31.63 33.44 35.47 37.75 40.33 43.29 46.69 50.6754.01 57.23;plotx;y;'k.';'markersize';15;axis0 100 20 60;grid;hold ons=0:10:100;z=splinex;y;s;plots;y;v=trapzx;y;fprintf'v=%.2f\n';v结果分析：上图为加速度变化图..上图为计算结果;求得火箭在100秒时速度约为4168.95m/s..纪录：梯形法可以推广解决许多已知微分求积分的其他问题..6题目：实验七：圆周率的计算与数值积分目的:计算曲线弧长闭曲线周长可使用微元法;ds=sqrtdx^2+dy^2;在转化微积分问题;累加即可得到结果..问题：实验七上机练习题第三题：计算椭圆想x^2/4+y^2=1的周长;使结果具有五位有效数字..问题分析：编程：s=0;dx=0.001;for x=0:0.001:1.999dy=1.-x+0.001.^2/4-1.-x.^2/4;ds=sqrtdx.^2+dy.^2;s=s+ds;ends=4*s;fprintf'the length is'fprintf'%.4f';s结果分析：上图为计算结果;给定椭圆的周长约为9.1823五位有效数字纪录：计算不规则曲线弧长;可使用微元法;划分为若干小的看做直角三角形;利用勾股定理解决..7题目：实验九人口预测与数据拟合目的：掌握一些曲线拟合的方法;了解曲线拟合常用函数..问题：用电压U=10v的电池给电容器充电;t时刻的电压Vt=U-U-V0exp-t/τ;其中V0是电容器的初始电压;τ是充电常数;由所给数据确定V0和τ..t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;问题分析：题中已给出函数关系式;为指数函数曲线拟合;将所给函数式整理可得标准的exp形函数曲线;从而便于解决..编程：t=0.5 1 2 3 4 5 7 9;V=3.64 3.52 2.74 1.78 1.34 1.01 0.57 0.37;plott;V;'k.';'markersize';20;axis0 10 0 4;grid;hold onpause0.5n=8;a=sumt1:n;b=sumt1:n.*t1:n;c=sumlogV1:n;d=sumt1:n.*logV1:n;A=n a;a b;B=c;d;p=invA*Bx=0:10;y=expp1+p2*x;plotx;y;'r-';'linewidth';2结果分析：上图为电压与时间关系图..上图为计算结果;即U-V0=1.4766;所以V0=8.5234;-1/τ=-0.2835;所以τ=3.5273纪录：曲线拟合的一个重难点是选择合适的曲线函数;才能提高拟合度..8题目：实验七圆周率的计算与数值积分目的：拓展圆周率的各种计算方法;掌握其他数值的近似计算方法..问题：实验七练习2：计算ln2的近似值精确到10的-5次方（1）利用级数展开的方法来计算（2）利用梯形法计算（3）利用抛物线法问题分析：级数展开;梯形法;抛物线法是常见的近似运算方法..编程：1级数展开的方法clc;clear;n=0;r=1;p=0;k=-1;while r>=0.1e-5n=n+1;k=k*-1;p1=p+k/n;r=absp1-p;fprintf'n=%.0f;p=%.10f\n';n;p1;p=p1;end2梯形法clc;clear;f=inline'1./x';x=1:0.1:2;y=fx;p=trapzx;y;fprintf'p=%.6f\n';p3抛物线法clc;clear;f=inline'1./x';a=1;b=2;n=1;z=quadf;a;b;fprintf'z=%.10f\n';z结果分析：（1）级数展开的方法（2）梯形法3抛物线法纪录：级数展开法;梯形法;抛物线法;计算近似值时应合理利用..梯形法和抛物线法不易提高精确度;级数展开法可以提高精确度..9题目：实验八河流流量估计与数据插值目的：掌握求插值多项式的方法;并利用此计算近似值..问题：已知y=fx的函数表如下x=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382;求四次拉格朗日插值多项式;并由此求f0.596问题分析：利用所给函数表可计算拉格朗日插值多项式..编程：function p=lagrangex;yL=lengthx;a=onesL;for j=2:La:;j=a:;j-1.*x';endx=inva*y';for i=1:Lpi=xL-i+1;endx=0.40 0.55 0.65 0.80 0.90 1.05;y=0.41075 0.57815 0.69675 0.88811 1.02652 1.25382; plotx;y;'k.';'markersize';15axis0 2 0 2grid;hold on;p=lagrangex;y;t=0:0.1:1.5;u=polyvalp;t;plott;u;'r-'a=polyvalp;0.596结果分析：上图为所求结果;估算值和插值多项式..纪录：插值多项式是一项十分实用的方法..10题目：求正整数n的阶乘：p=1*2*3*…*n=n;并求出n=20时的结果目的：练习使用循环变量解决数学问题问题：对程序：Clear;clc;n=20;p=1;for i=1:np=p*i;fprintf’i=%.0f;p=%.0f\n’;i;pend进行修改使它：利用input命令对n惊醒赋值问题分析：题中给出程序中“n=20”修改;使用input命令；讲题中的输出命令放出循环之外..编程：clear;clc;n=input'n=';p=1;for i=1:np=p*i;endfprintf'i=%.0f;p=%.0f\n';i;p结果：n=20i=20;p=2432902008176640000>>结果分析：使用input命令可以实现人机对话;使用户自由赋值；输出语句在程序中的位置对输出的结果有很大的影响；在循环内部可以在计算过城中不断输出结果；在循环之外则可以控制只输出最后结果..11题目：对于数列｛√2｝;n=1;2;…;求当其前n项和不超过1000时的n的值及合的大小..目的：运用条件循环解决文帝个项数的循环程序求解;问题：对程序：clear;clc;n=0;s=0;while s<=1000n=n+1;s=s+sqrtn;fprintf’n=%.0f;s=%.4f\n’;n;send问题分析：题中所给程序中的限制变量为上次循环之后的s;导致s超过上限后仍有一次的循环;若把循环变量改为这次的s;则可以避免这种情况的发生..编程：clear;clc;n=0;s=0;while s+sqrtn<=1000n=n+1;s=s+sqrtnfprintf'n=%.0f;s=%.4f\n';n;send结果：……s =970.8891n=128;s=970.8891s =982.2469n=129;s=982.2469s =993.6487n=130;s=993.6487>>结果分析：从结果中可以看出;最后一步为我们需要的答案;从这道题我们可以得出循环变量对一道编程的重要性..。
软件技术(MATLAB)实验手册
实验课需要我们干什么？一、实验课前的知识点预习。
二、准备实验课要演练的例题。
三、实验课后的总结。
MA TLAB实验一MA TLAB软件环境一、实验目标1.熟悉MA TLAB软件的运行环境以及界面构成。
2.熟悉MA TLAB环境中命令窗口的基本操作。
3.熟悉变量的相关操作。
4.熟悉MA TLAB中基本运算函数。
二、实验内容及要求1.运行MA TLAB软件，观察MA TLAB的界面组成。
并了解各个部分的功能。
（标题栏+菜单栏+工具栏+命令窗口+历史命令窗口+当前目录窗口+工作空间）。
2.命令窗口主要用于进行MA TLAB的各种命令操作。
在该窗口中，可以运行MA TLAB的指令、函数、表达式以及显示除了图形外的所有运行结果。
3.变量的赋值以及引用。
4.MA TLAB基本运算函数的演练。
三、实验步骤（略）四、训练习题1.书上所有例题和课后习题。
2.ppt上的学生演练习题。
五、测试习题1.x=2;MA TLAB实验二深入学习MA TLAB环境设置一、实验目标1．熟悉MA TLAB软件环境设置（包括搜索路径和工作空间等设置）。
2．熟悉变量的存储与读取操作。
3．熟悉.m文件。
二、实验内容及要求1．查找和设置搜索路径，并且测试。
2．在工作空间中对变量进行各种操作。
3．使用save和load命令存储与读取变量。
4．区分.m文件中的命令文件和函数文件。
三、实验步骤（略）四、训练习题1．如何设置和更改matlab的初始化当前目录？2．演练课本中对应章节的例子。
3．独立完成课堂布置的实验作业：add.m 和diaoyongjiafa.m文件。
4．测试并体验.m文件的两大分类的区别和相似。
发现并解决.m文件中的若干问题（函数文件中是否可以定义多个函数？主函数子函数可以互相调用么？外部函数可以调用子函数么？子函数可以调用其他函数么？）。
五、测试习题1．精度计算的问题。
看如下的代码，分析：x=0;for n=1:82x=x+0.01endx==0.82MA TLAB实验三MA TLAB编程基础一、实验目标2．熟悉MA TLAB三大程序结构。
MATLAB实验报告
MATLAB实验报告MATLAB实验报告姓名：学号：专业班级：指导⽼师：地点：⽇期：2012.12.17⽬录实验⼋............................................... - 3 -⼀、实验内容...................................... - 3 -⼆．实验结果：.................................... - 3 - 三．实验代码：.................................... - 4 - 实验九............................................... - 5 - 实验⼗⼀............................................. - 7 - ⼀．实验内容...................................... - 7 - ⼆．实验结果...................................... - 7 - 三．实验代码...................................... - 7 - 实验⼗⼆............................................. - 7 - ⼀．实验内容...................................... - 7 - ⼆．实验结果...................................... - 8 - 三．实验代码：................................... - 11 -实验⼋⼀、实验内容2、将100个学⽣5门功课的成绩存⼊矩阵P中，进⾏如下处理：（1）分别求每门课的最⾼分、最低分及相应学⽣序号。
（2）分别求没门功课的平均分和标准⽅差。
MATLAB)课后实验答案[1](最新整理)
](https://img.taocdn.com/s3/m/2bd1298ccf84b9d529ea7a94.png)
5. 下面是一个线性方程组：
1112340.95x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥
实验三 选择结构程序设计
1. 求分段函数的值。
222603560523
1x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩
且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
解：M 文件如下：
2. 输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。
其中90分~100分为A，80分~89分为B，79分~79分为C，60分~69分为D，60分以下为E。
要求：
(1) 分别用if语句和switch语句实现。
(2) 输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
解：M文件如下
3. 硅谷公司员工的工资计算方法如下：
(1) 工作时数超过120小时者，超过部分加发15%。
(2) 工作时数低于60小时者，扣发700元。
(3) 其余按每小时84元计发。
试编程按输入的工号和该号员工的工时数，计算应发工资。
解：M文件下
运算结果如下；
5. 若两个连续自然数的乘积减1是素数，则称这两个边疆自然数是亲密数对，该素数是亲密素数。
例如，2×3-1=5，由于5是素数，所以2和3是亲密数，5是亲密素数。
求[2,50]区间内：
(1) 亲密数对的对数。
(2) 与上述亲密数对对应的所有亲密素数之和。
解：
M文件：。
(完整版)MATLAB)课后实验答案[1]
![(完整版)MATLAB)课后实验答案[1]](https://img.taocdn.com/s3/m/d27ebaa958fb770bf78a55b8.png)
实验一 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin 851z e =+(2) 21ln(2z x =，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e a z a a --+=++=--L (4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.5 解：4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
解：(1) 结果：(2). 建立一个字符串向量例如：ch='ABC123d4e56Fg9';则要求结果是：实验二 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A O S +⎡⎤=⎢⎥⎣⎦。
解: M 文件如下；5. 下面是一个线性方程组：1231112340.951110.673450.52111456x x x ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎣⎦(1) 求方程的解。
(2) 将方程右边向量元素b 3改为0.53再求解，并比较b 3的变化和解的相对变化。
(3) 计算系数矩阵A 的条件数并分析结论。
解： M 文件如下： 123d4e56g9实验三 选择结构程序设计1. 求分段函数的值。
2226035605231x x x x y x x x x x x x ⎧+-<≠-⎪=-+≤<≠≠⎨⎪--⎩且且及其他用if 语句实现，分别输出x=-5.0,-3.0,1.0,2.0,2.5,3.0,5.0时的y 值。
MATLAB实验讲义
MATLAB实验讲义目录实验大纲 (2)实验一/二 MATLAB的基础操作 (3)实验三 MATLAB运算基础(一) (3)实验四 MATLAB运算基础(二) (4)实验五循环结构程序设计（一） (5)实验六循环结构程序设计（二） (5)实验七 MATLAB的绘图操作(一) (6)实验八 MATLAB的绘图操作(二) (7)实验九函数和文件（一） (7)实验十函数和文件（二） (7)实验十一线性代数中的数值计算问题 (8)实验十二 MATLAB函数库的运用（一） (9)实验十三 MATLAB函数库的运用（二） (10)《MATLAB》课程实验教学大纲课程名称：MATLAB（MATLAB）课程编号：16072327课程性质：选修实验总学时：27实验室名称：电子设计自动化一、课程简介：本课程是电气工程及其自动化、自动化、电力工程与管理专业本科生的学科基础选修课，它在线性代数、信号分析和处理、控制系统设计和仿真等方面有着广泛的应用。
主要是学习MATLAB的语法规则、基本命令和使用环境，使学生掌握MATLAB的基本命令和基本程序设计方法，提高使用该语言的应用能力，具有使用MATLAB语言编程和调试的能力，以便为后续多门课程使用该语言奠定必要的基础。
二、课程实验目的与要求：1.基本掌握MATLAB在线帮助功能的使用、熟悉MATLAB运行环境和MATLAB语言的主要特点，掌握MATLAB语言的基本语法规则及基本操作命令的使用，学会M文件的建立和使用方法以及应用MATLAB实现二维和三维图形的绘制方法，具有使用MATLAB语言编程和调试的能力。
2.初步掌握MATLAB在电路和信号与系统中的应用。
3.能根据需要选学参考书，查阅手册，通过独立思考，深入钻研有关问题，学会自己独立分析问题、解决问题，具有一定的创新能力。
三、主要仪器设备及台（套）数：计算机50台、MATLAB软件五、主要参考书目：1.《MATLAB及在电子信息课程中的应用》陈怀琛、杨吉斌编著，电子工业出版社，2002年1版2.《MATLAB7.0编程基础》王家文、王皓、刘海等；机械工业出版社，2005年7月3.《MATLAB教程——基于6.x版本》张志涌、徐彦琴等；北京航空航天大学出版，2001年4月出版实验一/二 MATLAB的基础操作一、实验目的1、掌握MATLAB的启动和退出。
Matlab数学实验报告
实验一 Matlab基本操作1．实验课程名称数学实验2．实验项目名称Matlab基本操作3．实验目的和要求了解Matlab的基本知识，熟悉其上机环境，掌握利用Matlab进行基本运算的方法。
4．实验内容和原理内容：三角形的面积的海伦公式为：area=)s-sa--)()(s(csb其中： s=(a+b+c)/2原理：将一般数学问题转化成对应的计算机模型并进行处理的能力。
了解Matlab的基本功能，会进行简单的操作。
5．主要仪器设备计算机与Windows 2000/XP系统；Matlab等软件。
6．操作方法与实验步骤步骤：（1）在M文件编辑窗口输入以下程序，并以文件名”area_helen.m”保存：a= input(‘a=‘) ； b= input(‘b=‘) ； c= input(‘c=‘) ；s= (a+b+c)/2;area=sqrt (s* (s-a) * (s-b) * (s-c))（2）在命令窗口输入文件名“area_helen”，按回车键，即可运行上面的程序，输入三边长，立即可得三角形面积（3）第二题在命令窗口输入b=6;a=3;c=a*b,d=c-2*b(4) 按回车键，即可运行上面的程序7．实验结果与分析<1> a=3; b=4; c=5;时，aera=6 当a为3，b为4，c为5时，s=6，aera=6<2> c= 18,d=6，a为3，b为6时，c=18，d=6实验二 Matlab的数值计算1．实验课程名称数学实验2．实验项目名称Matlab的数值计算3．实验目的和要求了解一些简单的矩阵、向量、数组和多项式的构造和运算方法实例，懂得编写简单的数值计算的Matlab程序。
熟悉一些Matlab的简单程序，会用Matlab的工具箱，懂得Matlab的安装和简单的使用。
4．实验内容和原理内容：从函数表：)1(),5.0(),2( ,0x 1x 021x 1x f(x) 32-⎪⎩⎪⎨⎧≤≤<>+=f f f x x求设)1(),2( ,1211)(2-⎩⎨⎧≤>+=f f x xx x x f 求设 原理：利用矩阵、向量、数组、和多项式的构造和运算方法，用常用的几种函数进行一般的数值问题求解。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、实验目的
彩色图像处理
使用MatLab 软件对图像进行彩色处理。使学生通过实验熟悉使用MatLab软件进行图像 彩色处理的有关方法，并体会到图像彩色处理技术以及对图像处理的效果。
二、实验要求
要求学生能够完成彩色图像的分析，能正确讨论彩色图像的亮度、色调等性质；会对彩 色图像进行直方图均衡， 并能正确解释均衡处理后的结果； 能够对单色图像进行伪彩色处理、 利用多波长图像进行假彩色合成、进行单色图像的彩色变换。
T＝91，T＝130 ，T＝43，比较效果，见教材 P92 3.得到图像 I 的大小用函数： [M,N]=size(I)，M 为行，N 为列数 I=imread('blood1.bmp');figure; subplot(221);imshow(I);title('原图形'); %subplot(122);imhist(I);title('直方图'); I=double(I);I1=double(I);I2=double(I);I3=double(I); [M,N]=size(I); T=91; for i=1:M for j=1:N if I(i,j)<T I1(i,j)=0; else I1(i,j)=255; end end end subplot(222);imshow(uint8(I1));title('T=91'); T=130; for i=1:M for j=1:N if I(i,j)<T I2(i,j)=0; else I2(i,j)=255;
的变换函数为-sin( )，绿色分量G的变换函数为-cos( );，蓝色分量B的变换函数为sin( )； 显示变换曲线及变换合成的彩色图像并观察彩色变换图像的色调与原始图像灰度之间 的关系； 将RGB的变换公式至少互换一次(例如R与G互换)，显示变换曲线、变换结果并观察处理 的结果。 (6)打印全部结果并进行讨论。 利用MATLAB软件实现彩色图像处理的程序：
原图形 T=91
T=130
T=43
（二） 采用迭代式阈值选择方法对图像 blood.bmp 进行二值化处理， 并和(一)中处理的结果 比较。 f=imread('blood1.bmp'); subplot (121); imshow(f); title('原始图像'); f=double(f); T=(min(f(:))+max(f(:)))/2; done=false; i=0 while ~done R1=find(f<=T); R2=find(f>T); Tnew=(mean(f(R1))+mean(f(R2)))/2 %求 2 个灰度均值的平均值 done=abs(Tnew-T)<1 % abs 为求绝对值
二、实验要求
要求学生能够自行评价各主要算子在无噪声条件下和噪声条件下的分割性能。能够 掌握分割条件(阈值等)的选择。 完成规定图像的处理并要求正确评价处理结果， 能够从理论 上作出合理的解释。
三、实验内容与步骤
MATLAB 的图像处理工具箱中提供的 edge 函数可以实现检测边缘的功能，其语法格式 如下： BW = edge(I, ‘sobel’) BW = edge(I, ‘sobel’, direction) BW = edge(I, ‘roberts’) BW = edge(I, ‘log’) 这里 BW = edge(I, ‘sobel’)采用 Sobel 算子进行边缘检测。BW = edge(I, ‘sobel’, direction) 可以指定算子方向，即： direction = ‘horizontal’，为水平方向； direction = ‘vertical’，为垂直方向； direction = ‘both’，为水平和垂直两个方向。 BW = edge(I, ‘roberts’)和 BW = edge(I, ‘log’)分别为用 Roberts 算子和拉普拉斯高斯算子 进行边缘检测。 (1)使用 Roberts 算子的图像分割实验 调入并显示图像 room.tif 中图像；使用 Roberts 算子对图像进行边缘检测处理； Roberts 算子为一对模板：
T=Tnew; i=i+1; end f(R1)=0; f(R2)=1; subplot(122);imshow(f); title('迭代阈值二值化图像');
原始图像 迭代阈值二值化图像
（三）Otsu 法阈值分割 I = imread('lajiao.bmp'); level = graythresh(I); BW = im2bw(I,level); figure, imshow(BW)
⎧0 ⎪ g ( x, y ) = ⎨ ⎪ ⎩255 阈值可取为
f ( x, y ) < T
f ( x, y ) ≥ T
end end end subplot(223);imshow(uint8(I2));title('T=130'); T=43; for i=1:M for j=1:N if I(i,j)<T I3(i,j)=0; else I3(i,j)=255; end end end subplot(224);imshow(uint8(I3));title('T=43');
相应的矩阵为：rh = [0 1；-1 0]； rv = [1 0；0 -1]；这里的 rh 为水平 Roberts 算子，rv 为垂直 Roberts 算子。分别显示处理后的水平边界和垂直边界检测结果；用“欧 几里德距离”和“街区距离”方式计算梯度的模， 并显示检测结果； 对于检测结果进行二值 化处理，并显示处理结果； 提示：先做检测结果的直方图，参考直方图中灰度的分布尝试确定阈值；应反复调节
七、实验图片
实验十
图像的阈值分割
一、实验目的：
1、进一步了解 MATLAB 关于图像处理的相关指令。 2、掌握图像阈值分割的基本原理和基本方法。 3、能够用所学方法实现图像分割。
二、实验内容
(一) 采用人工选择阈值法编写程序，把图像 blood.bmp 分割为二值图像。尝试不同阈值来 分割图像，比较分割效果，选出最满意的阈值。 （提示： 1.先做图像的直方图，参考直方图中灰度的分布尝试确定阈值；应反 复调节阈值的大小，直至二值化的效果最为满意为止。 2.运算之前要先转化为双精度：I＝double(I)，运算之后再转化为二 进制数表示：I＝uint8(I) 变换式如下：
三、实验内容与步骤
(1) 彩色图像的分析 调入并显示彩色图像flower1.tif ； 拆分这幅图像，并分别显示其R，G，B分量； 根据各个分量图像的情况讨论该彩色图像的亮度、色调等性质。 (2) 彩色图像的直方图均衡 接内容(1)； 显示这幅图像的R，G，B分量的直方图，分别进行直方图均衡处理，并显示均衡后的直 方图和直方图均衡处理后的各分量； 将处理完毕的各个分量合成彩色图像并显示其结果； 观察处理前后图像的彩色、亮度、色调等性质的变化。 (3) 假彩色处理 调入并显示红色可见光的灰度图像vl_red.jpg、绿色可见光的灰度图像vl_green.jpg 和蓝色可见光的灰度图像vl_blue.jpg；以及近红外灰度图像infer_near.jpg和中红外灰度 图像infer_mid.jpg； 以图像vl_red.jpg为R；图像vl_green.jpg为G；图像vl_blue.jpg为B，将这三幅图像组 合成可见光RGB彩色图像； 分别以近红外图像infer_near.jpg和中红外图像infer_mid替换R分量，形成假彩色图 像； 观察处理的结果，注意不同波长红外线图像组成图像的不同结果 (4) 伪彩色处理1：灰度切片处理 调入并显示灰度图像head.jpg； 利用MATLAB提供的函数对图像在8~256级的范围内进行切片处理， 并使用hot模式和cool 模式进行彩色化； 观察处理的结果。 (5) 彩色变换(选做) 调入并显示灰度图像Lenna.jpg； 使用不同相位的正弦函数作为变换函数，将灰度图像变换为RGB图像。其中红色分量R
figure,imshow(cut_2),colormap(cool) %显示图ห้องสมุดไป่ตู้ cut_2,并使用 cool 模型彩色化
ice(交互彩色编辑)函数的参数参照书后面的注释。
四、实验仪器与软件
1．计算机； 2．MATLAB 程序； 3．移动式存储器（软盘、U 盘等） 。 4．记录用的笔、纸。
五、实验报告要求
三、实验仪器与软件
1．计算机； 2．MATLAB 程序； 3．移动式存储器（软盘、U 盘等） 。 4．记录用的笔、纸。
四、实验报告要求
1．叙述实验过程； 2．提交实验的原始图像和结果图像。
实验十一
一、实验目的
图像分割（边缘检测）
使用MatLab 软件进行图像的分割。使学生通过实验体会一些主要的分割算子对图像处 理的效果，以及各种因素对分割效果的影响。
rgb_image=imread('flower1.tif'); fR=rgb_image(:,:,1); fG=rgb_image(:,:,2); fB=rgb_image(:,:,3); figure(1),imshow(fR) figure(2),imshow(fG) figure(3),imshow(fB) %实现 rgb 图像转化为 NTSC 彩色空间的图像 yiq_image=rgb2ntsc(rgb_image); fY=yiq_image(:,:,1); fI=yiq_image(:,:,2); fQ=yiq_image(:,:,3); figure(4),imshow(fY) figure(5),imshow(fI) figure(6),imshow(fQ) fR=histeq(fR,256); fG=histeq(fG,256); fB=histeq(fB,256); RGB_image=cat(3,fR,fG,fB); %将直方图均衡化后的彩色图像合并 figure,imshow(RGB_image) f1=imread('v1_red.jpg'); f2=imread('v1_green.jpg'); f3=imread('v1_blue.jpg'); f4=imread('infer_near.jpg'); ture_color=cat(3,f1,f2,f3); figure,imshow(ture_color) false_color=cat(3,f4,f2,f3); figure,imshow(false_color) f=imread('head.jpg'); cut_1=imadjust(f,[0.0925 0.5],[0.0925 0.5]);%提取灰度在 16-128 之间的像素 cut_2=imadjust(f,[0.5 1],[0.5 1]); figure,imshow(cut_1),colormap(hot) (选做) f=imread('Lenna.jpg'); g=ice('image',f); %通过 ice(交互彩色编辑)函数对图像进行变换 %提取灰度在 128-256 之间的像素 %显示图像 cut_1,并使用 hot 模型彩色化 %显示由红、绿、蓝三幅图合成的彩色图 %用近红外图像代替 R 分量 %显示由近红外、绿、蓝三幅图合成的假彩色图 %观察处理后的彩色图色度，亮度参照前面 %对彩色图像的分量进行直方图均衡化 %图像 flower1.tif 的亮度 %图像 flower1.tif 的色调 %图像 flower1.tif 的饱和度 %读取图像 flower1.tif %获取图像的红色分量 %获取图像的绿色分量 %获取图像的蓝色分量 %分别显示图像