我有道matlab编程有点问题,哪位高手帮我看看,大恩不言谢!dt=0.001;t=0:dt:0.6;fp=30;R=(1-2*(pi*fp*(k*dt-1/(2*fp))).^2)*exp(-(pi*fp*(k*dt-1/(2*fp))).^2);dx=10;dz=10;t0=0.2;x=-500:dx:500;z=-500:dz:500;m=length(x);n=length(z);l=length(t
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/03 12:46:07
![我有道matlab编程有点问题,哪位高手帮我看看,大恩不言谢!dt=0.001;t=0:dt:0.6;fp=30;R=(1-2*(pi*fp*(k*dt-1/(2*fp))).^2)*exp(-(pi*fp*(k*dt-1/(2*fp))).^2);dx=10;dz=10;t0=0.2;x=-500:dx:500;z=-500:dz:500;m=length(x);n=length(z);l=length(t](/uploads/image/z/13582343-47-3.jpg?t=%E6%88%91%E6%9C%89%E9%81%93matlab%E7%BC%96%E7%A8%8B%E6%9C%89%E7%82%B9%E9%97%AE%E9%A2%98%2C%E5%93%AA%E4%BD%8D%E9%AB%98%E6%89%8B%E5%B8%AE%E6%88%91%E7%9C%8B%E7%9C%8B%2C%E5%A4%A7%E6%81%A9%E4%B8%8D%E8%A8%80%E8%B0%A2%21dt%3D0.001%3Bt%3D0%3Adt%3A0.6%3Bfp%3D30%3BR%3D%281-2%2A%28pi%2Afp%2A%28k%2Adt-1%2F%282%2Afp%29%29%29.%5E2%29%2Aexp%28-%28pi%2Afp%2A%28k%2Adt-1%2F%282%2Afp%29%29%29.%5E2%29%3Bdx%3D10%3Bdz%3D10%3Bt0%3D0.2%3Bx%3D-500%3Adx%3A500%3Bz%3D-500%3Adz%3A500%3Bm%3Dlength%28x%29%3Bn%3Dlength%28z%29%3Bl%3Dlength%28t)
我有道matlab编程有点问题,哪位高手帮我看看,大恩不言谢!dt=0.001;t=0:dt:0.6;fp=30;R=(1-2*(pi*fp*(k*dt-1/(2*fp))).^2)*exp(-(pi*fp*(k*dt-1/(2*fp))).^2);dx=10;dz=10;t0=0.2;x=-500:dx:500;z=-500:dz:500;m=length(x);n=length(z);l=length(t
我有道matlab编程有点问题,哪位高手帮我看看,大恩不言谢!
dt=0.001;
t=0:dt:0.6;
fp=30;
R=(1-2*(pi*fp*(k*dt-1/(2*fp))).^2)*exp(-(pi*fp*(k*dt-1/(2*fp))).^2);
dx=10;dz=10;
t0=0.2;
x=-500:dx:500;
z=-500:dz:500;
m=length(x);
n=length(z);
l=length(t);
v=1500;
a=v*dt/dx;
b=v*dt/dz;
for i=1:m;
for j=1:n;
if i==51,j==51
for k=0:round(1/fp/dt);
t=k*dt;
U(i,j,k)=R;
end
else
for k=0:1;
U(i,j,k)=0;
end
end
end
end
for k=0:l-1;
for i=2:m-1;
for j=2:n-1;
U(i,j,k+1)=2*(1-a.^2-b.^2)*U(i,j,k)+a.^2*(U(i+1,j,k)+U(i-1,j,k))+b.^2*(U(i,j+1,k)+U(i,j-1,k))-U(i,j,k-1);
end
end
imshow(U)
? Attempted to access (1,1,0); index must be a positive integer or logical.
我有道matlab编程有点问题,哪位高手帮我看看,大恩不言谢!dt=0.001;t=0:dt:0.6;fp=30;R=(1-2*(pi*fp*(k*dt-1/(2*fp))).^2)*exp(-(pi*fp*(k*dt-1/(2*fp))).^2);dx=10;dz=10;t0=0.2;x=-500:dx:500;z=-500:dz:500;m=length(x);n=length(z);l=length(t
注意,你这块有问题
for k=0:round(1/fp/dt);
t=k*dt;
U(i,j,k)=R;
Matlab中数组的下标是从1开始的,因此不能让k从0开始!