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

Question

clear; %clear all variables from memoryclose all; %close all windowsclc; %clear command windowdisp('Wave Animation 1'); %display the title%INPUTS%wave parametersfrequency = 1e6;velocity = 3e3;%x parametersx_end = 3e-02;%TO BE ENTERED%x_step = 1.5e-04;%TO BE ENTERED%%time parameterst_step = 5e-08;%TO BE ENTERED%t_points = 100;%TO BE ENTERED%%PROGRAMx = [0 : x_step : x_end]; %make a vector of x positionst = [0 : t_step : t_step * t_points]; %make a vector of times%calculate w and k from inputsw = 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)));

最后的结果为：