这个程序里主要有2个地方有误。
(1)不应该使用符号函数。
原先程序中,
%calculate w and k from inputs
w = sym('pi*2e6');%TO BE ENTERED%
k = sym('(2/3)*pi*1000');%TO BE ENTERED%
中使用了sym符号函数,但是后面的程序将这两个变量是当做数值计算的,因此不用符号化,改为
%calculate w and k from inputs
w = pi*2e6;%TO BE ENTERED%
k = (2/3)*pi*1000;%TO BE ENTERED%
(2)最后的for循环中,cos函数的计算有问题,
原先程序
for jj=1:length(t)
clf;
m=zeros(1,201);
[mx,mm]=meshgrid(x,m);
mmx=sqrt(mx.^2+mm.^2);
surf(mx,mm,cos(k*mmx-w*jj(t)));
pause(0.01);
end;
中,surf的这条语句
surf(mx,mm,cos(k*mmx-w*jj(t)));
中jj(t)根本就不对,jj是一个数字,怎么可能后面有下标,且t本身就不是一个整型,而是浮点型,
改为
surf(mx,mm,cos(k*mmx-w*t(jj)));
最后的结果为: