大侠,感谢您帮我解决了这个程序,还有些疑问,望不吝赐教clear;clcsyms nT=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);endplot(x,dou
来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/27 22:18:14
大侠,感谢您帮我解决了这个程序,还有些疑问,望不吝赐教clear;clcsyms nT=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);endplot(x,dou
大侠,感谢您帮我解决了这个程序,还有些疑问,望不吝赐教
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;clcsyms nT=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);endplot(x,dou
可以出图,只是运行时间很长,建议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