机械社区

 找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 2939|回复: 8

螺杆转子刀具设计(2) matlab 解法

[复制链接]
发表于 2015-6-4 11:00:10 | 显示全部楼层 |阅读模式
螺杆转子刀具设计 数学模型7 e. C8 ]" J  c1 c3 F+ E

" [/ Q9 Z% x! ~3 j9 E. L解如下方程:1679999.7*sin(x)-1679997.5*cos(y)*sin(x)-1499997.5*sin(y)*cos(x)-75000*sin(x)^3-179999.4000005*y*cos(y)*cos(x)+179999.4000005*y*sin(y)*sin(x)+75000*cos(y)*cos(x)*sin(x)=0. D6 B! k0 f: a3 l( W
求当x=-1.361时  y=?  ( }9 _: z) a% P
>> syms x y
+ ?2 h/ G7 B2 D. x6 Q>> f=1679999.7*sin(x)-1679997.5*cos(y)*sin(x)-1499997.5*sin(y)*cos(x)-75000*sin(x)^3-179999.4000005*y*cos(y)*cos. |/ M, C: p4 c
1 ?: u2 W3 c- S- ~: [
(x)+179999.4000005*y*sin(y)*sin(x)+75000*cos(y)*cos(x)*sin(x)& _" T$ S$ Z2 j' r% V9 P

1 G. y/ x+ X" k5 D1 [' o5 qf =  S4 Z5 }" H/ r8 @8 b: ~6 ]

3 H$ ^% W( Q& z3 I(7215543768789811*sin(x))/4294967296 - (2999995*cos(x)*sin(y))/2 - (3359995*cos(y)*sin(x))/2 - 75000*sin(x)^3 + 75000*cos(x)*cos(y)
1 [- i2 p# v1 N4 r( c/ Z% Z# k$ {' u, H* [2 l. d
*sin(x) - (6184732290414159*y*cos(x)*cos(y))/34359738368 + (6184732290414159*y*sin(x)*sin(y))/343597383683 p5 @  |9 {( n: ~( q
, ~' T( V1 m, @( J5 R
>>  subs(f, x, -1.361)# E! m$ J3 F+ h! ]+ L
, k4 P0 U6 U4 @0 K( M/ Z) N, q5 y
ans =
' ]* w  N8 W/ ^8 h5 N5 Q
9 H3 x0 n% ?3 t& a; J8 }(3359995*sin(1361/1000)*cos(y))/2 - (2999995*cos(1361/1000)*sin(y))/2 - (7215543768789811*sin(1361/1000))/4294967296 + 75000*sin
& H$ J5 [. V+ T. n+ y- z3 p0 `) |+ c4 h* S
(1361/1000)^3 - 75000*cos(1361/1000)*sin(1361/1000)*cos(y) - (6184732290414159*y*cos(1361/1000)*cos(y))/34359738368 -
; N, r8 z& Q. Y
$ k2 D8 A, T" E4 Y(6184732290414159*y*sin(1361/1000)*sin(y))/34359738368
& q# o3 D$ Q7 {: d3 M7 ~
1 r1 B3 K; n+ N% Z5 R
0 L8 n, x5 L* x0 A4 V. h9 T0 X>>  a=0; b=1;! \% r5 `# r9 V7 O3 O0 |9 w5 b
eps1=1e-8;eps2=1e-8;
' n) Q/ w! d8 q0 z N=300;
0 y+ f3 r. ]8 H- V8 G, `7 Af=@ (x) ((3359995*sin(1361/1000)*cos(x))/2 - (2999995*cos(1361/1000)*sin(x))/2 - (7215543768789811*sin(1361/1000))/4294967296 + 75000*sin(1361/1000)^3- 75000*cos(1361/1000)*sin(1361/1000)*cos(x) - (6184732290414159*x*cos(1361/1000)*cos(x))/34359738368 -(6184732290414159*x*sin(1361/1000)*sin(x))/34359738368);
  l0 ~7 M6 T6 D- a7 I  }Hfun=@Bisection;7 g  a) y; P. ?* S* c( F; W  y% U
[k,x,f_value]= feval(Hfun, f,a,b,eps1,eps2,N);1 f4 c: F' X) Y+ b/ _

/ b% D& y3 e+ F+ Z运行结果
7 J3 c1 a4 D5 V* ck ,a ,b ,x, f
  W" f* |* ]# h8 x$ m% S  1, 0.000000000, 1.000000000, 0.500000000,-352805.622314164,
4 U5 \( g7 ^; @" p. L  2, 0.000000000, 0.500000000, 0.250000000,-92968.331400711,$ r. q  N% u9 B1 z5 E
  3, 0.000000000, 0.250000000, 0.125000000,-4146.780462183,' H8 F  L, J4 i
  4, 0.000000000, 0.125000000, 0.062500000,29178.936016433,* N" h7 G8 C; F4 j
  5, 0.062500000, 0.125000000, 0.093750000,13458.364371693," l! o& M0 B$ z1 |: X# u9 o1 M
  6, 0.093750000, 0.125000000, 0.109375000,4890.124291139,3 s6 o) p5 n: g
  7, 0.109375000, 0.125000000, 0.117187500,430.092693351,0 N4 O2 b4 q/ G" G, k" x3 |$ P* o
  8, 0.117187500, 0.125000000, 0.121093750,-1843.759437944,. ^* ]0 v. o+ W# V1 r) t
  9, 0.117187500, 0.121093750, 0.119140625,-703.184646645,* K' c( r. O8 d3 A$ G8 U* G2 g
10, 0.117187500, 0.119140625, 0.118164063,-135.633470273,; j. L9 S4 n( |( u( V0 l7 c8 H
11, 0.117187500, 0.118164063, 0.117675781,147.457778641,
, A- X1 \1 C0 b( ?' v 12, 0.117675781, 0.118164063, 0.117919922,5.969190901,/ F% `! i1 {6 W. T- R; E8 P. L
13, 0.117919922, 0.118164063, 0.118041992,-64.817881139,
$ V9 m: |  x8 E2 A- Z 14, 0.117919922, 0.118041992, 0.117980957,-29.420780403,
$ m# D) Q6 [* l, M' d 15, 0.117919922, 0.117980957, 0.117950439,-11.724903562,: T5 |) H* ~2 C) g! R
16, 0.117919922, 0.117950439, 0.117935181,-2.877633532,) ~1 a( B) R: b' w* D% O  q
17, 0.117919922, 0.117935181, 0.117927551,1.545834384,
0 z# x: X. s& J: Y6 R* h* K/ P7 E3 q 18, 0.117927551, 0.117935181, 0.117931366,-0.665885649,
( j+ n0 g6 }; B9 A 19, 0.117927551, 0.117931366, 0.117929459,0.439977849,
) ]' L6 \2 c, B8 w6 F1 y! S 20, 0.117929459, 0.117931366, 0.117930412,-0.112953030,3 [" w# C0 O' t0 g
21, 0.117929459, 0.117930412, 0.117929935,0.163512627,
" y2 x  l4 c/ y6 j 22, 0.117929935, 0.117930412, 0.117930174,0.025279853,
3 W  D6 R7 d/ ~4 `% P3 E) _ 23, 0.117930174, 0.117930412, 0.117930293,-0.043836575,6 B$ V% S2 T0 A
24, 0.117930174, 0.117930293, 0.117930233,-0.009278357,2 P' B$ ~- ^( }# J0 x
25, 0.117930174, 0.117930233, 0.117930204,0.008000749,* _& B8 n9 F( f- F
26, 0.117930204, 0.117930233, 0.117930219,-0.000638804,
+ u, V3 E/ Z  l# o& ~/ e7 g 27, 0.117930204, 0.117930219, 0.117930211,0.003680972,# Q: @- G0 W$ L. J0 `6 T
>> 8 j! C5 l" R) x! K

' V0 a& a1 E" L7 d结论 当x=-1.361时  y=0.117930219
! D. `% P" B# s- M8 O5 \按此方法   可依次  得出150对x与y的值
. K, N: u3 b% ~
回复

使用道具 举报

发表于 2015-6-4 12:41:51 | 显示全部楼层
西交大学生??问邢老师
回复 支持 反对

使用道具 举报

发表于 2015-6-4 19:12:53 | 显示全部楼层
你这个是代数方程?
回复 支持 反对

使用道具 举报

发表于 2015-6-4 19:45:16 | 显示全部楼层
算法?

点评

二分法 还有一个程序 % Bisection.m function[k,x,f_value]=Bisection(f,a,b,eps1,eps2,N) fprintf('k ,a ,b ,x, f\n') for k=1:N x=(a+b)/2; f_value=f(x); fprintf( '%3d, %10.9f, %10.9   发表于 2015-6-5 08:36
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-5 08:37:59 | 显示全部楼层
cosxuan 发表于 2015-6-4 19:45
* R1 H& u& v( M0 |. i1 k算法?

5 n5 [0 A' r" }$ Z: H二分法        . x1 W# a/ n8 k# T
还有一个程序7 ]) P7 M8 \: W9 `0 c2 u
% Bisection.m* `  M+ G/ C, P9 l; X
function[k,x,f_value]=Bisection(f,a,b,eps1,eps2,N)0 m$ E* W  Y: |9 \4 H
fprintf('k ,a ,b ,x, f\n')
; [4 N$ t1 d% f( h# t; wfor k=1:N1 b5 h6 u7 u+ G& F
    x=(a+b)/2;$ ~0 H: o9 f- u; K' |1 h: r
    f_value=f(x);
9 ~1 J* @0 ~2 g  e1 c! P    fprintf( '%3d, %10.9f, %10.9f, %10.9f,%10.9f,\n'...( n4 _$ M) {# r
      ,k ,a ,b ,x, f_value)
: n7 p8 @- x1 w# S4 a  B    if abs(f_value)< eps1||0.5*(b-a)<eps2' W# v, M, W, x
    return  W5 J+ O8 {: e8 u- ^% t
    else3 {- k: Y0 v0 b1 o
        if f(x)*f(a)<0: P" }: {1 R9 w+ F
          b=x;
' @: y. u0 t8 ?6 o# @        else' y3 g" B1 y$ h' V
          a=x;
1 y" x6 P6 v- }+ s9 Y) ?* e- ~        end& u' Q1 w1 t; x& s3 D
        if k== N& |. X+ J6 Q; V! T7 V/ v  [1 ~
          warning ('算法超出最大迭代数!')   
6 S0 [/ e$ n9 h5 o; Y0 D6 c6 N        end; E: m" z+ ^6 |8 I
    end
' I7 E& p0 U2 Q( k* ?7 Z6 @
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-5 08:40:05 | 显示全部楼层
单车居士 发表于 2015-6-4 12:41
! i2 t  ?! z& c/ m# S/ k/ S6 b西交大学生??问邢老师

8 ?. M8 x& K2 Q这点问题   自己就可以解决
5 u4 k, W# e* i8 H, r) S
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-5 08:46:07 | 显示全部楼层
明月山河 发表于 2015-6-4 19:12 ! m7 Q' S% a, Y# \+ e( g4 a% M1 w0 U- M
你这个是代数方程?
2 [* C) C2 n+ O1 A1 f, o
是的呀      不过matlab写出来的方程  括号很多  不好看    0 x: p- }# ?0 p4 B& u
& i  ~( K$ h  N: a6 S  ]
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-6 18:55:35 | 显示全部楼层
明月山河 发表于 2015-6-4 19:12
8 Z+ F1 d& x5 n3 ]3 I7 ]你这个是代数方程?

# ~- h: E7 k0 l2 ?/ U2 x% U- M- V- c  J, w; [# Z  }0 s, c) A
是超越方程      6 }! W" B' }6 p1 w8 u; f& C
0 |- f5 k6 t6 t7 k! d+ f
* d0 e/ N( N) }5 j* I% H7 s& X
[chāo yuè fāng chéng
0 ^; L' z: U5 N8 m/ L2 U

点评

用数值方法,超越和代数差别不大。  发表于 2015-6-8 08:35
回复 支持 反对

使用道具 举报

发表于 2015-6-6 20:19:50 来自手机 | 显示全部楼层
提示: 作者被禁止或删除 内容自动屏蔽
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 注册会员

本版积分规则

小黑屋|手机版|Archiver|机械社区 ( 京ICP备10217105号-1,京ICP证050210号,浙公网安备33038202004372号 )

GMT+8, 2025-4-4 09:37 , Processed in 0.075919 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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