Matlab求数值解
发布网友
发布时间：2022-04-29 17:53
共1个回答
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
热心网友
时间：2023-10-26 12:13
第一步，将已知值分别赋值给变量，即
A_Acr0=0.1:0.01:0.25;
k0=1.2:0.1:1.4;
第二步，使用solve函数求解Ma值。即
Ma=solve(A_Acr==(((1+0.5*(k-1)*Ma^2)*2/(k+1))^(0.5*(k+1)/(k-1)))/Ma;
第三步，使用for双循环语句，得到所有的Ma值。即
for i=1:length(A_Acr0)
for j=1:length(k0)
sysm Ma
Ma=solve(eq==A_Acr)
Ma0(i,j)=Ma(1);
end
end
第四步，显示A/Acr、 k、 Ma值，即
disp(A/Acr)
disp(k)
disp(Ma)追问
以下是我的目的
Beta = 1 - 1./A_Acr
Beta - Ma 曲线图：
所以我需要解出来的Ma值除了数值解之外，还需要是在区间[0, 1]的实数解，你的方法我之前试过了，一般情况下会得到多个root解，vpa后，结果还是带有复数的。
