为什么出不来图像，貌似只有坐标系？

trl-008

clear;
w=0:0.01:1.27;
h=0:0.01:1.27;
for i=1:128
    for j=1:128
Y(i,j)=2484016669471539./4096.*exp((-2760./353.*w(i).^4-500./353).*(81-h(j))).*sin(510322903951331./309485009821345068724781056.*(-6708298677659438208253952.*w(i).^4+176380388135181248731545600-22483186824942380862406656.*w(i).^8).^(1./2).*(81-h(j))).*cos(w(i).*(2602697954573681./2199023255552-1667./100.*h(j)))./(-6708298677659438208253952.*w(i).^4+176380388135181248731545600-22483186824942380862406656.*w(i).^8).^(1./2);
end
end
surf(w,h,Y)
hold on
grid on
xlabel('w')
ylabel('h')
zlabel('y')
title('被积函数图像')

eight

原帖由 trl-008 于 2007-5-22 10:24 发表

                               
登录/注册后可看大图

clear;
w=0:0.01:1.27;
h=0:0.01:1.27;
for i=1:128
    for j=1:128
Y(i,j)=2484016669471539./4096.*exp((-2760./353.*w(i).^4-500./353).*(81-h(j))).*sin(510322903951331./3094850098213450687247810 ...

估计你的Y的值太小了，10^{-51}数量级，近似为0了

trl-008

这样啊，谢谢啦

trl-008

0.0622    0.0677    0.0699    0.0689    0.0645    0.0568    0.0461    0.0329    0.0177
    0.0463    0.0503    0.0519    0.0511    0.0477    0.0420    0.0340    0.0242    0.0130
    0.0068    0.0072    0.0072    0.0068    0.0062    0.0053    0.0041    0.0028    0.0015
   -0.0362   -0.0396   -0.0413   -0.0409   -0.0386   -0.0342   -0.0279   -0.0201   -0.0109
   -0.0606   -0.0660   -0.0684   -0.0674   -0.0632   -0.0557   -0.0453   -0.0324   -0.0174
   -0.0540   -0.0585   -0.0602   -0.0590   -0.0549   -0.0481   -0.0389   -0.0276   -0.0148
   -0.0199   -0.0210   -0.0211   -0.0201   -0.0182   -0.0155   -0.0121   -0.0083   -0.0043
    0.0240    0.0268    0.0284    0.0287    0.0276    0.0249    0.0206    0.0151    0.0083
    0.0549    0.0601    0.0625    0.0619    0.0582    0.0515    0.0420    0.0301    0.0162
    0.0572    0.0619    0.0636    0.0623    0.0580    0.0508    0.0409    0.0290    0.0154
    0.0304    0.0322    0.0323    0.0308    0.0279    0.0238    0.0186    0.0128    0.0066
   -0.0104   -0.0124   -0.0139   -0.0148   -0.0149   -0.0140   -0.0121   -0.0091   -0.0051
   -0.0433   -0.0478   -0.0500   -0.0499   -0.0472   -0.0421   -0.0345   -0.0247   -0.0133
   -0.0520   -0.0563   -0.0580   -0.0568   -0.0529   -0.0463   -0.0372   -0.0262   -0.0137
   -0.0340   -0.0360   -0.0361   -0.0345   -0.0313   -0.0267   -0.0209   -0.0142   -0.0072
   -0.0016   -0.0005    0.0007    0.0020    0.0030    0.0036    0.0037    0.0032    0.0020
    0.0266    0.0296    0.0314    0.0317    0.0303    0.0271    0.0223    0.0160    0.0084
    0.0367    0.0398    0.0410    0.0402    0.0374    0.0327    0.0262    0.0182    0.0092
    0.0271    0.0286    0.0287    0.0274    0.0248    0.0210    0.0163    0.0110    0.0053
    0.0070    0.0066    0.0058    0.0047    0.0035    0.0023    0.0012    0.0004    0.0000
   -0.0106   -0.0121   -0.0131   -0.0134   -0.0130   -0.0118   -0.0098   -0.0069   -0.0034
   -0.0174   -0.0189   -0.0195   -0.0191   -0.0177   -0.0154   -0.0122   -0.0082   -0.0037
   -0.0136   -0.0143   -0.0143   -0.0136   -0.0122   ......................
算了一下，数量级没那么小。
估计是Y那个坐标轴取得刻度有问题。
怎么设置？

eight

原帖由 trl-008 于 2007-5-22 13:08 发表

                               
登录/注册后可看大图

0.0622    0.0677    0.0699    0.0689    0.0645    0.0568    0.0461    0.0329    0.0177
    0.0463    0.0503    0.0519    0.0511    0.0477    0.0420    0.0340    0.0242    0.0130
    0.0068    0 ...

你这串数的前面还有一个乘数的

trl-008

不好意思，太粗心了！！
再问一个问题，想这种原函数可以积分吗？
谢谢！！:@P

trl-008

clear;
w=1.27:0.01:3;
h=0:0.01:1.73
for i=1:174
    for j=1:174
Y(i,j)=606449382195.20000000000000000000.*exp((-7.8186968838526912181303116147309.*w(i).^4-1.4164305949008498583569405099150).*(1.-1.*h(j))).*sinh(.16489422355089917061349150141643e-11.*(6708298677659438545710450.*w(i).^4-176380388135181260069835000.+22483186824942380862406656.*w(i).^8).^(1./2).*(1.-1.*h(j))).*cos(w(i).*(16.670000000000000000000000000000-16.670000000000000000000000000000.*h(j)))./(6708298677659438545710450.*w(i).^4-176380388135181260069835000.+22483186824942380862406656.*w(i).^8).^(1./2);
end
end
surf(w,h,Y)
hold on
grid on
xlabel('w')
ylabel('h')
zlabel('y')
title('被积函数图像')  
%%%%%%%%%%%%%%%%%%%
是不同的呀！算出Y(i,j)是复数，是不是这个原因？导致程序运行后出现错误：
??? Error using ==> surf
X, Y, Z, and C cannot be complex.

Error in ==> myint2 at 9
surf(w,h,Y)
那这种情况该怎么处理？谢谢了

eight

原帖由 trl-008 于 2007-5-23 19:00 发表

                               
登录/注册后可看大图

clear;
w=1.27:0.01:3;
h=0:0.01:1.73
for i=1:174
    for j=1:174
Y(i,j)=606449382195.20000000000000000000.*exp((-7.8186968838526912181303116147309.*w(i).^4-1.4164305949008498583569405099150). ...

help real