问题我已经贴到百度空间了,请大神帮我看看,来日一定报答
来源:学生作业帮助网 编辑:作业帮 时间:2024/09/29 23:34:21
问题我已经贴到百度空间了,请大神帮我看看,来日一定报答
问题我已经贴到百度空间了,请大神帮我看看,来日一定报答
问题我已经贴到百度空间了,请大神帮我看看,来日一定报答
楼上的程序写的还不错,但是,建议在提交答案之前至少应该做一些基本的验证工作.毕竟,提供一个错误的答案比起没有答案不见得会好到哪里去.
我对这个问题的直观反应是,两个曲面的交线应该会比较平滑,而不应该像楼上给出来的这样杂乱无章.不妨通过一小段绘图的代码验证一下:
N = 25;x=linspace(0.52, 0.53, N);
y=linspace(-1e-9, 0, N);
z=linspace(1, 20000, N);
[x,y,z]=meshgrid(x,y,z);
F1=-2.*pi.*0.05415.*0.0000002.*sin(x).*sin((5.*pi./6)+x)-4./3.*pi.* ...
0.0000002.^3.*2000.*z.*9.78+2./3.*pi.*z.*9.78.*0.0000002.^3.* ...
(1820-1000).*(cos(x)).^3-2./3.*pi.*z.*9.78.*0.0000002.^3.* ...
(1820+1000)-pi.*z.*9.78.*0.0000002.^2.*(y+0.0000002.*cos(x)).* ...
(1820-1000).*(sin(x)).^2;
F2=-sin((5.*pi./6)+x).*besselk(0,(z.*9.78.*(1820-1000)./0.05415).^ ...
0.5.*0.0000002.*sin(x))+(z.*9.78.*(1820-1000)./0.05415).^0.5.* ...
y.*besselk(1,(z.*9.78.*(1820-1000)./0.05415).^0.5.*0.0000002.*sin(x));
clf;
p = patch(isosurface(x,y,z,F1,0,z));
set(p,'facecolor','w','EdgeColor','flat');
p = patch(isosurface(x,y,z,F2,0,z));
set(p,'facecolor','interp','EdgeColor','interp');
view(-10,35);
grid on;
xlabel('x')
ylabel('y')
zlabel('z')
绘制得到的曲面图形如下.
从图中我们很直观地初步得到几点结论:
(1)方程1的曲面对y几乎不相关.这一点其实比较容易理解,因为方程中y只出现在下面这一项中
(y+0.0000002.*cos(x))我们把这个式子的两项绝对值做一下对比:
max(abs(y(:)./(0.0000002.*cos(x(:)))))可以得到,最大值仅为0.0058,所以y在这个方程中完全可以忽略.
(2)方程1的x和z之间几乎是线性关系.
(3)两个曲面的交线近似是一条很平滑的空间曲线.
接下来,我们尝试画出这条曲线.代码如下:
syms x y zeq1=-2.*pi.*0.05415.*0.0000002.*sin(x).*sin((5.*pi./6)+x)-4./3.*pi.* ...
0.0000002.^3.*2000.*z.*9.78+2./3.*pi.*z.*9.78.*0.0000002.^3.* ...
(1820-1000).*(cos(x)).^3-2./3.*pi.*z.*9.78.*0.0000002.^3.* ...
(1820+1000)-pi.*z.*9.78.*0.0000002.^2.*(y+0.0000002.*cos(x)).* ...
(1820-1000).*(sin(x)).^2;
eq2=-sin((5.*pi./6)+x).*besselk(0,(z.*9.78.*(1820-1000)./0.05415).^ ...
0.5.*0.0000002.*sin(x))+(z.*9.78.*(1820-1000)./0.05415).^0.5.* ...
y.*besselk(1,(z.*9.78.*(1820-1000)./0.05415).^0.5.*0.0000002.*sin(x));
Z=solve(eq1,z);
eq = subs(eq2, z, Z);
hold on
h = ezplot(eq, [0.52 0.53], [-1e-9 0]);
X = get(h,'XData');
Y = get(h,'YData');
view(-10,35)
grid on
z = subs(Z, {x y}, {X Y});
x(z>20000)=NaN;
y(z>20000)=NaN;
z(z>20000)=NaN;
set(h, 'XData', X, 'YData', Y, 'ZData', z, 'linewidth', 2);
title('')
p=findall(gca,'type','patch');
a=0.2;
set(p,'Facealpha',a,'EdgeAlpha',a);
得到的绘图结果如下:
对上面的代码,简单说明几点:
(1)可以单独运行,也可以接在前面的代码后面运行(单独运行时不画曲面);
(2)
(3)为了突出显示交线,把两个曲面设为半透明,如果不想要这两个曲面,可以直接在把代码的最后两行换成:
delete(p)axis tight
即可删除曲面而只保留曲线.
以上程序在2007b上测试通过,其它环境暂未测试.
全部代码已作为附件上传,希望对楼主有帮助,如有问题可追问.
问题我已经贴到百度空间了,请大神帮我看看,来日一定报答
求大神帮我看看这到
请大神们帮我看看这个字有什么问题
求大神帮我看看文章有没有问题.
请帮我看看,
全部,已经做过的请帮我看看对不对.
请大神帮我看看这个鞋标正不正,
请大神帮我看看,这是什么动物.
请大神帮我看看,这道题这么证明可以么?
麻烦帮我看看,线性代数,基与基的转换问题题如图:请给我解析一下该如何转换,一二问就好,有具体过程更好!一二题已经解决了,麻烦帮看看第三题,
选错哦, 请帮我看看我的问题追加, 谢谢
江湖救急,运放问题LM385 懂电路的大神请帮我看看请帮我看一下,D4是一个放大电路,把0-3.3放大3倍,到0-9.9V D5是一个减法移位电路,放大平移到-8V到+8V, D6就是个跟随器. 先帮我看看这样给电
哪个数学大神帮我看看究竟哪里错了,
请帮我看看这几道题.
请帮我看看这是啥玉
麻烦帮我写一篇英文?题目是”my spare time”请各位可以帮我重新写一份,我也看了百度上的文章了,我同学已经用过了.谢谢,急用.周四之前要用,今天已经周一了.谢谢.不要太长,100到120
各路大神帮我看看这个英语
大神帮我看看什么意思?