作业帮 > 综合 > 作业

matlab程序 解方程组

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/05/22 08:01:06
matlab程序 解方程组
syms a b Rsd Rsp
f1=72.8*(1+cos(a))-2*(sqrt(51*Rsp)+sqrt(21.8*Rsd));
f2=50.8*(1+cos(b))-2*(sqrt(2.3*Rsp)+sqrt(48.5*Rsd));
A=input('输入:A=');
B=input('输入:B=');
for i=1:length(A)
a=A(i);
b=B(i);
f11=subs(f1);
f22=subs(f2);
[Rsp1,Rsd1]=solve(f11,f22);
Rsp0=double(Rsp1);
Rsd0=double(Rsd1);
disp(['Rsp' num2str(i) ' = '])
disp(Rsp0)
disp(['Rsd' num2str(i) ' = '])
disp(Rsd0)
disp(['R' num2str(i) ' = '])
disp(Rsp0+Rsd0)
end
要求 1.输出格式尽量整齐,有warning 尽量帮我屏蔽,最好是能把
R
Rsp
Rsd
做的可以整齐的导入excel,直接在程序中写了更好!
2.
如果您对解方程有更好的想法请您说一下,或者这么是否有问题!
3.大致可以就行,我给200分,给100,成功加100
syms a b Rsd Rsp
Ral=72.8;
Ralp=51;
Rald=21.8;
Rbl=50.8;
Rblp=2.3;
Rbld=48.5;
f1=Ral*(1+cos(a))-2*(sqrt(Ralp*Rsp)+sqrt(Rald*Rsd));
f2=Rbl*(1+cos(b))-2*(sqrt(Rblp*Rsp)+sqrt(Rbld*Rsd));
A=input('输入:A=');
B=input('输入:B=');
n=length(A);
base=cell(3,n);
warning off
for i=1:n
a=A(i);
b=B(i);
f11=subs(f1);
f22=subs(f2);
[Rsp1,Rsd1]=solve(f11,f22);
Rsp0=double(Rsp1);
Rsd0=double(Rsd1);
R0=Rsp0+Rsd0;
base{i,1}=[Rsp0;0];
base{i,2}=[Rsd0;0];
base{i,3}=[R0;0];
% disp(['Rsp' num2str(i) ' = '])
% disp(Rsp0)
% disp(['Rsd' num2str(i) ' = '])
% disp(Rsd0)
% disp(['R' num2str(i) ' = '])
% disp(Rsp0+Rsd0)
end
base1=cell2mat(base);
disp(base1)
mydir=uigetdir('d:','选择一个excel保存目录');%%%选择文件的保存位置
%%%可以选择已存在文件夹,也可以新建文件夹,也可以
%%%选择D盘以外的文件夹,这里的d只起参考作用
name='jie';%%%此处改名字
xlswrite([mydir '\' name '.xls'],base1);%数据存储至D盘根目录下
warning on