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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 光机热集成分析时采用Zernike多项式对变形镜面拟合的问题！（最小二乘法）
51/1返回列表
查看: 1222  |  回复: 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美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 函评这一关首先干掉一大半，今年卷上天了 +10 mycaas 2024-05-07 13/650 2024-05-07 10:45 by pizi1024
[教师之家] 双非本科，学生又懒又笨。双一流本科好一些？入职门槛和考核要求比双非本科高多少？ +4 河西夜郎 2024-05-05 4/200 2024-05-07 10:22 by Quakerbird
[有机交流] DMF怎么除水 +10 爱吃香菜的小黄� 2024-05-01 15/750 2024-05-07 09:49 by zyp0009928
[论文投稿] 审稿人一收一拒，编辑给了修 +10 njuyck 2024-04-30 23/1150 2024-05-07 09:27 by BusyGer
[基金申请] 特聘教授，算高级职称吗？青年基金还没结题的情况下还能报面上吗？ +12 aasahr 2024-05-05 14/700 2024-05-07 09:22 by vmuch520
[硕博家园] 大家能不能给点意见啊 +12 是小蚂蚁啦 2024-05-01 23/1150 2024-05-07 09:07 by 是小蚂蚁啦
[基金申请] AI派发是根据关键词、基金题目、研究者研究方向匹配？ 200+3 flllgend 2024-05-04 10/500 2024-05-06 23:53 by dxcharlary
[论文投稿] 想请问下IEEE Trans最后一轮，你们等了多久？ +4 DD4Dr 2024-05-03 9/450 2024-05-06 23:19 by topedit
[教师之家] 领导好安排自己做杂事。在成果上给领导挂名，请其不要经常安排自己做杂事 +7 河西夜郎 2024-05-03 7/350 2024-05-06 17:56 by lyfbangong
[硕博家园] 去美国的华人导师底下读博后会很难出站吗？ +7 Aries-FZU 2024-05-03 7/350 2024-05-06 17:46 by jmmchen
[考博] 2024还有没有补录的博士 +7 crazy peng 2024-04-30 8/400 2024-05-06 15:37 by 愿你快乐啊
[有机交流] 酯化和酯交换 40+4 whw818 2024-05-04 5/250 2024-05-06 14:24 by 世界之源
[考博] 25年春秋博士申请 +4 梦想研究生 2024-05-04 7/350 2024-05-06 11:24 by 梦想研究生
[催化] 孔径分布怎么分析？ +3 林佳树 2024-05-02 7/350 2024-05-05 14:38 by 醉舞清风
[论文投稿] 想问下电池发在nano letters含金量高不高 +6 未来富婆蛙 2024-04-30 11/550 2024-05-05 10:53 by eemdl
[考博] 西安交通大学转化医学研究院蔡文课题组招收2025年入学的申请-考核制博士生 +3 caiwen2016 2024-05-03 3/150 2024-05-05 09:39 by 宋mlk
[基金申请] 关于公布2024年度国家自然科学基金项目申请初审结果的通告 (金币+2) +12 hdzw9071 2024-04-30 18/900 2024-05-04 17:49 by felicity6056
[考博] 北京师范大学人工智能学院2024年第二批博士研究生招生报名 +4 taoto 2024-05-02 5/250 2024-05-04 15:46 by LYS1200
[论文投稿] 第一篇论文终于进去审稿状态了 +5 lizhengke06 2024-05-01 8/400 2024-05-03 09:14 by nono2009
[高分子] 高分子聚合物有机材料方向的合作/招聘科研助理 +3 wwwzg 2024-05-01 3/150 2024-05-02 21:38 by Dr. Zhang
信息提示
关闭
请填处理意见
关闭