UC知道MATLAB 图象具体数据的读取
来源：百度知道 编辑：UC知道 时间：2024/06/08 07:33:57
具体的MATLAB程序
clear
us=5 %cm-1
ua=0.2 %cm-1
k=sqrt(3.*ua.*(ua+us))
d=10 %cm
x=-0.2*d: 0.2:1.2*d
z=0:0.2:d
[X,Z]= meshgrid(x,z) %生成二维矩阵
P=Z.^2.*exp(-k.*(sqrt(X.^2+Z.^2)+sqrt((d-X).^2+Z.^2))).*(k.*sqrt(X.^2+Z.^2)+1).*(k.*sqrt((d-X).^2+Z.^2)+1)./[(X.^2+Z.^2).^(3/2).*((d-X).^2+ Z.^2).^(3/2)]
%pcolor(X,Z,P) %函数pcolor是画二维平面像素图的函数
[C,h] =contour(X,Z,P,60) %40为线条数, 画X轴正半轴的二维等高图
clabel(C,h)%标注轮廓线，l轮廓线示数随线条数变化。
xlabel('x/cm') %标注X轴
ylabel('z/cm') %标注Z轴
title('us=5cm') %标注标题
shading flat %设置着色模式
想要具体得到在X=5的时候对应图象Z的值,哪位达人可以给个答复,谢谢了
哦，谢谢您.不过我想知道的是图做出来后,图上坐标x=5的时候,对应图上的Z的坐标.
clear
us=5 %cm-1
ua=0.2 %cm-1
k=sqrt(3.*ua.*(ua+us))
d=10 %cm
x=-0.2*d: 0.2:1.2*d
z=0:0.2:d
[X,Z]= meshgrid(x,z) %生成二维矩阵
P=Z.^2.*exp(-k.*(sqrt(X.^2+Z.^2)+sqrt((d-X).^2+Z.^2))).*(k.*sqrt(X.^2+Z.^2)+1).*(k.*sqrt((d-X).^2+Z.^2)+1)./[(X.^2+Z.^2).^(3/2).*((d-X).^2+ Z.^2).^(3/2)]
%pcolor(X,Z,P) %函数pcolor是画二维平面像素图的函数
[C,h] =contour(X,Z,P,60) %40为线条数, 画X轴正半轴的二维等高图
clabel(C,h)%标注轮廓线，l轮廓线示数随线条数变化。
xlabel('x/cm') %标注X轴
ylabel('z/cm') %标注Z轴
title('us=5cm') %标注标题
shading flat %设置着色模式
想要具体得到在X=5的时候对应图象Z的值,哪位达人可以给个答复,谢谢了
哦，谢谢您.不过我想知道的是图做出来后,图上坐标x=5的时候,对应图上的Z的坐标.
最末加几行就行。
clc;clear
us=5 %cm-1
ua=0.2 %cm-1
k=sqrt(3.*ua.*(ua+us))
d=10 %cm
x=-0.2*d: 0.2:1.2*d;
z=0:0.2:d;
[X,Z]= meshgrid(x,z); %生成二维矩阵
P=Z.^2.*exp(-k.*(sqrt(X.^2+Z.^2)+sqrt((d-X).^2+Z.^2))).*(k.*sqrt(X.^2+Z.^2)+1).*(k.*sqrt((d-X).^2+Z.^2)+1)./[(X.^2+Z.^2).^(3/2).*((d-X).^2+ Z.^2).^(3/2)];
%pcolor(X,Z,P) %函数pcolor是画二维平面像素图的函数
[C,h] =contour(X,Z,P,60); %40为线条数, 画X轴正半轴的二维等高图
clabel(C,h)%标注轮廓线，l轮廓线示数随线条数变化。
xlabel('x/cm') %标注X轴
ylabel('z/cm') %标注Z轴
title('us=5cm') %标注标题
shading flat %设置着色模式
figure
surf(X,Z,P)
z5=z(find(x==5))
Z5=Z(find(X==5))
结果：
z5 =
7
Z5 =
0
0.2
0.4
0.6
0.8
1
1.2
1.4
1.6
1.8
2
2.2
2.4