为什么出不来图像，貌似只有坐标系？
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('被积函数图像') 原帖由 trl-008 于 2007-5-22 10:24 发表 http://www.chinavib.com/forum/images/common/back.gif
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了
回复 #2 eight 的帖子
这样啊，谢谢啦回复 #2 eight 的帖子
0.0622 0.0677 0.0699 0.0689 0.0645 0.0568 0.0461 0.0329 0.01770.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那个坐标轴取得刻度有问题。
怎么设置？ 原帖由 trl-008 于 2007-5-22 13:08 发表 http://www.chinavib.com/forum/images/common/back.gif
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 ...
你这串数的前面还有一个乘数的
回复 #5 eight 的帖子
不好意思，太粗心了！！再问一个问题，想这种原函数可以积分吗？
谢谢！！:@P
回复 #5 eight 的帖子
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)
那这种情况该怎么处理？谢谢了 原帖由 trl-008 于 2007-5-23 19:00 发表 http://www.chinavib.com/forum/images/common/back.gif
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
页:
[1]