作业帮 > 综合 > 作业

使用matlab画图 y(x,r)=∑∑sin(ax)cos(br)*exp(-(a^2+b^2)) a=[2,5,7,

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/17 09:24:39
使用matlab画图 y(x,r)=∑∑sin(ax)cos(br)*exp(-(a^2+b^2)) a=[2,5,7,12,23,34,45] b=[1,5,7,9,14,16,17]
这个问题解不出的话 求一个相似的问题也行
方程中a,b为数组的元素
[x,r]=meshgrid(linspace(-2*pi,2*pi));
y=x-x;
for a=[2,5,7,12,23,34,45]
for b=[1,5,7,9,14,16,17]
y=y+sin(a*x).*cos(b*r)*exp(-(a^2+b^2));
end
end
mesh(x,r,y)
xlabel('x');ylabel('r');zlabel('y')
再问: 我的方程里面存在贝塞尔函数 besselj(0,br)能不能直接加进去啊??y=x-x是啥意思 不好意思啊 大哥 我是小白
再答: 可以的,注意写成b*r即可,还有如果是乘除或次方的话不要忘记加点即可。 x-x就是0了呀,生成一个0矩阵,和zeros(size(x))一样
再问: 如果a=[0,2,5~~~] 并且 y(x,r)=∑∑sin(ax)/a*cos(br)*exp(-(a^2+b^2)) b=0时候 sin(ax)/a就为1了 这样的话 程序能运行出来么 ,分母为0会报错,怎么改程序~~~~谢谢你了 我给你加分
再答: sin(ax)/a当a趋于0时应该等于x而不是1,你把a=0改成a=eps好了,这样结果就差不多了a=[eps,2,5~~~]
再问: 恩 真好 还有个问题如果我的a值不是有限个数,a为acota=10的解,b为bcotb=20的解 那怎么表示啊??大哥,能计算么
再答: 级数除非能找到收敛区间内的表达式否则不可能用循环语句一个个加上去的。
再问: 明白了 辛苦了 分数加没了!!~~