是什么原因啊,想画一个cos函数的三维图,然后能一帧一帧的沿X轴正向移动，请帮帮忙

Question

clear; %clear all variables from memory
close all; %close all windows
clc; %clear command window

disp('Wave Animation 1'); %display the title

%INPUTS

%wave parameters
frequency = 1e6;
velocity = 3e3;

%x parameters
x_end = 3e-02;%TO BE ENTERED%
x_step = 1.5e-04;%TO BE ENTERED%

%time parameters
t_step = 5e-08;%TO BE ENTERED%
t_points = 100;%TO BE ENTERED%

%PROGRAM

x = [0 : x_step : x_end]; %make a vector of x positions
t = [0 : t_step : t_step * t_points]; %make a vector of times

%calculate w and k from inputs
w = sym('pi*2e6');%TO BE ENTERED%
k = sym('(2/3)*pi*1000');%TO BE ENTERED%

figure(2);
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;

Answer 1

这个程序里主要有2个地方有误。（1）不应该使用符号函数。原先程序中，%calculate w and k from inputsw = sym('pi*2e6');%TO BE ENTERED%k = sym('(2/3)*pi*1000');%TO BE ENTERED%中使用了sym符号函数，但是后面的程序将这两个变量是当做数值计算的，因此不用符号化，改为%calculate w and k from inputsw = 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)));最后的结果为：