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

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

是空间平面还是曲面?

[复制链接]

该用户从未签到

尚未签到

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

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

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

×
最近研究一个保护问题需要绘制空间图形:3 S! ?$ C1 x; c/ G, b
(2)x>=0 & y>4x     k=1.1# p& Q7 p; l! l
(3)y>=-4x &x<0: S9 _& l& i# ~
k=-1.6x/y-0.04x*24*9+1.1
" \  f, A$ i& B& t(4)y<-4x7 u+ U' G/ u# j
K=-0.1y/x+0.01y*24*9+1.1. j2 A$ G( ?/ R& a9 E
在x,y,k空间作出的图形可能是一个平面吗?- G1 M: F! w, `4 Z2 p7 f3 f; T4 B
还是一个有凸凹的曲面。我个人觉得因为这是一个二元一次方程,故应该是一个空间平面。
! @6 ]- D; f5 ?* S恳求大家发表下看法,谢谢 。4 I& h8 ?: @" O. l$ b5 y  {# @
我用MATLAB编程如下:
. X' j5 z) h; Y" B+ [x = -50 : 1 : 10;/ p7 g: v) X* E# }+ }. [
y = 0 : 1 : 45;% y5 S8 X% v. D4 W  H
nx = length(x);5 M: z6 D. k* B7 u
ny = length(y);7 D: Z& Z: E  H6 \$ i
[X,Y] = meshgrid(x,y);7 K6 Y8 ]2 [- r# k/ {8 j
for n1 = 1 : ny
" z, ^" P' C  L% H0 a5 m    for n2 = 1 : nx9 v% p9 @) m7 `( R, j
%        if Y(n1,n2) < 4.*X(n1,n2)! n- f0 y; X: l5 V& |; t
%             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;# J/ t/ X  m' F
%         end
7 B" |$ E. B8 U6 E0 J$ R6 G        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )
, f# c2 x& y4 k6 v( `! D) Q+ ^            K(n1,n2) = 1.1;
/ f; P1 }( O- r$ T9 C$ I$ _) @        end# r" g7 b, S' N9 ]' `
        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )# @  W/ c! a& q
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;
% \! F0 n; \* A        end
1 U. ?8 m* M. B. ]# x: A        if Y(n1,n2) < ( -4.0*X(n1,n2) )
  {$ q9 y8 r+ ]) a" G# p% r9 p            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;
. H2 `  j/ j0 n        end; A  ~9 P/ B9 e( a" O
    end5 i0 n- V+ z" q
end4 K! A- J- l& U# I  {

0 a+ E2 {3 G5 ]# @5 zsubplot(2,1,1);) z2 i' m1 u% i. N3 m
mesh(X,Y,K);
" z; ^* {8 h2 _& J% colormap([0 0 1])% 用蓝色绘图0 X- F8 ^( c3 x. y6 [7 C$ R
xlabel('--x--');ylabel('--y--');zlabel('--k--');
5 r9 q! ]. N6 O! p6 o) O3 }! ctitle('全部画出的图形');
1 M0 D$ U! g4 Mhold on
. H9 O7 m0 x4 `7 H. i. `
0 e  D/ r2 o. X2 Q%将不画图部分处理为零,与第一张全图比较看看
: F2 r8 v7 W$ x8 \! F% Lfor n1 = 1 : ny
) x; q, C2 V, e5 G' {: [0 C2 q    for n2 = 1 : nx* N& V* ?* C9 t& U8 W; h- W
        if Y(n1,n2) < ( 4*X(n1,n2) )' m( l0 @% ~/ `) `% W9 ?, D; x1 E
             Y(n1,n2)=0;X(n1,n2)=0;K(n1,n2) = 0;7 [" Y+ ?* n; U1 ?! B
        end5 {9 x$ E0 A5 F2 |
        if ( X(n1,n2)>=0 ) & ( Y(n1,n2)>4*X(n1,n2) )' ~. a. r: J/ U" j9 y7 @. N
            K(n1,n2) = 1.1;; M# Y+ v* n- q
        end
3 i8 L& j  h. ]) Y/ r        if ( Y(n1,n2) >= -4*X(n1,n2) ) & ( X(n1,n2)<0 )* G+ o6 A5 E( p. l: f" O
            K(n1,n2) = (-1.6*X(n1,n2))/Y(n1,n2) - (24*9*0.04)*X(n1,n2) + 1.1;
8 r) O% s: e- V- d5 t        end+ ^0 e/ ^! [  F1 ^$ ?
        if Y(n1,n2) < ( -4.0*X(n1,n2) ), T- p* F$ [" o% k# f  u
            K(n1,n2) = (-0.1*Y(n1,n2))/X(n1,n2) + (24*9*0.01)*Y(n1,n2) + 1.1;
4 C( W/ r, ?. _% \) _        end
5 J0 _' Y" U  L7 B9 C0 ~' [    end
" p* z+ l" D, @! q0 W; C! xend1 O. _) N) V8 N" p
! e' R) e: V5 }9 V
subplot(2,1,2);
! ]* f9 \' }% @' J$ O8 u2 qmesh(X,Y,K);" B) ^% H7 c: m/ g1 F% v
colormap([1 0 0])
! Z" g/ Z. x4 p: {  Ixlabel('--x--');ylabel('--y--');zlabel('--k--');
2 [; u0 Z9 ?& N! G. P' G  b( gtitle('Y<4x 时不画,处理为零的图形');. {& I8 a- i1 c: t/ `6 T5 U5 Z% A' N

2 T7 z; N8 i! Z0 p  M' l+ }) ~0 i" [0 N出来图形为:“我不知道怎么贴上来”,就是一个空间平面。, s+ @3 e$ z, m5 m! A
请大家发表下看法和意见,多谢了啊
"真诚赞赏,手留余香"
还没有人打赏,支持一下
楼主热帖
帖文化:【文明发帖 和谐互动】 社区精神:【创新、交流、互助、共享】
您需要登录后才可以回帖 登录 | 立即加入

本版积分规则

招聘斑竹

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

GMT+8, 2025-2-23 05:48

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

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