数学建模社区-数学中国

标题: mathematica如何把不同取值范围的曲线画到同一个坐标系 [打印本页]

作者: zhb    时间: 2012-7-5 20:01
标题: mathematica如何把不同取值范围的曲线画到同一个坐标系
求解求解,,,,本人刚开始学,,这个暑假要做课题,,,希望大家帮帮忙,
作者: 风筝11    时间: 2012-7-5 21:18
PlotRange->All
作者: zhb    时间: 2012-7-5 22:19
风筝11 发表于 2012-7-5 21:18
. ~; F3 X6 c! ?5 i" o' jPlotRange->All

. F' a0 u1 q* e- |4 @8 g2 l- R谢谢哈哈!!!!
作者: zhb    时间: 2012-7-5 22:22
风筝11 发表于 2012-7-5 21:18
, z# g: d6 r( M, l* ZPlotRange->All
: D+ G, f% f0 i5 r3 b) D# l3 [7 A
那个。。我试了下,感觉不对的么。。。
作者: 风筝11    时间: 2012-7-5 22:25
每个都加 给我 看看你的代码吧
作者: 风筝11    时间: 2012-7-5 22:26
把几张图显示在一张图是用Show
作者: zhb    时间: 2012-7-5 22:39
风筝11 发表于 2012-7-5 22:26 1 \' l8 ?& d3 c
把几张图显示在一张图是用Show
4 ~+ j) w+ c( T7 O* X  Y" S
这个我也试过了不行哎。。因为画的图有根号,所以两个函数取值范围不一样,超过了这个范围就提示错误。真心不知道该怎么弄。
作者: 风筝11    时间: 2012-7-5 23:38
把代码给我 看看啊
作者: zhb    时间: 2012-7-5 23:53
风筝11 发表于 2012-7-5 23:38
* c) h; r- E* H+ L2 |  D# r# X* x把代码给我 看看啊

; @- v2 S/ Y' {( ]  h7 |a11 = 3.609*10^7; a12 =
' B. p+ M: E2 Z, _# e% P4 g    2.284*10^8; a111 = 1.437*10^8; a112 = 6.630*10^8; a123 = -3.151*10^9; Q11 \( e# ]* J: X+ N8 Z
= 9.151*10^-2; Q12 = -4.013*10^-2; Q44 = 7.674*10^-2; S11 = 7.963*10^-12; S12 \
" N/ Z$ ]% V. ~* x7 K= -3.351*10^-12; S44 = 25.716*10^-12;
8 ]4 z" a7 A' |& y3 o1 ma11' = a11 + 0.5*((Q11^2 + Q12^2)S11 - 2Q11*Q12*S12)/(S11^2 - S12^2);
' X9 t0 E' e3 l9 c+ E' l! R+ ca33' = a11 + Q12^2/(S11 + S12);/ [& `7 f' C4 O0 ]% w# h& _
a12' = a12 - 1/(S11^2 - S12^2)*((Q11^2 - Q12^2)S12 - 2*Q11*Q12*S12) + \/ Y$ G" _, R, g6 ^# a+ E  z
Q44^2/(2S44);2 @! }/ j  H0 ?/ K" \
a13' = a12 + Q12(Q11 + Q12)/(S11 + S12);
2 ]4 \$ T- t9 }' c" R) Qa1' = a1 - S(Q11 + Q12)/(S11 + S12);
2 O. e) Y3 a% u% O( E9 |* Ha3' = a1 - 2*S*Q12/(S11 + S12);
; I6 X4 ?1 z. U: z5 d/ oa1 = (25 - 360)/(2*8.85*10^-12*3.701*10^5);& g! Q3 P+ R4 p
P3' = 1/(3a111)*(-a3' + Sqrt[a33'^2 - 3a111*a3']);9 w/ {6 l" U2 M
P1' = 1/(6a111 + 6a112)*((-2a11' - 7 j7 S$ q3 K; C* g9 ^2 I* A
    a12') - Sqrt[(2a11' + a12')^2 - 12a1'*(a111 + a112)]);5 b. D* V$ d' Q3 g$ b
Plot[Sqrt[P3'], {S, -0.016, 0.005}]
" i2 J+ N" F: {2 P* hPlot[Sqrt[-P1'], {S, -0.005, 0.016}]! o' E( n; ?, W5 [# t3 P
  {8 w9 r9 c! h# `  [0 b
2 \3 u, {( I: S8 H3 Q6 r9 L

5 [6 g  ~" R* Q3 f, l
3 P0 P& u. z. d7 x9 B) c8 C- k% o4 Y5 {- O  ?9 U2 t# n
把两条合并。
作者: zhb    时间: 2012-7-6 00:23
zhb 发表于 2012-7-5 23:53
8 x9 f) l$ `1 N. g8 V# ~3 J. S6 @a11 = 3.609*10^7; a12 =
$ i9 i8 ]2 p4 s! @    2.284*10^8; a111 = 1.437*10^8; a112 = 6.630*10^8; a123 = -3.151*10^9 ...

, P" B4 s! ~+ V4 l还有啊,
3 R5 u# ~- U3 u4 aContourPlot3D[(x^2 + (9 y^2)/4 + z^2 - 1)^3 - x^2 z^3 - (9 y^2 z^3)/2 q2 x# S& o0 P% _9 U
   80 == 0, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {z, -1.5, 1.5}, 6 k8 d/ b0 v6 @9 b6 e) L9 X# f% Z" U. P- B: ~* @. H
PlotPoints -> 30, Axes -> False, Lighting -> False, ( }$ M1 W  T5 @# p. J$ z
) I2 |7 R; p8 Q; JContourStyle -> {RGBColor[1, 0.0, 0.0]}]这个是画心形图的。。。我的mathematica5.0好像画不出。。。
作者: 风筝11    时间: 2012-7-6 16:13
zhb 发表于 2012-7-5 23:53 8 s9 x- C7 {6 e, @
a11 = 3.609*10^7; a12 = ' A3 D! N2 T; I* [( X% w9 Q) a
    2.284*10^8; a111 = 1.437*10^8; a112 = 6.630*10^8; a123 = -3.151*10^9 ...
5 \* M( A& t4 T% P' Q, b
C:\Users\Administrator\Desktop
作者: 风筝11    时间: 2012-7-6 16:15
03 02 01 6 z* T. ^" K8 J( O$ j* s5 n
您的代码有点乱 我不知道是不是这样。。。
作者: 风筝11    时间: 2012-7-6 16:18
我的能画 8.0的  可能是版本问题
2 w9 }/ D6 u+ g* E好像这些语句不太好 运行时会出现橙色的信息# G, x% K8 L- e% u( b/ b/ f! a
1 M  f9 i6 v$ Z
04.png
作者: 风筝11    时间: 2012-7-6 16:25
前面那3个图顺序颠倒了
作者: zhb    时间: 2012-7-6 21:44
风筝11 发表于 2012-7-6 16:15 , R* p- h. U4 s: L* x5 N
您的代码有点乱 我不知道是不是这样。。。
' t. U: {- P8 E( `3 {# B
高手的说,哈哈哈,不错不错,你是学数学的么?
作者: zhb    时间: 2012-7-6 21:46
风筝11 发表于 2012-7-6 16:25 8 J& n# O- e* [; t# ^$ Y' m. |4 ]3 R
前面那3个图顺序颠倒了

& F5 o1 Q! Y- @7 ~/ w. a- u。。。话说你还是位女生!给力的!!
作者: 风筝11    时间: 2012-7-6 23:05
过奖了~~因为我也百度过那个心: [/ S6 H% u5 _4 T( p* |
我学通信的。。因为想参加建模竞赛,所以学了点mathematica 那些老师上课都讲过。。
作者: zhb    时间: 2012-7-6 23:06
风筝11 发表于 2012-7-6 23:05 ) [& |9 a! m/ o1 E: u
过奖了~~因为我也百度过那个心
, q1 V8 D1 A8 v/ O$ V& u我学通信的。。因为想参加建模竞赛,所以学了点mathematica 那些老师上课都 ...
  _/ T( q0 z6 d4 h1 Y) w
干嘛不学matlab呢?感觉那个建模更有效,你也是大二的么?
作者: 风筝11    时间: 2012-7-6 23:07
zhb 发表于 2012-7-6 21:44
; y0 y! N: h2 F) g9 m高手的说,哈哈哈,不错不错,你是学数学的么?

# A, x/ ?9 t% ^' G/ R  g) K您的代码应该要 整齐一点这样的话有错比较好发现,看上去也比较舒服
作者: zhb    时间: 2012-7-6 23:08
风筝11 发表于 2012-7-6 23:07 0 m2 f8 O2 N) \
您的代码应该要 整齐一点这样的话有错比较好发现,看上去也比较舒服
& u" i' s5 v9 K4 c5 A2 {( j: x
恩恩,是的,今天把代码整理了。。不要用“您”这个称呼,感觉怪怪的。
作者: 风筝11    时间: 2012-7-6 23:09
zhb 发表于 2012-7-6 23:08 * h# ~1 K8 s5 F  U: v% Q
恩恩,是的,今天把代码整理了。。不要用“您”这个称呼,感觉怪怪的。
+ R& U  p9 a3 N6 W
有时候会用这称呼 嘿嘿
作者: 风筝11    时间: 2012-7-6 23:10
zhb 发表于 2012-7-6 23:06 0 [: J: R+ A: c
干嘛不学matlab呢?感觉那个建模更有效,你也是大二的么?
. j' ]1 I- Q# B- C4 P' J5 p
我们老师讲的是mathematica      matlab也懂一些 大一大二都学过基础,可能暑假看吧 是大二 不过现在看来应该是大三了
作者: zhb    时间: 2012-7-6 23:11
风筝11 发表于 2012-7-6 23:09 ( [' A4 @9 s  g! ~4 ?  r
有时候会用这称呼 嘿嘿
0 Q) i, @* B: B9 I
哈哈,感觉这论坛不错的,能认识彼此都在学习软件的人,
作者: zhb    时间: 2012-7-6 23:12
风筝11 发表于 2012-7-6 23:10
! U/ l" K  p: H( D# M5 C/ J我们老师讲的是mathematica      matlab也懂一些 大一大二都学过基础,可能暑假看吧 是大二 不过现在看来 ...
% i5 P3 \# z) y! Q
真好,我也是准大三了。
作者: 风筝11    时间: 2012-7-6 23:12
zhb 发表于 2012-7-6 23:11 3 J& L, n: j9 C
哈哈,感觉这论坛不错的,能认识彼此都在学习软件的人,

' t2 \0 o' i1 i& p9 \2 P) L* O* L还行吧 这个论坛也有很多不太好的地方
作者: zhb    时间: 2012-7-6 23:14
风筝11 发表于 2012-7-6 23:12 + m2 ~: C# M# ^
还行吧 这个论坛也有很多不太好的地方
3 i# ]6 ?& G* T. }9 B
怎么说?是帖子太少了么?感觉人员不太够。
作者: 风筝11    时间: 2012-7-6 23:15
zhb 发表于 2012-7-6 23:14
+ L0 \0 a* s: \# s怎么说?是帖子太少了么?感觉人员不太够。

0 o) h- W$ {+ b就是搜索不太方便 还有就是隔60秒才能再发帖子 觉得不是很好
作者: zhb    时间: 2012-7-6 23:16
风筝11 发表于 2012-7-6 23:15
" I# C7 O7 V; A( u% ]8 [/ t就是搜索不太方便 还有就是隔60秒才能再发帖子 觉得不是很好
/ ~! b% a+ ]8 f! u$ Q
是的是的,我学校内外的论坛感觉也是很多的权限啊什么的。
作者: 风筝11    时间: 2012-7-6 23:19
zhb 发表于 2012-7-6 23:16 0 r7 h+ O; x: U. ^
是的是的,我学校内外的论坛感觉也是很多的权限啊什么的。
! p8 m2 ?6 n; N0 |5 e: v& r3 V% V7 `
可能不想有太多杂帖 而是注重交流学习吧
作者: zhb    时间: 2012-7-6 23:21
风筝11 发表于 2012-7-6 23:19
; i% T0 i) R& E可能不想有太多杂帖 而是注重交流学习吧
% C9 H% V8 l6 C* j* ~1 f" p
这么说有理的,唉,暑假在学校做些课题。
作者: zhb    时间: 2012-7-6 23:22
风筝11 发表于 2012-7-6 23:19
& `2 _& q5 k0 Z3 [/ A, W! z可能不想有太多杂帖 而是注重交流学习吧

, R) B& U$ `# _, Q7 x数学建模都感觉不是很有兴趣所有没有参加。
作者: sycxyc    时间: 2012-7-7 15:44
zhb 发表于 2012-7-6 00:23
# V! b7 M- q; `! e: S还有啊,
! q0 B1 D1 ]' k, kContourPlot3D[(x^2 + (9 y^2)/4 + z^2 - 1)^3 - x^2 z^3 - (9 y^2 z^3)/  @  ?9 w3 ~/ }3 V
   80 == 0, {x, -1 ...

6 N3 H+ ]) e# ^9 ~; C5.0代码要有些变化
: K! w2 H  p$ |! v4 K- [# X
7 t7 _. X  J/ i7 [( @- _<< Graphics`ContourPlot3D`+ {: T2 l3 J2 q( F# _6 N/ y: P1 g
然后再
5 F! |' k- p$ u; z# }ContourPlot3D[(x^2 + (9 y^2)/4 + z^2 - 1)^3 - x^2 z^3 - (9 y^2 z^3)/80, {x, -1.5, 1.5}, {y, -1.5, 1.5}, {z, -1.5, 1.5}, PlotPoints -> {9, 9},ContourStyle -> {RGBColor[1, 0.0, 0.0]}]
作者: zhb    时间: 2012-7-7 17:50
sycxyc 发表于 2012-7-7 15:44
1 w9 L  ?0 Q$ u( N- a3 k- h( I5.0代码要有些变化5 r( g1 T& S6 O: W" d
2 L9 g- P$ ?6 S  j* y  g2 K1 M
{9, 9},ContourStyle -> {RGBColor[1, 0.0, 0.0]}]

! B# z6 h+ U" k+ b1 m( r  ^哇,高手啊,谢啦!!!
作者: qingxiong    时间: 2012-7-9 15:59
跟着学习了一下……
作者: zhb    时间: 2012-7-9 20:34
qingxiong 发表于 2012-7-9 15:59
9 q, K% v! c; q: E6 a+ m" D跟着学习了一下……

; e) G& ^/ K4 Q$ D, `同学习。。。。。。。。。。。。。。。。




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5