24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yuding_sun

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 368.3
帖子: 40
在线: 11.6小时
虫号: 3159600
注册: 2014-04-24
性别: GG
专业: 光子与光电子器件

[求助] 做光热热声实验得到的数据用1stopt进行非线性回归分析存在一些问题

图片是理论公式，求D的值
数据如下
f(Hz) phi
1 103.1
11 56.68
21 26.05
31 6.01
41 10.96
51 -21.8
61 -30.97
71 -37.41
81 -40.6
91 -45
101 -50.2
111 -50.7
121 -53.98
131 -55.01
141 -55.01
151 -59.04
161 -60.95
171 -60.95
181 -66.04
191 -63.44
201 -63.44
211 -69.45
221 -69.45
231 -66.8
241 -66.8
251 -66.8
261 -71.57
271 -71.57
281 -71.57
291 -71.57
301 -71.57
我用1stopt15做的代码如下
Title "f";
Parameter D;
Variable f,phi;
Function
tan(phi)=(1-3/(2*sqrt(3.14*f/D)*2)*(sin(sqrt(3.14*f/D)*2)+1/2*(exp(sqrt(3.14*f/D)*2)-exp(-sqrt(3.14*f/D)*2)))/(cos(sqrt(3.14*f/D)*2)+1/2*(exp(sqrt(3.14*f/D)*2)+exp(-sqrt(3.14*f/D)*2))))/(3/(2*sqrt(3.14*f/D)*2)*(sin(sqrt(3.14*f/D)*2)-1/2*(exp(sqrt(3.14*f/D)*2)-exp(-sqrt(3.14*f/D)*2)))/(cos(sqrt(3.14*f/D)*2)+1/2*(exp(sqrt(3.14*f/D)*2)+exp(-sqrt(3.14*f/D)*2))));
Data;
//f phi
1 103.1
11 56.68
21 26.05
31 6.01
41 10.96
51 -21.8
61 -30.97
71 -37.41
81 -40.6
91 -45
101 -50.2
111 -50.7
121 -53.98
131 -55.01
141 -55.01
151 -59.04
161 -60.95
171 -60.95
181 -66.04
191 -63.44
201 -63.44
211 -69.45
221 -69.45
231 -66.8
241 -66.8
251 -66.8
261 -71.57
271 -71.57
281 -71.57
291 -71.57
301 -71.57
算出结果D为90多，用教授的matlab程序跑出来是107
1stopt出的图像很奇怪
用SPSS做的更扯，D算的是10。
请各位高手帮忙处理下，用1stopt和matlab都行，要能出图，算到110附近就好。

捕获.JPG

1.JPG

回复此楼

» 猜你喜欢

求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有274人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有19人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

非线性回归如何操作（应力应变曲线数据如何进行非线性回归）求指导！已经有5人回复

1楼2014-04-24 10:57:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yuding_sun

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 368.3
帖子: 40
在线: 11.6小时
虫号: 3159600
注册: 2014-04-24
性别: GG
专业: 光子与光电子器件

引用回帖:

2楼: Originally posted by dingd at 2014-04-24 13:31:29
Title "f";
Parameter D;
Variable f,phi;
Function phi=180/pi*Arctan((1-3/(2*sqrt(3.14*f/D)*2)*(sin(sqrt(3.14*f/D)*2)+1/2*(exp(sqrt(3.14*f/D)*2)-exp(-sqrt(3.14*f/D)*2)))/(cos(sqrt(3.14*f ...

我的程序是哪里出问题了呢

赞一下

回复此楼

3楼2014-04-24 23:32:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yuding_sun

捐助贵宾 (初入文坛)

应助: 0 (幼儿园)
金币: 368.3
帖子: 40
在线: 11.6小时
虫号: 3159600
注册: 2014-04-24
性别: GG
专业: 光子与光电子器件

谢谢，只是不知道为什么我代码的问题在哪里

赞一下

回复此楼

4楼2014-04-24 23:34:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 yuding_sun 的主题更新

返回列表