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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 光机热集成分析时采用Zernike多项式对变形镜面拟合的问题！（最小二乘法）
51/1返回列表
查看: 1238  |  回复: 4
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者老铁双击999将赠送您 20 个金币

老铁双击999

新虫 (小有名气)

[求助] 光机热集成分析时采用Zernike多项式对变形镜面拟合的问题！（最小二乘法）已有1人参与

如题：


最近在做镜头表面的光机热耦合分析，先在workbench中获得镜头某一面的变形数据，然后在Matlab中编写了Zernike多项式拟合程序，采用的是最小二乘法，但是获得的Fringe Zernike 37项系数非常大，将37项系数带入到最小二乘方程中却能很好的吻合，我看别人采用的正交方法得到的37项Zernike系数都非常小。求虫友查看问题所在。。。。另外哪位朋友有现成的光机热集成分析程序愿意分享一下，十分感激！！！


附带自己做的matlab采用最小二乘法的Zernike拟合程序

下面是求解程序：

程序一：求解程序

clear all;clc
format long;
S= xlsread('C:\Users\zzz\Desktop\A面.xlsx');
x=S(:,2);%拟合点X坐标
y=S(:,3);%拟合点Y坐标
z=S(:,4);%拟合点Z坐标
DX=S(:,5);%拟合点X坐标误差
DY=S(:,6);%拟合点Y坐标误差
DZ=S(:,7);%拟合点Z坐标误差
W=(DZ+z);%拟合点实际位置
%W=0.00001*(DZ+z);%拟合点实际位置
for i=1:1:length(x) %一共有多少个点
Z1(i)=(x(i)+DX(i))^0;
Z2(i)=2*(x(i)+DX(i));
Z3(i)=2*(y(i)+DY(i));
Z4(i)=1.732*(2*(x(i)+DX(i))^2+2*((y(i)+DY(i))+DY(i))^2-1);
Z5(i)=2*sqrt(6)*(x(i)+DX(i))*((y(i)+DY(i))+DY(i));
Z6(i)=sqrt(6)*((x(i)+DX(i))^2-((y(i)+DY(i))+DY(i))^2);
Z7(i)=sqrt(8)*(3*(x(i)+DX(i))^2+3*((y(i)+DY(i))+DY(i))^2-2)*((y(i)+DY(i))+DY(i));
Z8(i)=sqrt(8)*(3*(x(i)+DX(i))^2+3*((y(i)+DY(i))+DY(i))^2-2)*(x(i)+DX(i));
Z9(i)=sqrt(8)*(3*(x(i)+DX(i))^2-((y(i)+DY(i))+DY(i))^2)*((y(i)+DY(i))+DY(i));
Z10(i)=sqrt(8)*(3*(x(i)+DX(i))^2-((y(i)+DY(i))+DY(i))^2)*(x(i)+DX(i));
Z11(i)=sqrt(5)*(6*((x(i)+DX(i))^2+((y(i)+DY(i))+DY(i))^2)^2-6*((x(i)+DX(i))^2+(y(i)+DY(i))^2)+1);
Z12(i)=sqrt(10)*(4*(x(i)+DX(i))^2+4*(y(i)+DY(i))^2-3)*((x(i)+DX(i))^2-(y(i)+DY(i))^2);
Z13(i)=2*sqrt(10)*(4*(x(i)+DX(i))^2+4*(y(i)+DY(i))^2-3)*(x(i)+DX(i))*(y(i)+DY(i));
Z14(i)=sqrt(10)*((x(i)+DX(i))^4-6*(x(i)+DX(i))^2*(y(i)+DY(i))^2+(y(i)+DY(i))^4);
Z15(i)=4*sqrt(10)*((x(i)+DX(i))^2-(y(i)+DY(i))^2)*(x(i)+DX(i))*(y(i)+DY(i));
Z16(i)=sqrt(12)*(10*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^2-12*((x(i)+DX(i))^2+(y(i)+DY(i))^2)+3)*(x(i)+DX(i));
Z17(i)=sqrt(12)*(10*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^2-12*((x(i)+DX(i))^2+(y(i)+DY(i))^2)+3)*(y(i)+DY(i));
Z18(i)=sqrt(12)*(5*((x(i)+DX(i))^2+(y(i)+DY(i))^2)-4)*((x(i)+DX(i))^2-3*(y(i)+DY(i))^2)*(x(i)+DX(i));
Z19(i)=sqrt(12)*(5*((x(i)+DX(i))^2+(y(i)+DY(i))^2)-4)*((x(i)+DX(i))^2-3*(y(i)+DY(i))^2)*(y(i)+DY(i));
Z20(i)=sqrt(12)*((x(i)+DX(i))^5-10*(x(i)+DX(i))^3*(y(i)+DY(i))^2-(x(i)+DX(i))*(y(i)+DY(i))^4);
Z21(i)=sqrt(12)*(5*(x(i)+DX(i))^4*(y(i)+DY(i))-10*(x(i)+DX(i))^2*(y(i)+DY(i))^3+(y(i)+DY(i))^5);
Z22(i)=sqrt(7)*(20*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^3-30*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^2+12*((x(i)+DX(i))^2+(y(i)+DY(i))^2)-1);
Z23(i)=sqrt(14)*(30*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^2-40*((x(i)+DX(i))^2+(y(i)+DY(i))^2)+12)*(x(i)+DX(i))*(y(i)+DY(i));
Z24(i)=sqrt(14)*(15*((x(i)+DX(i))^2+(y(i)+DY(i))^2)^2-20*((x(i)+DX(i))^2+(y(i)+DY(i))^2)+6)*((x(i)+DX(i))^2+(y(i)+DY(i))^2);
Z25(i)=sqrt(14)*(24*((x(i)+DX(i))^2+(y(i)+DY(i))^2)-20)*((x(i)+DX(i))^2+(y(i)+DY(i))^2)*(x(i)+DX(i))*(y(i)+DY(i));
Z26(i)=sqrt(14)*(6*((x(i)+DX(i))^2+(y(i)+DY(i))^2)-5)*((x(i)+DX(i))^4-8*(x(i)+DX(i))^2*(y(i)+DY(i))^2+(y(i)+DY(i))^4);
Z27(i)=sqrt(14)*(6*(x(i)+DX(i))^4-24*(x(i)+DX(i))^2*(y(i)+DY(i))^2+6*(y(i)+DY(i))^4)*(x(i)+DX(i))*(y(i)+DY(i));
Z28(i)=sqrt(14)*((x(i)+DX(i))^6-17*(x(i)+DX(i))^4*(y(i)+DY(i))^2+17*(x(i)+DX(i))^2*(y(i)+DY(i))^4-(y(i)+DY(i))^6);
Z29(i)=-20*(x(i)+DX(i)).^3.*(y(i)+DY(i))+20*(x(i)+DX(i)).*(y(i)+DY(i)).^3+24*(x(i)+DX(i)).^3.*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)-24*(x(i)+DX(i)).*(y(i)+DY(i)).^3.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2);
Z30(i)=10*(x(i)+DX(i)).^3-30*(x(i)+DX(i)).*(y(i)+DY(i)).^2-30*(x(i)+DX(i)).^3.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+90*(x(i)+DX(i)).*(y(i)+DY(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+21*(x(i)+DX(i)).^3.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2-63*(x(i)+DX(i)).*(y(i)+DY(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2;
Z31(i)=30*(x(i)+DX(i)).^2.*(y(i)+DY(i))-10*(y(i)+DY(i)).^3-90*(x(i)+DX(i)).^2.*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+30*(y(i)+DY(i)).^3.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+63*(x(i)+DX(i)).^2.*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2-21*(y(i)+DY(i)).^3.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2;
Z32(i)=-10*(x(i)+DX(i)).^2+10*(y(i)+DY(i)).^2+60*(x(i)+DX(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)-60*(y(i)+DY(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)-105*(x(i)+DX(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2+105*(y(i)+DY(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2+56*(x(i)+DX(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3-56*(y(i)+DY(i)).^2.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3;
Z33(i)=-20*(x(i)+DX(i)).*(y(i)+DY(i))+120*(x(i)+DX(i)).*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)-210*(x(i)+DX(i)).*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2+112*(x(i)+DX(i)).*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3;
Z34(i)=5*(x(i)+DX(i))-60*(x(i)+DX(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+210*(x(i)+DX(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2-280*(x(i)+DX(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3+126*(x(i)+DX(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^4;
Z35(i)=5*(y(i)+DY(i))-60*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)+210*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2-280*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3+126*(y(i)+DY(i)).*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^4;
Z36(i)=-1+30*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2)-210*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^2+560.*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^3-630*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^4+252*((x(i)+DX(i)).^2+(y(i)+DY(i)).^2).^5;
zzz=[Z1 Z2 Z3 Z4 Z5 Z6 Z7 Z8 Z9 Z10 Z11 Z12 Z13 Z14 Z15 Z16 Z17 Z18 Z19 Z20 Z21 Z22 Z23 Z24 Z25 Z26 Z27 Z28 Z29 Z30 Z31 Z32 Z33 Z34 Z35 Z36];
end
A=reshape(zzz,[533,36]);%workbench导出的数据节点一共有533个
[k,res]=qrls(A,W);



程序二：调用程序qrls

function [k,res]=qrls(A,b)
[Q,R]=qr(A);f=Q'*b;
k=R\f; res=norm(b-A*k);


求得的36项系数：

6049659667.54738
7844884566.74876
216762438.556726
10365554433.7699
1871031.02626269
14658506300.6731
17582.7919732397
-7782607390.75034
-3749945.53506522
-23358639536.5575
2285591.27181637
4865182595.10850
0
15715532381.5909
-12453299292.5750
-6235411170.63858
-42288675.1348851
-8139305890.99716
-60808217.2942579
46463032.5780004
156832138.352915
-4530360121.35739
618375927.349488
4691493849.88505
988932936.019116
2124450161.35459
-2496678077.98561
3411993916.61092
-7880863374.96611
-4671269106.14148
-83169252.5029032
-1024932189.16105
1388708787.87335
1017089783.55429
1210281.41041704
87836356.2172928

求虫友帮忙看看问题，另外有用正交法做过光机热集成分析的朋友还望分享下方法！！！十分感谢！！
回复此楼

» 猜你喜欢
1楼2019-06-21 22:22:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keaSSy

至尊木虫 (著名写手)
你能不能上传一个TXT或者m文件，这个看起来太费劲了

发自小木虫IOS客户端
赞一下 回复此楼
你必须非常努力才能看起来毫不费力
2楼2019-06-21 22:33:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

keaSSy

至尊木虫 (著名写手)
首先，W有问题，而且你这个程序其实最好可以先做归一化处理，然后还有就是其实你这样拟合出来的应该是没有去处刚体位移的，也不是不行，但是应该先做径向修正，要不然求出来的不是矢高，存在的问题可能比较多
赞一下 回复此楼
你必须非常努力才能看起来毫不费力
3楼2019-06-21 22:50:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

哲理的山路

新虫 (初入文坛)
老铁现在会了不，求教啊～

发自小木虫Android客户端
赞一下 回复此楼
4楼2019-09-24 18:19:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LYC112233

新虫 (初入文坛)

【答案】应助回帖

可以请教一下楼主关于zernike拟合的问题吗？
赞一下 回复此楼
5楼2020-06-29 20:17:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 老铁双击999 的主题更新
51/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E05青基有几个评审 +4 KYXY123 2024-05-28 4/200 2024-05-28 19:25 by popt2t
[基金申请] E10开始送了，希望有好运 +5 sail 2024-05-27 5/250 2024-05-28 18:36 by 芝小芝
[有机交流] 液相纯度高，但产品析不太出来 10+4 cui19236 2024-05-27 6/300 2024-05-28 17:09 by chu-yue
[基金申请] 山东省基金：作为项目负责人承担面上项目数量累计不得超过2项 +13 未岸亡灵 2024-05-23 26/1300 2024-05-28 17:01 by xuechenli
[教师之家] 在大地上我们只过一生---看完我的阿勒泰上头了好几天，完结那天晚上几乎失眠 +4 瞬息宇宙 2024-05-27 4/200 2024-05-28 15:55 by chemhua
[论文投稿] EI学报，一审返修后，为啥不再送审，直接终审中? 10+3 qweasd12345 2024-05-27 5/250 2024-05-28 14:55 by topedit
[论文投稿] 真急着毕业，CPB主编终审17天了，邮件催稿了两次，就是一点动静没有 5+3 kkkk夏 2024-05-28 5/250 2024-05-28 14:43 by topedit
[教师之家] 被惯着的学生终究要吃大亏 +27 535743368 2024-05-24 29/1450 2024-05-28 11:44 by Pingfan010
[有机交流] D-阿拉伯糖-1，4-内酯的合成 10+4 Leeu55 2024-05-24 10/500 2024-05-28 11:34 by disguiseme
[考博] 24年博士招生 +4 abinit432 2024-05-27 5/250 2024-05-28 10:47 by 牛奶小面包
[电化学] 2200mA锂离子电池设计 +5 小乐987 2024-05-23 5/250 2024-05-28 10:02 by 打工艺术家
[考博] 24或25申博 +4 Jacob- 2024-05-22 7/350 2024-05-27 21:21 by 安塔瓦拉多
[教师之家] 女博士高校择业三天之内签合同，求支招 +38 chengmy19 2024-05-23 55/2750 2024-05-27 18:02 by 好男人￥
[基金申请] 科研之友阅读量近一周增加了200多。 +13 hdzw9071 2024-05-24 14/700 2024-05-27 17:29 by leeeru
[教师之家] 研究生在毕业答辩时挂了，遗憾 +7 spacr2022 2024-05-27 10/500 2024-05-27 12:45 by 懒羊羊2024
[硕博家园] 周日 +6 1加油哦棒 2024-05-26 9/450 2024-05-27 10:30 by hahamyid
[论文投稿] 共同通讯谁投稿？ +3 qvhm2609 2024-05-23 10/500 2024-05-26 13:25 by bobvan
[硕博家园] 人生 +15 暮色恋伊人 2024-05-22 15/750 2024-05-26 08:23 by elainzai
[论文投稿] materials letter +4 烟雨盛世 2024-05-24 4/200 2024-05-25 02:07 by 1872075
[有机交流] 苯磺酰氯与醇羟基反应 5+4 杨怼怼? 2024-05-22 12/600 2024-05-23 15:59 by mrzhl1986
信息提示
关闭
请填处理意见
关闭