| 查看: 1248 | 回复: 14 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
飞扬2282荣誉版主 (著名写手)
|
[交流]
【求助】Matlab最小二乘法简单编程代码？
|
||
|
求Matlab7编写程序代码，实现如下简单计算！ (uj-ui)v-(vj-vi)u=vi*uj-vj*ui ( i=1,2,....N/2; j=i+N/2 ) 其中： ui，uj，vi，vj，N均为已知数值，求u,v? [ Last edited by 飞扬2282 on 2009-4-8 at 15:43 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
上海工程技术大学 激光智能制造课题组 2027级博士研究生招生
已经有6人回复
-
课题组招2027级博士 上海工程技术大学 激光智能制造方向
已经有6人回复
-
现代”学阀”该如何界定
已经有14人回复
-
求合成方法
已经有8人回复
-
师弟论文见刊大半年才想起来申请专利，还能抢救一下吗？
已经有4人回复
-
有机合成以后会不会被AI改变？做科研的虫友怎么看
已经有8人回复
-
我的奶奶
已经有3人回复
★ ★ ★ ★
sunxiao(金币+4,VIP+0):谢谢参与，欢迎常来仿真编程版 4-9 23:13
sunxiao(金币+4,VIP+0):谢谢参与，欢迎常来仿真编程版 4-9 23:13
|
我老是出点错，麻烦版主退10个金币给楼主吧，前面几个人也费了不少劲，呵呵 我把程序改了一下 function y=solveuv(u,v,N) %Av-Bu=C ind=floor(N/2); A=u(ind+1:ind*2)-u(1:ind); B=v(ind+1:ind*2)-v(1:ind); C=v(1:ind).*u(ind+1:ind*2)-v(ind+1:ind*2).*u(1:ind); %[A -B]*[v u]'=C y=([A -B])\C; % y(1)=u % y(2)=v y=[y(2) y(1)]'; >> N=6; u=[1 3 3 4 5 6]'; v=[1.5 2.5 3.5 4.5 5.5 6.5]'; y=solveuv(u,v,N) y = 5.0000 5.5000 输出第一个值是u，第二个值是v，我只是根据你的式子来看的，没有考虑实际意思，感觉应该符合你的方程了。 第1个问题： ind+1:ind不就是ind+1:ind+ind,也就是ind+1:ind*2吗，一样的。 第2个问题： *是矩阵运算，.*是数组运算 比如： >> a=[1 2 ;3 4] a = 1 2 3 4 >> a.*a ans = 1 4 9 16 >> a*a ans = 7 10 15 22 随便找本MATLAB入门级的书先讲的都是这个 |
11楼2009-04-09 18:28:54
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来仿真编程版 4-8 23:44
飞扬2282(金币+8,VIP+0):如何执行命令求出u，v的值呢？函数执行不了啊 4-9 10:15
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来仿真编程版 4-8 23:44
飞扬2282(金币+8,VIP+0):如何执行命令求出u，v的值呢？函数执行不了啊 4-9 10:15
|
function y=solveuv(u,v,N) %Av-Bu=C ind=floor(N/2); A=u(ind+1:ind*2)-u(1:ind); B=v(ind+1:ind*2)-v(1:ind); C=v(1:ind).*u(ind+1:ind*2)-v(ind+1:ind*2).*u(1:ind); %[A B]*[u v]'=C y=([A B])\C; |
2楼2009-04-08 18:46:42
飞扬2282
荣誉版主 (著名写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.315
- 金币: 10773.4
- 散金: 1185
- 红花: 17
- 帖子: 1897
- 在线: 709.6小时
- 虫号: 661639
- 注册: 2008-11-25
- 性别: GG
- 专业: 核技术及其应用
- 管辖: 学术会议
3楼2009-04-09 10:15:23
4楼2009-04-09 10:48:44