求助MATLAB，以下是详细问题

Question

测得平板表面x=1:5,y=1:3网格点处的温度分别为：
82 81 80 82 84
79 63 61 65 81
84 84 82 85 86
试作出平板表面的温度分布曲面z=f(x,y)在网格x=1:0.1:5,y=1:0.1:3上的图形。

Answer 1

这是对5x3数据直接绘制得到的图，

还可以用插值的方法去作更为平滑的图。matlab下的命令为：

tmp = [82  81  80  82 84;

79  63  61  65  81;

84  84  82  85  86;]

figure

meshc(1:5,1:3,tmp);

[xi,yi] = meshgrid(1:0.1:5,1:0.1:3)

[xi,yi] = meshgrid(1:0.1:5,1:0.1:3);

zi3 = interp2(x,y,tmp,xi,yi,'bicubic');

figure

meshc(1:0.1:5,1:0.1:3,zi3)

figure

surf(xi,yi,zi3);

Answer 2

用法：
morsecode=morse('matlabsky baidu wacs5',0) %编码
srccode=morse(morsecode,1) %解码

morsecode(2)='.'; %修改码
srccode=morse(morsecode,1) %错误的morse码进行解码

=============以下是m文件==========
function Dst=morse(Src,Flag)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%Author : wacs5
%%Email : wacs5@126.com
%%Date : 20101031
%%Function : morse(莫斯密码)编码和解码
%%Argument : Src为待编码/解码的字符串
%% : Flag=0为编码，Flag=1为解码
%%Return : Dst为完成的编码/解码的字符串
%%Called Mtb : disp,find,isempty,isequal,length,nargin,strrep,warning
%%Example:
%% : morsecode=morse('matlabsky baidu wacs5',0)
%% : srccode=morse(morsecode,1)
%% : morsecode(2)='.'; srccode=morse(morsecode,1) %错误的morse码
%%See also :
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if nargin==1
Flag=0; %缺省为编码
end
Dst=[];

SrcAlpha=['A':'Z','1':'9','0']; %原始字母
EncAlpha=; %相应的morse码

if Flag==0
disp('Encoding'); %编码
for i=1:length(Src)
if Src(i)>='a' && Src(i)<='z' %小写转大写
Src(i)=Src(i)-32;
end
index=find(SrcAlpha==Src(i)); %查找待编码字母在SrcAlpha中的位置
if isempty(index) %不在字典中
Dst=[Dst,' ']; %两个空格。因为字母后面已经有一空格
else
Dst=[Dst,EncAlpha,' ']; %每个morse码后加一空格
end
end
else
disp('Decoding'); %解码
if Src(end)~=' ' %最后一个不是空格，添加一个空格，仅为了方便编程
Src(end+1)=' ';
end
SpaceIndex=find(Src==' '); %找出所有空格位置
SpaceIndex=[0,SpaceIndex]; %将0加入，也仅仅是为了方便编程
for i=1:length(SpaceIndex)-1
Code=Src(SpaceIndex(i)+1:SpaceIndex(i+1)-1);%获得morse密文
if isempty(Code) %连续两个空格
Dst=[Dst,' ']; %译文也为空格
else
found=0; %是否在morse编码字典中存在
for j=1:length(EncAlpha)
if (isequal(EncAlpha,Code)) %在morse编码字典中找到
found=1;
break;
end
end
if found==0 %morse编码字典中没有找到，有误
Dst=[Dst,'*'];
warning(['Input Code frag "',Code,'" wrong']);
else
Dst=[Dst,SrcAlpha(j)]; %相应的原文
end
end
end
Dst=strrep(Dst,' ',' '); %多空格合为一个空格
end