版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

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

回复此楼

» 猜你喜欢

求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有138人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收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的回帖

回帖置顶 ( 共有1个 )

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.5小时
虫号: 291104
注册: 2006-10-28

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
fegg7502: 金币+1, 3ks 2014-04-25 08:23:20
yuding_sun: 回帖置顶 2014-04-25 09:03:33
yuding_sun(dbb627代发): 金币+30 2014-04-25 09:50:34

CODE:

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/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

均方差(RMSE):4.73865362084477
残差平方和(SSE):696.099982288702
相关系数(R): 0.994572026107018
相关系数之平方(R^2): 0.989173515114618
确定系数(DC): 0.986599889046332
卡方系数(Chi-Square): 68.4609625135357

参数最佳估算
---------- -------------
d 111.245186098327

f1.jpg

赞一下(1人)

回复此楼

2楼2014-04-24 13:31:29

已阅回复此楼关注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的回帖