How to plot a function for different values of input parameter in a function

1 次查看（过去 30 天）
University Glasgow
University Glasgow 2022-9-13
评论： University Glasgow 2022-9-13
Hi,
Please, how do I create a procedure/function such that it carries \xi as an input and use it to plot r for different values of \xi. I have attached my code.

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2022-9-13
Just define xi as a column vector of values, e.g.
xi = [0.03 0.06].'; % activity parameter
  11 个评论
显示 9更早的评论
Torsten
Torsten 2022-9-13
Ran in:
%% initialization
clear
clc
close all
%% Model parameters
global k1 eta1 alpha3 gamma1 d N Phi xi h A B C G
k1 = 6*10^(-12); % elastic constant
eta1 = 0.0240; % viscosity
xi = [-0.02 0.01 0.0].'; % activity parameter
alpha3 = -0.001104; % viscosity
gamma1 = 0.1093; % viscosity
Theta = 0.0001;
d = 0.0002;
%% Extract r(q) data from maple
%rqdata= xlsread('rqdata.xlsx');
%rq= rqdata(:,2);
%qsize = rqdata(:,1);
%% r(q) from matlab
alpha=1-alpha3^2/gamma1/eta1;
q=0:0.001:(5*pi);
r=q-(1-alpha)*tan(q)+(alpha3*xi*alpha/eta1)./(4*k1*q.^2/d^2-alpha3*xi/eta1).*tan(q);
% Plotting r(q)
%figure
%hold on
%plot(qsize,rq, 'r', MarkerSize=2)
plot(q,r, 'b--', MarkerSize=2)
xlabel('t (seconds)')
ylabel('\theta(d/2,t)(rad)')
legend('maple', 'matlab')
axis([0 5*pi -20 20])
University Glasgow
University Glasgow 2022-9-13
Thank you Torsten. This is actually working fine. But I have a different problem now. This is just part of my codes. I didn't wanna post all code. I thought using this as an instance could solve my problem.
In the attached file,I also wanna plot theta_sol and v_sol for xi = [0.03 0.06].'; but I received an error: Error in Case1_Ijuptilk_130822>lcode1 (line 147): rhsode(N,1) = (G/(h^3))*(-Phi + 3*theta(1) - 3*theta(2) + theta(3)) + (C/(h^2))*(v(2) -2*v(1)+ 0)
Please could you me check it out.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Annotations 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by