机械社区

 找回密码
 注册会员

QQ登录

只需一步,快速开始

搜索
查看: 2933|回复: 8

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

[复制链接]
发表于 2015-6-4 11:00:10 | 显示全部楼层 |阅读模式
螺杆转子刀具设计 数学模型) m# C; B  M- M" k: g
: ~9 \5 V5 A" ?
解如下方程: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
; X( u& z# J- @1 [求当x=-1.361时  y=?  
; D- T( U! l* o- x) T: {/ f6 Z% A" s>> syms x y
4 k$ [+ [- a5 R* E- P4 Q+ ]+ h>> 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)*cos0 T7 P0 c  I( M" c: ^
; J: b' G" Q3 k/ y
(x)+179999.4000005*y*sin(y)*sin(x)+75000*cos(y)*cos(x)*sin(x)
& o- b: Z9 c, ]' d7 a7 H# v4 h
) a% @  m; G' o0 O$ mf =
& p2 {9 i% u. s" F2 m7 X) I3 X( T* I8 A: [; S# `
(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)( A8 d% {8 }7 t" q( k

/ m; n: N+ G% {5 C/ w*sin(x) - (6184732290414159*y*cos(x)*cos(y))/34359738368 + (6184732290414159*y*sin(x)*sin(y))/34359738368% u! c* n/ r1 {  q( {, K
1 @; x: l, ^7 v3 X0 {
>>  subs(f, x, -1.361)
8 h: @: B6 ]( j" Z& h
  X2 a* o  l' \8 B7 s: ~ans =
- N5 q7 g) A: \4 B4 A$ c2 h9 R0 ~
. D" Y4 @# [; ]+ N1 m2 x: e(3359995*sin(1361/1000)*cos(y))/2 - (2999995*cos(1361/1000)*sin(y))/2 - (7215543768789811*sin(1361/1000))/4294967296 + 75000*sin
4 p6 I; \" X1 L9 Z* _  O. t* \6 Q5 t0 }
(1361/1000)^3 - 75000*cos(1361/1000)*sin(1361/1000)*cos(y) - (6184732290414159*y*cos(1361/1000)*cos(y))/34359738368 - 6 R% s; K# x, u1 Q1 M  d5 D' J

% [, a/ j  z4 X2 }! G/ a(6184732290414159*y*sin(1361/1000)*sin(y))/34359738368( q+ w( n3 {$ J& a2 N

5 G3 B5 p% x: R9 s! U# S) {% g& m3 M2 Q2 d* ]
>>  a=0; b=1;" i: D( r. g2 b' l$ ]7 N
eps1=1e-8;eps2=1e-8;+ X' D. G3 C% d# y% i! T
N=300;3 u9 E9 Q& a7 q1 ]; D
f=@ (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);
4 \8 j0 D+ g6 f  Q$ F2 g" j$ SHfun=@Bisection;
9 l4 c. \/ [8 Z% y[k,x,f_value]= feval(Hfun, f,a,b,eps1,eps2,N);
! J( C. N! F, x: r6 o2 ]" q
: `$ F) x4 s0 _运行结果  q( ^/ T; h" n/ z  P7 y
k ,a ,b ,x, f/ r. Y5 E  s5 V8 f. n8 p4 w$ o
  1, 0.000000000, 1.000000000, 0.500000000,-352805.622314164,
( z8 L) W1 h+ g7 h! r, Q2 A  2, 0.000000000, 0.500000000, 0.250000000,-92968.331400711,
/ t. a/ y9 n( u4 M9 ?7 J; u  3, 0.000000000, 0.250000000, 0.125000000,-4146.780462183,
( F3 c8 k# q( y/ m+ F- y  4, 0.000000000, 0.125000000, 0.062500000,29178.936016433,+ o3 ?5 V' F3 {. [( x1 ]$ U6 j
  5, 0.062500000, 0.125000000, 0.093750000,13458.364371693,/ w0 `) S8 q: p- \* s
  6, 0.093750000, 0.125000000, 0.109375000,4890.124291139,
( u$ c# \6 S* O2 c  k6 _9 v  7, 0.109375000, 0.125000000, 0.117187500,430.092693351,
& J8 k. g8 {1 J9 B  8, 0.117187500, 0.125000000, 0.121093750,-1843.759437944,1 l1 G, P; t! y- `
  9, 0.117187500, 0.121093750, 0.119140625,-703.184646645,9 Q$ }/ Q) c5 U) S( |' l
10, 0.117187500, 0.119140625, 0.118164063,-135.633470273,
. _) y: S6 }+ J* u5 | 11, 0.117187500, 0.118164063, 0.117675781,147.457778641,
: J5 ]+ L% _( X% h* O) E 12, 0.117675781, 0.118164063, 0.117919922,5.969190901,
' Y, M; G1 b' }, s7 x+ s/ H( n 13, 0.117919922, 0.118164063, 0.118041992,-64.817881139,1 D7 G: S) N/ X+ L( f1 Z; y
14, 0.117919922, 0.118041992, 0.117980957,-29.420780403,  ^( h. o* |# i' w8 n: X
15, 0.117919922, 0.117980957, 0.117950439,-11.724903562,8 D4 O& v! X" c4 l, j
16, 0.117919922, 0.117950439, 0.117935181,-2.877633532,9 T% |7 Q: }+ e# s* }
17, 0.117919922, 0.117935181, 0.117927551,1.545834384,% z, v) P: E! {6 M! d5 P
18, 0.117927551, 0.117935181, 0.117931366,-0.665885649,
: ?! p! |  Z! H# o 19, 0.117927551, 0.117931366, 0.117929459,0.439977849,
0 C+ [+ }% |8 R% C 20, 0.117929459, 0.117931366, 0.117930412,-0.112953030,( m, E/ L& M' W3 Q
21, 0.117929459, 0.117930412, 0.117929935,0.163512627,! T: d7 p+ Y+ Z! V( j. K
22, 0.117929935, 0.117930412, 0.117930174,0.025279853,
4 Y1 }; r1 J5 K, q0 a 23, 0.117930174, 0.117930412, 0.117930293,-0.043836575,- n5 |0 @4 `. |4 K
24, 0.117930174, 0.117930293, 0.117930233,-0.009278357,5 u! J$ X; P! B7 G5 W
25, 0.117930174, 0.117930233, 0.117930204,0.008000749,; E+ i: L( w/ j$ |# u
26, 0.117930204, 0.117930233, 0.117930219,-0.000638804,$ k% E6 J" @# q" r7 u) A, g
27, 0.117930204, 0.117930219, 0.117930211,0.003680972,/ M3 \0 C# V. t% J# I8 M% I0 @
>>
# _% I( P# v2 i: r
; P0 n& o# s% A' F结论 当x=-1.361时  y=0.117930219& a5 K, B* F( N1 P
按此方法   可依次  得出150对x与y的值6 L! p, T( T6 Z1 m- ?
回复

使用道具 举报

发表于 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 0 w& U9 j0 _# u( w  E( q
算法?
6 @6 x# u9 Q2 d/ v$ R+ ^
二分法        , W% _! J+ Q: ~4 G# X% }9 ~& p
还有一个程序
; Y) \, f1 `9 R$ p2 K( a0 n# K% Bisection.m. w/ X1 W& y$ U
function[k,x,f_value]=Bisection(f,a,b,eps1,eps2,N)' I$ x" k4 H: R* Z' j
fprintf('k ,a ,b ,x, f\n')
# z) E# J4 M- X4 P! `1 \4 L3 Qfor k=1:N
  F- x. P) k5 [    x=(a+b)/2;0 c% L2 K2 b0 Q- ~& m
    f_value=f(x);5 ~5 L& y! G% x( s
    fprintf( '%3d, %10.9f, %10.9f, %10.9f,%10.9f,\n'...
; n# i- l) v) }5 R      ,k ,a ,b ,x, f_value)
! |- I4 N# S2 O& I6 x- a0 u    if abs(f_value)< eps1||0.5*(b-a)<eps2
. |: z% h+ t# s' b1 T    return
6 R: ?2 w2 q+ q    else' _, e: s. p; g. W3 s
        if f(x)*f(a)<0+ H  p9 |- Y7 J9 T
          b=x;
/ z. |2 M# C8 J( O        else
) @! D( H# V4 N- A: A9 F: f          a=x;" A, ^* O5 f, h/ r2 D) Y; v
        end
$ f( E# ]8 t4 ~, ?: R        if k== N* E- q- H7 U1 D0 S$ H
          warning ('算法超出最大迭代数!')   
! {! H9 X4 r: u  S3 w1 N        end
) G6 [" y# q+ q! @    end5 o! B: S' z, U8 J7 s
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-5 08:40:05 | 显示全部楼层
单车居士 发表于 2015-6-4 12:41 9 P7 J( T- G7 ~7 B2 v8 m
西交大学生??问邢老师

$ y3 E- x* u. @3 X9 a这点问题   自己就可以解决
$ l2 P% f$ s( v. U9 `
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-5 08:46:07 | 显示全部楼层
明月山河 发表于 2015-6-4 19:12
! W8 F/ M! B! K- I4 P你这个是代数方程?
" G0 r) F: z& g/ ^, \" U6 b* e
是的呀      不过matlab写出来的方程  括号很多  不好看    8 N6 h. V5 b1 l6 g8 a+ S1 P1 r
) M  I) c9 V) c$ Y7 d7 [8 w' B
回复 支持 反对

使用道具 举报

 楼主| 发表于 2015-6-6 18:55:35 | 显示全部楼层
明月山河 发表于 2015-6-4 19:12
" Q- W6 y2 G7 I# O你这个是代数方程?

, ^, x% D  Y5 }$ [  J* M0 f0 Z- x$ ^6 j' q) x3 D2 \
是超越方程      9 j! ]4 \( s# E
+ T& y1 c6 |0 y) s; ]' @3 O$ `
' q$ X  G# V1 ], f( @; [, @! j
[chāo yuè fāng chéng
8 a: K, W, f/ i

点评

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

使用道具 举报

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

使用道具 举报

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

本版积分规则

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

GMT+8, 2025-4-2 13:40 , Processed in 0.070164 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

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