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

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

是空间平面还是曲面?

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
最近研究一个保护问题需要绘制空间图形:7 ]) J3 P) |) b8 r& |
(2)x>=0 & y>4x     k=1.1* F6 s/ B) Z8 b6 v: i) \
(3)y>=-4x &x<0: }# d1 Y) U' X& G$ ]% W6 L; f
k=-1.6x/y-0.04x*24*9+1.1
6 ^; c6 |3 T+ j6 n(4)y<-4x
1 Q. h! [8 P% \K=-0.1y/x+0.01y*24*9+1.1
6 V) m& T- v+ B& c在x,y,k空间作出的图形可能是一个平面吗?- i- h3 b5 D8 G+ l5 @% O2 q
还是一个有凸凹的曲面。我个人觉得因为这是一个二元一次方程,故应该是一个空间平面。% C4 J0 e: r8 T! h; V+ K
恳求大家发表下看法,谢谢 。
5 q" P% B. J- I0 p* ~+ l2 v2 p我用MATLAB编程如下:
4 s/ c  p7 b4 N% l4 ~x = -50 : 1 : 10;% _& Q1 Z1 t. L9 U
y = 0 : 1 : 45;
  b5 L8 s9 R7 W, H7 W8 B7 b  mnx = length(x);* W' R) W7 Z; _9 i
ny = length(y);9 s" e% W! z& F, t5 x- p* p
[X,Y] = meshgrid(x,y);
8 x2 r/ O5 K0 X% z+ gfor n1 = 1 : ny6 s: U3 O& v$ i
    for n2 = 1 : nx
0 }: }) C+ A) _, k! w' {%        if Y(n1,n2) < 4.*X(n1,n2)
. A) T& b- `- w( s6 }%             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;6 ~# {6 R# L; g) s7 ], [4 Z
%         end) S5 k3 Q: s  K! q
        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )
( L) B- z# v$ l9 z. X7 {  n5 X            K(n1,n2) = 1.1;; U; Q6 i# ^" k( |' N/ Q
        end  K! Y" j5 }, `$ u1 Z" X
        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )9 e$ B$ m8 H/ v6 c0 m
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;. i* S7 O; v8 o0 C; G) A  K
        end/ \8 M  g7 U# G& [) a* v" Z9 B/ O/ F
        if Y(n1,n2) < ( -4.0*X(n1,n2) )
# V3 v- {) {+ |            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;
$ o* v0 i) H2 B        end0 V( y9 n9 X9 k& W1 ^. P
    end
+ d; ^3 k5 D. v+ t( v$ {end! y0 [) g& z% a/ L  k& `
  Z7 Q( F0 b2 H
subplot(2,1,1);7 H. w0 y( \4 p) {- B( s! V6 L
mesh(X,Y,K);7 P$ C$ q- U0 j2 V) n
% colormap([0 0 1])% 用蓝色绘图
7 B  \0 o$ s3 K2 pxlabel('--x--');ylabel('--y--');zlabel('--k--');* x' X8 d6 o6 ?: D3 A  |: N
title('全部画出的图形');
6 M0 F+ l% h$ E1 R  g: Qhold on
3 ?& p8 o1 |3 |5 x0 B
5 ~1 l8 @! ^; P3 [$ s! G%将不画图部分处理为零,与第一张全图比较看看+ I5 M# ]$ O) w& Q6 ^0 G+ ^
for n1 = 1 : ny
; A4 r$ O0 O+ Z' Y, x; n    for n2 = 1 : nx
0 I" x' @4 g, @, J$ b$ f        if Y(n1,n2) < ( 4*X(n1,n2) )
* W$ R' f3 W! [/ F' p             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;2 L1 X* ]) k: x1 u! z
        end$ Z/ X: g; G! q$ _# n
        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )
% r" ?; s9 |: q% G: }# j            K(n1,n2) = 1.1;4 J1 k/ K0 k8 _/ e& a6 F
        end
1 ~* j1 q7 w  ?$ _        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )2 _% z8 e9 b4 w% e; s
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;9 P& j6 \) T1 i! g
        end
  q# i* |) Z5 f* M" R        if Y(n1,n2) < ( -4.0*X(n1,n2) )( v+ e8 J6 ~. y. T& h' V
            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;
8 }$ {/ i: y6 C+ `        end
' r  A* v5 g- W$ ^- u4 k    end
! t. S6 X- H) S6 g/ ?$ K7 r3 Gend3 _; _7 Q) o1 D" m+ a

/ f5 ]( d5 l& ?/ A6 c' ^0 msubplot(2,1,2);
" d  S/ C4 b0 a! V0 Y0 ?mesh(X,Y,K);0 \0 d2 k; F) h2 ~6 j
colormap([1 0 0]) % k' u7 i% R6 {( I3 r+ |# b% O
xlabel('--x--');ylabel('--y--');zlabel('--k--');2 {  a) H7 S/ L' f0 r
title('Y<4x 时不画,处理为零的图形');$ E* ^& g$ F* i- [

( r' `; l, h3 l! P" {" B出来图形为:“我不知道怎么贴上来”,就是一个空间平面。
: P7 X6 u( q5 M% j+ ]! H请大家发表下看法和意见,多谢了啊
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

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

GMT+8, 2026-4-30 18:51

Powered by Discuz! X3.5 Licensed

© 2001-2026 Discuz! Team.

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