18楼: Originally posted by wurongjun at 2015-07-22 20:34:19
定子面积程序:
function dingz
y=4*rombg(@fdz,0,pi/2,1e-2);
disp(num2str(y,6))
function ds1=fdz(t)
z1=6;
z2=z1-1;
R=48.78;
r=8.13;
Q=R-r;
e=7.05;
syms tau;
k=12.6;    %re
f=r/e;
a=sin(z1* ...

4楼: Originally posted by cooooldog at 2015-07-18 08:32:42
下面的曲线外轮廓所包围的面积用数值近似计算得到大约是
7936.86...

2楼: Originally posted by wurongjun at 2015-07-17 11:09:33
以前回答过这个类似问题!
quad函数里面调用的函数你没有写好,写成表达式啦!
要写成function形式的!
关键是你的问题里面其他参数!
这个新版的Matlab里面已经可以解决啦!
搜索匿名函数的用法!完全可以解决你的问 ...

3楼: Originally posted by cooooldog at 2015-07-17 11:18:14
1. 问问题自己先把问题尽可能描述清清楚楚很重要;

2. 定子曲线方程是一条闭合的曲线


但是转子曲线方程无法在0~2 pi之内闭合


3. 你自己的代码除了有一些已知量的赋值可参考之外,其它几乎毫无意义,而且t ...

5楼: Originally posted by 丫头丫头2014 at 2015-07-20 10:50:01
您好，我已经尽自己最大努力，尽量把问题说清楚了，写成word形式，麻烦您再看一下，如果哪块不清楚，我再解释，这个问题困扰了我很久，不知道该怎么解决，先谢谢您啦...

7楼: Originally posted by cooooldog at 2015-07-20 14:01:43
主要不清楚的地方是，您到底求定子，转子面积，还是求曲线的外轮廓包围的面积？

正如前面已经给出的两个曲线的图片一样

1. 定子对应的曲线是封闭的，但是有两个环绕在曲线内的圈，这个是\tau 在0-2 \Pi 正好 ...

8楼: Originally posted by 丫头丫头2014 at 2015-07-20 14:35:49
您之前给出的两个图形是不正确的 我在文本中附上了正确的图片，两个图形都是要0-2pi内闭合曲线围成的面积，定子方程在0-2pi是差一点点闭合的，这个是为什么，我也不是特别清楚，转子的方程，是依赖定子方程产生的， ...

8楼: Originally posted by 丫头丫头2014 at 2015-07-20 14:35:49
您之前给出的两个图形是不正确的 我在文本中附上了正确的图片，两个图形都是要0-2pi内闭合曲线围成的面积，定子方程在0-2pi是差一点点闭合的，这个是为什么，我也不是特别清楚，转子的方程，是依赖定子方程产生的， ...

ClearAll["Global`*"];
R = 4878/100;
r = 813/100;
z1 = R/r;
z2 = z1 - 1;
e = 705/100;
f = r/e;
re = 126/10;
\[Theta] = ArcTan[Sin[z1 \[Tau]]/(f + Cos[z1 \[Tau]])] - \[Tau];
\[CurlyPhi] = ArcSin[f Sin[\[Theta] + \[Tau]]] - \[Theta];
\[Psi] = z1/(z1 - 1) \[CurlyPhi];
curve01 = {(R - r) Sin[\[Tau]] + e Sin[z2 \[Tau]] -
     re Sin[\[Theta]], (R - r) Cos[\[Tau]] - e Cos[z2 \[Tau]] +
     re Cos[\[Theta]]} // FullSimplify;
curve02 = {curve01[[1]] Cos[\[CurlyPhi] - \[Psi]] -
     curve01[[2]] Sin[\[CurlyPhi] - \[Psi]] - e Sin[\[Psi]],
    curve01[[1]] Sin[\[CurlyPhi] - \[Psi]] +
     curve01[[2]] Cos[\[CurlyPhi] - \[Psi]] - e Cos[\[Psi]]} //
   Simplify;

base = ParametricPlot[Evaluate[curve01], {\[Tau], 0, 2 \[Pi]},
   Exclusions -> None, MaxRecursion -> 15, PlotPoints -> 1500,
   AxesStyle -> Arrowheads[.035],
   AxesLabel -> {Style[#, Black, Bold, 21 ,
        FontFamily -> "Times New Roman", Italic] &@"x",
     Style[#, Black, Thick, 21 , FontFamily -> "Times New Roman",
        Italic] &@"y"},
   TicksStyle -> Directive[Black, 14, FontFamily -> "Arial"],
   PlotRangePadding -> Scaled[.1], ImageSize -> 450] /.
  Line[x_] :> {Blue, Thickness[.001], Line[x]}


ParametricPlot[Evaluate[curve02], {\[Tau], 0, 2 \[Pi]},
  Exclusions -> None, MaxRecursion -> 15, PlotPoints -> 1500,
  AxesStyle -> Arrowheads[.035],
  AxesLabel -> {Style[#, Black, Bold, 21 ,
       FontFamily -> "Times New Roman", Italic] &@"x",
    Style[#, Black, Thick, 21 , FontFamily -> "Times New Roman",
       Italic] &@"y"},
  TicksStyle -> Directive[Black, 14, FontFamily -> "Arial"],
  PlotRangePadding -> Scaled[.1], ImageSize -> 450] /.
Line[x_] :> {Red, Thickness[.001], Line[x]}

NIntegrate[
   curve01[[2]] D[curve01[[1]], \[Tau]] -
    curve01[[1]] D[curve01[[2]], \[Tau]], {\[Tau], 0, 2 Pi}]/2 //
NumberForm[#, 10] &

NIntegrate[
   curve02[[2]] D[curve02[[1]], \[Tau]] -
    curve02[[1]] D[curve02[[2]], \[Tau]], {\[Tau],
    0, \[Tau] /. FindRoot[curve02[[2]] == 783/20, {\[Tau], 1.6 Pi}]},
   PrecisionGoal -> 6, Method -> "LocalAdaptive"]/2 //
NumberForm[#, 10] &