| 查看: 663 | 回复: 0 | ||
| 【悬赏金币】回答本帖问题，作者xmqiu将赠送您 5 个金币 | ||
xmqiu金虫 (小有名气)
|
[求助]
matlab for循环 解方程
|
|
|
我想求解一个方程，方程四元一次，系数有7个数，每次采用4个来解方程，x1,x2,x3,x4用完，然后保存求解的a b c d的值，接着跳到x2,x3,x4,x5有解出一组a b c d的值......可是我写的程序总是跳出错误使用 sym/vpasolve (line 145) Unable to find variables in equations. 方程式求解一个四元一次方程，每次用四个x的值解出一组a b c d。一组x1 x2 x3 x4解出一组，然后跳到x2 x3 x4 x5,再解出一组.....求教啊！！在线等！！ syms a b c d %e=[];f=[]; vr1list=[];vr2list=[];vr3list=[];vr4list=[]; for i=1:length(e) x(i)=e(i); % A[k+1]=B[k] end for j=1:length(f) y(j)=f(j); % A[k+1]=B[k] end for k=1:length(e) [a,b,c,d]=vpasolve(a+b*x(k)*exp(-c*x(k))+d*x(k)-y(k)==0,a+b*x(k+1)*exp(-c*x(k+1))+d*x(k+1)-y(k+1)==0,a+b*x(k+2)*exp(-c*x(k+2))+d*x(k+2)-y(k+2)==0,a+b*x(7)*exp(-c*x(k+3))+d*x(k+3)-y(k+3)==0); A=double([a b c d]) vr1=A(:,1);vr2=A(:,2);vr3=A(:,3);vr4=A(:,4); vr1list=[vr1list vr1]; vr2list=[vr2list vr2]; vr3list=[vr3list vr3]; vr4list=[vr4list vr4]; end |
回复此楼» 猜你喜欢
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有202人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
-
求助NH4V4O10晶体的CIF文件
已经有0人回复
-
英国全奖博士招聘-深度学习与量子物理
已经有0人回复
-
间接带隙半导体有效质量求助
已经有0人回复