设为首页收藏本站|繁體中文 快速切换版块

 找回密码
 立即加入
搜索
查看: 1817|回复: 0

是空间平面还是曲面?

[复制链接]

该用户从未签到

尚未签到

发表于 2007-12-6 11:20:04 | 显示全部楼层 |阅读模式

马上加入,结交更多好友,共享更多资料,让你轻松玩转电力研学社区!

您需要 登录 才可以下载或查看,没有账号?立即加入

×
最近研究一个保护问题需要绘制空间图形:; X  q: b9 u  s4 _
(2)x>=0 & y>4x     k=1.1
  \$ K6 h" ?& k% H1 i0 ~(3)y>=-4x &x<0
/ q  r1 Z: T) e. x$ s" I: @k=-1.6x/y-0.04x*24*9+1.1$ U4 G6 f- A, L5 |+ U# L, G3 o: W! h
(4)y<-4x0 }% i2 l9 w, i/ p
K=-0.1y/x+0.01y*24*9+1.14 E( F6 h; _" K$ F8 o7 k5 U" F
在x,y,k空间作出的图形可能是一个平面吗?
, y1 S) c9 B- w/ J5 S. d9 s还是一个有凸凹的曲面。我个人觉得因为这是一个二元一次方程,故应该是一个空间平面。
+ h( c' o( r: N+ B3 I! X7 k恳求大家发表下看法,谢谢 。% x$ Y* T  d/ S: v8 o
我用MATLAB编程如下:2 T$ V- R: Q7 z
x = -50 : 1 : 10;. R: ?. m# \; B3 `  w$ p9 T
y = 0 : 1 : 45;8 q+ v, U  W, r* |$ ?  C& S
nx = length(x);
6 E6 V3 e% S( r0 Gny = length(y);
0 z$ \7 b7 `' a, T* |+ M[X,Y] = meshgrid(x,y);9 H* ~& z0 ^$ L/ w+ {) r1 c
for n1 = 1 : ny0 w" b, L5 p7 b# G! o# W9 O
    for n2 = 1 : nx! B6 E6 m/ s+ n6 j' q  ^+ L
%        if Y(n1,n2) < 4.*X(n1,n2)/ g# d8 W- Q2 ^# H6 W. t+ ?
%             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;
) Q- u1 e2 Z+ |% Y7 N2 x3 |5 _%         end
. P3 Q2 \# j4 E+ s1 B        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )
! q/ p& }+ t4 _) t7 v            K(n1,n2) = 1.1;) g) W& \. J9 ]. C
        end
$ c  ~/ m9 {. X3 D        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )1 }, Z5 R8 r7 y
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;, q, v5 l, R5 N
        end' U* ~6 `* ?8 }( e+ o  \
        if Y(n1,n2) < ( -4.0*X(n1,n2) ), h& W8 o. `) ~6 l; {  h/ J! j
            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;* c3 y+ B6 L9 Q) F5 S/ y
        end0 U4 M) M. N! C1 p8 \0 ?
    end0 p; A# k* K" Q, _* }: E4 `' H
end0 a8 M: _" E- G& G& o& r
8 L# C5 H& \$ R; v/ _; l
subplot(2,1,1);( ]; ^5 H# V) s  V2 Y
mesh(X,Y,K);; u! U# C+ M% M3 b
% colormap([0 0 1])% 用蓝色绘图+ i% |8 K! E. P% a1 A
xlabel('--x--');ylabel('--y--');zlabel('--k--');8 A: E  Y; D8 p8 l% e3 l
title('全部画出的图形');  H) [& ^& p8 [
hold on % p. C1 k3 ^* C0 i. t4 [1 ?% F2 n

8 K1 a8 ~  D* W  [1 B%将不画图部分处理为零,与第一张全图比较看看
3 R9 {' |) z1 j  a4 y% Zfor n1 = 1 : ny
( s' Z* e4 {( m- n4 Y    for n2 = 1 : nx8 k/ I) E# o; g9 E% S1 a) e
        if Y(n1,n2) < ( 4*X(n1,n2) )
9 v  X; n/ M* t- G8 W; A5 B             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;7 O6 K, i9 G4 E& q% W) D) M; n1 @/ d" ]
        end8 n8 `& x8 r5 G$ w
        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )/ w* f7 c6 Y  V0 |  @* Y9 X: M3 p
            K(n1,n2) = 1.1;
1 Z+ v. X- P; |1 K        end
: D8 y3 p& J2 E        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )/ D+ V) ]6 o$ x1 u8 n) v
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;
8 L( Z# k+ i1 i2 w) U        end  ^' ~  |, q( l
        if Y(n1,n2) < ( -4.0*X(n1,n2) )0 }5 w1 [. `' a- E
            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;
* z# W" C7 q* {6 B3 I        end
. T, R6 m, p) B    end
6 Y) g/ A2 X2 r+ Kend3 L9 \: O6 M. o. W
/ F% B( |3 [% z* T0 |
subplot(2,1,2);
; ?! A) Y3 t* I& T2 x. Vmesh(X,Y,K);" W- `1 ]' `! P  B7 J! I
colormap([1 0 0]) 4 ^5 b5 o9 i$ g+ d
xlabel('--x--');ylabel('--y--');zlabel('--k--');4 J5 d  L* s0 r6 ~
title('Y<4x 时不画,处理为零的图形');
# a: P( l' @9 P" b% e# q9 P' a2 N" @' _6 D
出来图形为:“我不知道怎么贴上来”,就是一个空间平面。7 \) P' l3 c) Y( m  G
请大家发表下看法和意见,多谢了啊
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

小黑屋|手机版|APP下载(beta)|Archiver|电力研学网 ( 赣ICP备12000811号-1|赣公网安备36040302000210号 )|网站地图

GMT+8, 2025-5-23 18:30

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回复 返回顶部 返回列表