作业帮大侠,感谢您帮我解决了这个程序,还有些疑问,望不吝赐教
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/31 12:44:32
大侠,感谢您帮我解决了这个程序,还有些疑问,望不吝赐教
clear;clc
syms n
T=300;H=10;t=1;
x=0:1:10;
for ii=1:length(x)
w(ii)=symsum(((-1)^(n+1)*4*T/(2*n-1)/pi)*cos((2*n-1)*pi*x(ii)/2/H)*exp(-((2*n+1)/(2*H)*pi)^2*t),n,1,inf);
end
plot(x,double(w))
如果我想把步长改为0.01可以吗,我把x=0:1:10;改为x=0:0.01:10;
可是出错了
还有,如果让t以0.1为步长变化,吧t也作为变量,区间是【0,3】,那么能否在一张图上画出30根曲线
clear;clc
syms n
T=300;H=10;t=1;
x=0:1:10;
for ii=1:length(x)
w(ii)=symsum(((-1)^(n+1)*4*T/(2*n-1)/pi)*cos((2*n-1)*pi*x(ii)/2/H)*exp(-((2*n+1)/(2*H)*pi)^2*t),n,1,inf);
end
plot(x,double(w))
如果我想把步长改为0.01可以吗,我把x=0:1:10;改为x=0:0.01:10;
可是出错了
还有,如果让t以0.1为步长变化,吧t也作为变量,区间是【0,3】,那么能否在一张图上画出30根曲线
可以出图,只是运行时间很长,建议x=0:0.1:10
-----------------------------
clear;clc %我的能出图
syms n
T=300;H=10;
t=0.2:0.4:3; %步长0.1,运行时间会很长
x=0:1:10;
for jj=1:length(t)
for ii=1:length(x)
w(jj,ii)=symsum(((-1)^(n+1)*4*T/(2*n-1)/pi)*cos((2*n-1)*pi*x(ii)/2/H)*exp(-((2*n+1)/(2*H)*pi)^2*t(jj)),n,1,inf);
end
plot(x,double(w(jj,:)))
hold on
end
-----------------------------
clear;clc %我的能出图
syms n
T=300;H=10;
t=0.2:0.4:3; %步长0.1,运行时间会很长
x=0:1:10;
for jj=1:length(t)
for ii=1:length(x)
w(jj,ii)=symsum(((-1)^(n+1)*4*T/(2*n-1)/pi)*cos((2*n-1)*pi*x(ii)/2/H)*exp(-((2*n+1)/(2*H)*pi)^2*t(jj)),n,1,inf);
end
plot(x,double(w(jj,:)))
hold on
end