o****r 发帖数: 57 | 1 如图, 两个圆,半径都为L,两个圆心相距r. 已知r
问两个圆重合部位的面积。
试图用积分进行计算,积分区定位以一个圆心为中心的圆弧。
Matlab指令:
function y=area(l,L,r)
if l>=(L-r)
y=2*acos((l.^2+r.^2-L^2)/(2*l.*r))*l;
else
y=2*pi*l;
end
******************************************
L=10;
r=6;
quad(@area,0,L,[],[],L,r);%wrong result
235.2635
quad(@area,0,L-r,[],[],L,r)+quad(@area,L-r,L,[],[],L,r);%correct answer
195.9851
我想知道,为何必须分段积分才能得到正确的答案? function area应该如何定义才能
不用分段积分也能得到正确答案??
谢谢!!
补充说明:正确答案可以用简单的几何原理“直接扇形面积减三角形面积,再乘2”得出, 不一定要用这个积分方法。但我是在搞某个课题时候需要用 |
|
p******h 发帖数: 1783 | 2 把函数里面的if给去掉,换成逻辑判断
function y=area(l,L,r)
y=2*acos((l.^2+r.^2-L.^2)./(2*l.*r)).*l .*(l>=(L-r)) + 2*pi*l .*(l<(L-r));
end
不过算出来的结果有一项较小的虚数项,再在function里面取实数
就行了 |
|
h**********c 发帖数: 4120 | 3 如何处理这样的精度?
double r1 = 1.0*2.0 + 1.0*2.0 + 1.0 *2.0;
double r10 = r1/sqrt(1.0+1.0+1.0)/sqrt(2.0*2.0+2.0*2.0+2.0*2.0);
VC++ 2010 express, Vista 64 结果
r10 = 1.0000000000000002;
求acos(r10)的时候溢出,用float也一样戴一个尾巴。
这样的运算怎样避免溢出?
非express 版也有这个问题吗?
没有试别的compler. |
|
h**********c 发帖数: 4120 | 4 THNKS, the article is a little bit long.
Has someone tried to rewrite acos/atan sth? |
|
d****e 发帖数: 12 | 5 Stat@OSU和Math@Gatech,都是PhD
这两个钱差不多,都是24000左右
本人以后更想去业界发展
OSU的好处是可以拿个硕士走人找工作,但地理位置没有亚特兰大好
Gatech的话进去读计算数学。数学系与IE系有合办ACO的项目,但到时候不知道有没有
机会转到那个项目去读optimization。但是如果转不成,就读计算数学,也能接受,但
好像读计算对计算机编程要求比较高,本人对计算机不是很感兴趣。
我希望能够去个找工作前景比较好的地方,希望各位能给我指点下
谢谢~~~ |
|
a*********m 发帖数: 26 | 6 Dear friends,
Welcome aboard!
Each year after match day, the Student/Resident/Fellow Committee of ACAP
will host the "Welcome and Career Development Seminar". It will be a good
time for us to know and network with each other, and learn from our senior
residents/ attendings. This year, the seminar will be hosted at Flushing
Hospital on 3/24/2013. The topics will focus on "How to be a good intern"
and "GI/ Cardiology/ Hem-Onc Fellowship Application".
The following is the agenda for 3/24/2013 "Welc... 阅读全帖 |
|
a*********m 发帖数: 26 | 7 Dear friends,
Welcome aboard!
Each year after match day, the Student/Resident/Fellow Committee of ACAP
will host the "Welcome and Career Development Seminar". It will be a good
time for us to know and network with each other, and learn from our senior
residents/ attendings. This year, the seminar will be hosted at Flushing
Hospital on 3/24/2013. The topics will focus on "How to be a good intern"
and "GI/ Cardiology/ Hem-Onc Fellowship Application".
The following is the agenda for 3/24/2013 "Welc... 阅读全帖 |
|
a*********m 发帖数: 26 | 8 Dear friends,
Congratulations to all the CMGs matched this year!
This is a reminder for the "Welcome and Career Development Seminar" on 3/24/
2013 at Flushing Hospital, New York. Each year after match day, the Student/
Resident/Fellow Committee of ACAP will host the "Welcome and Career
Development Seminar". It will be a good time for us to know and network with
each other, and learn from our senior
residents/ attendings. This year, the seminar will be hosted at Flushing
Hospital on 3/24/2013. Th... 阅读全帖 |
|
u*******s 发帖数: 688 | 9 受到某位大神指点,开始有意无意关注Sara White。虽然没有读过她的长篇大作(汗
),但是经常看她的博客更新,觉得内容很有帮助。这里就是一篇指导药剂师(和其他
所有职业)增加自己在机构里的可视度的小短文。作为培养leadership的第一步,
visibility很重要。文中列举的都是简单可行的小tips。比如开会的时候大胆发言(可
以先从提问开始培养胆量),多为组织volunteer或者ask for highly visible
project,多参加学习机会,和多承认他人的成就。
下面是全文给大家参考。
How Visible Are You?
By Sara White posted Sun, Jan 25, 2015 09:33 AM
To be an effective leader for your area you must increase your visibility/
raise your profile at work. The Mind Tools Club offers the following advice.
Most people feel un... 阅读全帖 |
|
e**********n 发帖数: 359 | 10 对头,仔细想一下还是挺有意思的,不过还是不用解测地线方程,应为光线的方程不过是
ds = 0
有质量的东西才用测地线。
如果光线在x,z平面内,此方程变成
(z')^2+(x')^2 [1+Acos(w(z-t))] = 1
x',z'分别是对t的导数。用x'= 1近似,按常规应该够了。
z’=+/- A cos(w(z-t))
一开始取正号,w(z-t) = pi/2 时开始取负号直到z又等于0,被另一艘飞船收到。也可
继续传播下去,到w(z-t) = 3pi/2 时再改变一次符号。
计算一下近似的光程就可以得到往返的周期。发射和接收的时
候z'=A,所以镜子因该有一个偏转角 +\-arctan(A),符号根据路径选择。
具体路径不用求出,如果一定要算,只能把z当小量做微扰。
z(t)=Ae^{-A sin (wt)/w} \int_0^t e^{A sin(ws)/w} sin(ws) ds
大概就是这样了,好久不做习题引力波波速是光速都快忘了。
其实我原来说的也没错,只要把光线投影到x-y平面上,定义一个等效的折射率,结果是一样
的。只是这样不能得到发射和接收的角度。 |
|
b***k 发帖数: 2673 | 11
F(
Var
pi=2*acos(0.0) or pi=4*atan(1.0) |
|
s******y 发帖数: 14 | 12
现在被录的是Mathmatical Sci.不过具体方向还没定,他说是后来细分成传统应数、
PhD in ACO或Financial Math.都有可能;anyway, 谢谢你 |
|
w**d 发帖数: 2334 | 13 All following are in the paper:
D. Gottlieb, M.Y. Hussaini and S.A. Orzag, 1984
Theory and application of spectral methods.
(SIAM-CBMS, Philadelphia, PA, 1984) pp 1-94.
Given a smooth function f(x), to get its upto n-degree
polynomials interpolation.
Assume domain : [-1,1],
take Xj= cos(j*pi/n) as the interpolation points.
take Tj(x)=cos( j* acos(x)) as the n+1 basis.
Note Tj is just j-degree polynomial.
let f_n = A0*T0 + ... + An*Tn, then by definition of interpolation
A0*T0(X0) + ... |
|
A*******s 发帖数: 9638 | 14 I think cutting cost is a consensus for everyone, but how?
ACO just sounded too good to be true.
money! |
|
m******h 发帖数: 297 | 15 来自主题: Medicalpractice版 - 职业指导 请问楼主关于Affordable care的理解有没有资料的link?觉得lz是想谈accountable
care organization (ACO)。
"Obamacare就是要改变这个状况,美国没钱了,就得过紧日子。 所以Afford care act
(ACA)要做的是改变fee for service,代之为capitation或者salary。 比方说一个
医院一年得到联邦补助1个固定的数目,让你cover50万人口,花多了自己负责,花少了
就是利润。 医院得到了这些钱,它就不能像以前那样接受病人越多,做检查越多,利
润越高。 相反,它会尽一切力量减少ER visit, 比方说帮助Nursing home提供基本的
医疗服务来减少ER/hospital visit; 它会帮助hospice program做生意,来减少病入
膏肓的病人不必要的医疗开支; 它更会帮助urgent care这样的business来shunt它的
病人,减少ER的负担。。。。。。" |
|
|
t****t 发帖数: 158 | 17 来自主题: Medicalpractice版 - 职业指导 中國古代有本書叫 黃帝內經,說過,不治已病治未病,不治已亂治未亂。上醫治國,
中醫治人,下醫治病。
扯遠了,“ACO和BPCI强调的是各科之间的合作,觉得是一个进步啊, too
到底obamacare糟糕在哪呢? |
|
A*******s 发帖数: 9638 | 18 来自主题: Medicalpractice版 - 职业指导 你说的ACO并不需要obamacare才能实行。
obamacare说白了就是socialism。 capitation也好, 工资也好,目的是减少开支,但
客观上就是大锅饭,干多干少一个样。
医生没有了积极性,对病人是好还是坏? |
|
A*******s 发帖数: 9638 | 19 Let's see what happens next. |
|
A*******s 发帖数: 9638 | 20 When Obama finishes his term.
MGH,
★ 发自iPhone App: ChineseWeb 7.8 |
|
A*******s 发帖数: 9638 | 21 They try to dictate physician practice. But do they know why some physicians
write more MRIs? |
|
a**e 发帖数: 5094 | 22 来自主题: Medicalpractice版 - 感觉很不好 其实就是比看谁更流氓。保险公司敢不给钱,我们学校自己卖保险,从HMO到PPO,连
malpractice insurance都卖,而且公司设在什么virgin island。你搞ACO,我们有自
己的PCP chain store。 |
|
p*********e 发帖数: 32207 | 23 GT1@1994:恐怖Dauer
一个时代的结束也往往正是宣告另一个时代的到来。94年开始BPR组织了一
个世界性的GT系列赛,其中最高的一个组别正是本文的主角——GT1。同时
,ACO也在Le Mans大赛中增加了相应的组别。
首届BPR GT最终成了又一次德法战争。总共8站比赛,规则并不完全相同,4
小时,3小时,1000km以及500km的都有。Porsche以964平台的911 Turbo
S为基础改装了参赛的911 Turbo S LM,法国人则有装备Renault发动机的
Venturi 600 LM,结果双方各获得3站优胜,但Porsche还是排到了积分榜
的首位。另外两站的第一名则被改装了Turbo的Ferrari F40获得。
同年的Le Mans竞争则更加激烈。C1组取消后,新的C组比赛实际上更接近于
原来的低一级别的C2组,但毕竟是车身重量和空气动力学方面的优势都十分
明显的运动原型车,人们并不知道新生的GT1赛车实力是否能够向C组赛车挑
战。C组的主要参赛车辆有原C2组中的大牛,Toyota的C-V LM系列的94最新
款,还有都装备了Porsche发动机的Kr |
|
|
k*********r 发帖数: 758 | 25 ~{!>~} ~{TZ~} allegretto (hollywood sucks) ~{5D4sWwVPLa5=~}: ~{!?~}
~{6
~{PB0f5D~}lolita, ~{>]K5A=F,A,O8=Z6<:\O`KF!#6x:sCf5DDG8vSP~}jeremy irons~{!#~}
~{:\Ff9V5DJGN*J2C4C;SP~}A clockwork orange~{#?K5@4G0P)LlE
student center~{3T79#,N^RbVP7"OVG=IO5D5gJS2%5D~}MTV~{U}:CJG8v@V6STZQ]~}
ACO~{VP5DR;P)F,6N#,2;V*5@JGDD8v@V6S!#~} |
|