24小时热门版块排行榜    

查看: 2062  |  回复: 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的回帖

wdlgezi

金虫 (正式写手)

【答案】应助回帖

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

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

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你实测了多少个点,你就可以最小二乘多少个参数,所以N个实验数据最多可以确定到N/3阶。
3楼2015-01-21 12:08:32
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 fudange 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-25 01:47 by BZKMTicpDhFj
[基金申请] 明天应该可查了!? +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +7 scalable 2026-08-24 7/350 2026-08-24 23:01 by anata1209
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 0XLacIJUOj8D 2026-08-24 3/150 2026-08-24 21:59 by BZKMTicpDhFj
[基金申请] 人气不行了 +10 fansofjerry 2026-08-21 10/500 2026-08-24 21:03 by zhanghaozhu
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见