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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求教：Matlab利用最小二乘法进行参数估计（谐波方程）
51/1返回列表
查看: 1779  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fudange

木虫 (小有名气)

农民

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

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

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
路漫漫其修远兮，吾将上下而求索！
1楼 2015-01-21 08:43:07
已阅   回复此楼   关注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)的投影谱，原理类似二维傅里叶变换，但是反演核矩阵不同。这样做可以免除阶数的猜测，且运算属于凸规划，理论上保证收敛。但问题变为求解欠定线性方程，可能需要用到压缩感知算法。就好像你对时间域数据进行处理，为什么采用傅里叶变换而不采用多阶拟合？因为傅里叶变换可以给出全频率域信息，而多阶拟合只能在明确知道阶数以及大致知道所有阶的信息时才能给出较可靠的结果。
这个才是真的有用
赞一下(1人) 回复此楼
5楼2015-01-22 08:55:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fudange 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 材料080500调剂求收留 +8 一颗meteor 2026-03-13 8/400 2026-03-20 20:00 by zhukairuo
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 085600材料与化工 +6 安全上岸！ 2026-03-16 6/300 2026-03-20 15:07 by lbsjt
[考研] 295材料求调剂，一志愿武汉理工085601专硕 +4 Charlieyq 2026-03-19 4/200 2026-03-20 14:26 by 无懈可击111
[考研] 286分人工智能专业请求调剂愿意跨考！ +3 lemonzzn 2026-03-17 4/200 2026-03-20 11:04 by lemonzzn
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 294求调剂材料与化工专硕 +14 陌の森林 2026-03-18 14/700 2026-03-19 22:38 by 学员8dgXkO
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 6/300 2026-03-19 17:18 by fei626-918
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 一志愿西南交大，求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 材料，纺织，生物（0856、0710），化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 中科院材料273求调剂 +4 yzydy 2026-03-15 4/200 2026-03-16 15:59 by Gaodh_82
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
[考研] 本科南京大学一志愿川大药学327 +3 麦田耕者 2026-03-14 3/150 2026-03-14 20:04 by 外星文明
信息提示
关闭
请填处理意见
关闭