Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 1802  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fudange

木虫 (小有名气)

农民

[求助] 求教:Matlab利用最小二乘法进行参数估计(谐波方程) 已有4人参与

请教众大神们,小弟在利用谐波法模拟土壤温度的时候遇到一些问题。以下公式是通过变量分离法得出热传导方程的解析解:

其中,T(z,t)是实测的温度,T(z)是日均温度(上面有横线),z是水分深度,w是地球自转的角速度,为已知常数,t是记录温度的时间;n为谐波阶数,Ai和φi分别表示第i阶谐波的振幅和相位。公式实际上是实测温度T(z,t)随时间t变化的函数,现在需要用最小二乘法估计参数Ai,Bi以及φ的之。由于不确定函数的阶数,不知道该怎么编写命令代码,希望各路大神尽显神通,帮小弟解决燃眉之急!在此谢过各路大侠!

求教:Matlab利用最小二乘法进行参数估计(谐波方程)
公式.png
回复此楼

» 猜你喜欢

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

路漫漫其修远兮,吾将上下而求索!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你实测了多少个点,你就可以最小二乘多少个参数,所以N个实验数据最多可以确定到N/3阶。
3楼2015-01-21 12:08:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

wdlgezi

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
只是阶数不确定问题的话,你大可以给定一个。各次谐波整体呈衰减走势,n越大精确度越高

[ 发自小木虫客户端 ]
敢问路在何方
2楼2015-01-21 08:50:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

戴钢盔的猪头

木虫 (知名作家)

【答案】应助回帖

感谢参与,应助指数 +1
以上说的都对,但都不是普适方法。可以将后面非线性部分作为基展开,通过正则化反演得到A关于离散基exp(-Bz)*sin(iwt+phi-Bz)的投影谱,原理类似二维傅里叶变换,但是反演核矩阵不同。这样做可以免除阶数的猜测,且运算属于凸规划,理论上保证收敛。但问题变为求解欠定线性方程,可能需要用到压缩感知算法。就好像你对时间域数据进行处理,为什么采用傅里叶变换而不采用多阶拟合?因为傅里叶变换可以给出全频率域信息,而多阶拟合只能在明确知道阶数以及大致知道所有阶的信息时才能给出较可靠的结果。

[ 发自手机版 http://muchong.com/3g ]
4楼2015-01-21 12:34:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaodong1023

金虫 (初入文坛)

【答案】应助回帖

感谢参与,应助指数 +1
上说的都对,但都不是普适方法。可以将后面非线性部分作为基展开,通过正则化反演得到A关于离散基exp(-Bz)*sin(iwt+phi-Bz)的投影谱,原理类似二维傅里叶变换,但是反演核矩阵不同。这样做可以免除阶数的猜测,且运算属于凸规划,理论上保证收敛。但问题变为求解欠定线性方程,可能需要用到压缩感知算法。就好像你对时间域数据进行处理,为什么采用傅里叶变换而不采用多阶拟合?因为傅里叶变换可以给出全频率域信息,而多阶拟合只能在明确知道阶数以及大致知道所有阶的信息时才能给出较可靠的结果。
这个才是真的有用
5楼2015-01-22 08:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂,一志愿华南理工大学071005 +25 ioodiiij 2026-04-08 25/1250 2026-04-14 15:59 by zs92450
[考研] 245求调剂 +6 冰糖橘?汽水 2026-04-13 10/500 2026-04-14 10:49 by jyl0317
[考研] 085404 298分求调剂 +11 呼啦呼啦呼呼呼 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[基金申请] 2026 WR青拔 +3 冬日阳光CAS 2026-04-09 6/300 2026-04-13 18:40 by liuchb715
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 322求调剂 +6 123安康 2026-04-12 13/650 2026-04-12 15:51 by 123安康
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 299求调剂 +8 ZVVZ13 2026-04-08 8/400 2026-04-12 00:40 by 蓝云思雨
[考研] 化工调剂求导师收留!一志愿失利,踏实肯干,有植物提取科研经历 +20 yzyzx 2026-04-09 21/1050 2026-04-12 00:12 by 小小小小啦啦啦
[考研] 0859,337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 11:34 by caotw2020
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 机械专硕270求调剂,接受跨专业 +12 老师看看我吧aba 2026-04-09 14/700 2026-04-11 10:21 by laoshidan
[考研] 中药学调剂 初试324 +4 洋甘菊、 2026-04-10 6/300 2026-04-11 09:41 by gong120082
[考研] 085402通信工程调剂,有4项学科竞赛国奖(电赛国二),硕士研究生调剂自荐信。 +5 m永o不v言o弃m 2026-04-09 5/250 2026-04-11 09:33 by zhq0425
[考研] 一志愿东北大学控制工程085406数二英二385,求调剂 +8 Ezra_Zhang 2026-04-09 8/400 2026-04-11 09:15 by 猪会飞
[考研] 342电子信息专硕求调剂 +9 你让我怎么荔枝 2026-04-10 10/500 2026-04-11 08:33 by zhq0425
[考研] 调剂 +19 小张ZA 2026-04-10 20/1000 2026-04-10 22:08 by 猪会飞
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
信息提示
请填处理意见